FUNCTION : mixedajc2jac - converts a sum of quotients of jacprods to a common base
CALLING SEQUENCE : mixedjac2jac(jacexp,T)
PARAMETERS : jacexp - sum of quotients of theta functions written
in terms of JAC(a,b,infinity)
T - positive integer
SYNOPSIS : jac2series and jacprodmake are used to convert each quotient into a
a quotient with the same base. T is chosen large enough for conversion to work.
NOTE: Must contain at least two terms.
EXAMPLES :
> J1:=1-1/JAC(0,10,infinity)^2/JAC(0,120,infinity)^3/JAC(1,3,infinity)/JAC(2,6,infinity)
/JAC(3,10,infinity)/JAC(4,12,infinity)/JAC(8,20,infinity)^2*JAC(0,1,infinity)
*JAC(0,2,infinity)*JAC(0,5,infinity)^2*JAC(0,8,infinity)*JAC(4,10,infinity)
*JAC(6,20,infinity)*JAC(16,40,infinity)*JAC(120,360,infinity)^3;
> J2 := mixedjac2jac(J1,360);
%term %, 1, %of %, 2
%term %, 2, %of %, 2
2 2 2
J2 := 1 - JAC(5, 40, infinity) JAC(6, 40, infinity) JAC(14, 40, infinity)
2 2 /
JAC(15, 40, infinity) JAC(16, 40, infinity) / (JAC(3, 40, infinity)
/
2 3
JAC(7, 40, infinity) JAC(8, 40, infinity) JAC(12, 40, infinity)
JAC(13, 40, infinity) JAC(17, 40, infinity) JAC(20, 40, infinity))
> series(jac2series(J1-J2,1000),q,1000);
1000
O(q )
DISCUSSION : We able to convert given jacexp with mixed bases (divisors of 360)
to base 40.
SEE ALSO :