From: Raymond T. <rt...@us...> - 2004-11-29 18:33:48
|
Update of /cvsroot/maxima/maxima/tests In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv22473/tests Modified Files: rtesthyp.mac Log Message: Correct and update some of the expected results to match what maxima actually produces. (Mostly because maxima produces 4^(foo/2) instead of 2^foo in some cases.) Added some calls to expand to get slightly simpler results and to expand some of the Bessel functions into elementary functions. Index: rtesthyp.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtesthyp.mac,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- rtesthyp.mac 26 Nov 2004 15:57:56 -0000 1.1 +++ rtesthyp.mac 29 Nov 2004 18:33:32 -0000 1.2 @@ -8,23 +8,23 @@ /* A&S 13.6.1 */ hgfred([v+1/2],[2*v+1],2*%i*z); -gamma(1+v)*exp(%i*z)*(z/2)^(-v)*bessel_j(v,z)$ +4^(v/2)*bessel_j(v,z)*gamma(v+1)*exp(%i*z)/z^v$ /* A&S 13.6.3 */ hgfred([v+1/2],[2*v+1],2*z); -gamma(1+v)*exp(z)*(z/2)^(-v)*bessel_i(v,z)$ +4^(v/2)*bessel_i(v,z)*gamma(v+1)*exp(z)/z^v$ /* A&S 13.6.4 */ -hgfred([n+1],[2*n+2],2*%i*z); -gamma(n+3/2)*exp(%i*z)*(z/2)^(-n-1/2)*bessel_j(n+1/2,z)$ +expand(hgfred([n+1],[2*n+2],2*%i*z)); +4^(n/2+1/4)*bessel_j(n+1/2,z)*gamma(n+3/2)*z^(-n-1/2)*%e^(%i*z)$ /* A&S 13.6.5 */ -hgfred([-n],[-2*n],2*%i*z); -gamma(1/2-n)*exp(%i*z)*(z/2)^(n+1/2)*bessel_j(-n-1/2,z)$ +expand(hgfred([-n],[-2*n],2*%i*z)); +bessel_j(-n-1/2,z)*gamma(1/2-n)*4^(-n/2-1/4)*z^(n+1/2)*%e^(%i*z)$ /* A&S 13.6.6 */ hgfred([n+1],[2*n+2],2*z); -gamma(3/2+n)*exp(z)*(z/2)^(-n-1/2)*bessel_i(n+1/2,z)$ +bessel_i(n+1/2,z)*gamma(n+3/2)*z^(-n-1/2)*%e^z/4^((-n-1/2)/2)$ /* A&S 13.6.10 */ hgfred([a],[a+1],-z); @@ -35,11 +35,14 @@ exp(z)$ /* A&S 13.6.13 */ -hgfred([1],[2],-2*%i*z); +besselexpand:true; +true$ + +expand(hgfred([1],[2],-2*%i*z)); exp(-%i*z)*sin(z)/z$ /* A&S 13.6.14 */ -hgfred([1],[2],2*z); +expand(hgfred([1],[2],2*z)); exp(z)*sinh(z)/z$ /* A&S 13.6.17 */ @@ -54,6 +57,11 @@ hgfred([1/2],[3/2],-z^2); sqrt(%pi)/2/z*erf(z); +/* Hypergeometric function 2F1 tests +* +* See A&S 15.1 +* +*/ /* A&S 15.1.3 */ hgfred([1,1],[2],z); -log(1-z)/z; |