From: Raymond T. <rt...@us...> - 2005-02-27 02:53:31
|
Update of /cvsroot/maxima/maxima/tests In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv27409/tests Modified Files: rtest14.mac rtesthyp.mac Log Message: src/hyp.lisp: o Implement branch cuts for LEGF20 and LEGF24 tests/rtest14.mac o specint(t^u*bessel_j(v,a*t)*exp(-p*t),t) had the wrong expected result. We get this right now. tests/rtesthyp.lisp: o Add necessary assume's so we select the desired branch cut for the tests that use LEGF20 and LEGF24. Index: rtest14.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtest14.mac,v retrieving revision 1.18 retrieving revision 1.19 diff -u -d -r1.18 -r1.19 --- rtest14.mac 9 Feb 2005 22:01:37 -0000 1.18 +++ rtest14.mac 27 Feb 2005 02:53:21 -0000 1.19 @@ -381,7 +381,7 @@ radcan(specint(t^u*bessel_j(v,a*t)*exp(-p*t),t)); assoc_legendre_p(-u-1,-v,p/sqrt(p^2+a^2))*gamma(v+u+1) - /((p^2+a^2)^((u+1)/2)*(-1)^(v/2))$ + /((p^2+a^2)^((u+1)/2))$ /* * (25) Index: rtesthyp.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtesthyp.mac,v retrieving revision 1.24 retrieving revision 1.25 diff -u -d -r1.24 -r1.25 --- rtesthyp.mac 26 Feb 2005 20:33:07 -0000 1.24 +++ rtesthyp.mac 27 Feb 2005 02:53:21 -0000 1.25 @@ -211,23 +211,27 @@ * * Test for legf24 function. */ -radcan(hgfred([a,a+1/2],[c],z)*2^(1-c)*z^(c/2-1/2)*(1-z)^(1/2+a-c/2)/gamma(c)); +assume(zg1 > 1); +[zg1 > 1]; +radcan(hgfred([a,a+1/2],[c],zg1)*2^(1-c)*zg1^(c/2-1/2)*(1-zg1)^(1/2+a-c/2)/gamma(c)); /* Use A&S 8.2.1 to get this */ -assoc_legendre_p(c-2*a-1,1-c,1/sqrt(1-z))$ +assoc_legendre_p(c-2*a-1,1-c,1/sqrt(1-zg1))$ -2^u*(z^2-1)^(-u/2)*z^(v+u)/gamma(1-u)*hgfred([-v/2-u/2,1/2-v/2-u/2],[1-u],1-1/z^2); -assoc_legendre_p(v,u,z)$ +radcan(2^u*(zg1^2-1)^(-u/2)*zg1^(v+u)/gamma(1-u)*hgfred([-v/2-u/2,1/2-v/2-u/2],[1-u],1-1/zg1^2)); +assoc_legendre_p(v,u,zg1)$ /* A&S 15.4.12 * * F(a,b;a+b+1/2;z) = 2^(a+b-1/2)*gamma(a+b+1/2)*(-z)^(1/2-a-b)*P(a-b-1/2,1/2-a-b,sqrt(1-z)) * Test for legf20 function. */ -expand(ratsimp(2^(1/2-a-b)*(-x)^(1/2*(a+b-1/2))/gamma(a+b+1/2)*hgfred([a,b],[a+b+1/2],x))); -assoc_legendre_p(-b+a-1/2,1/2-a-b,sqrt(1-x))$ +assume(x2>1); +[x2 > 1]$ +expand(ratsimp(2^(1/2-a-b)*(-x2)^(1/2*(a+b-1/2))/gamma(a+b+1/2)*hgfred([a,b],[a+b+1/2],x2))); +assoc_legendre_p(-b+a-1/2,1/2-a-b,sqrt(1-x2))$ -2^m*(z^2-1)^(-m/2)/gamma(1-m)*hgfred([1/2+n/2-m/2, -n/2-m/2], [1-m], 1-z^2); -assoc_legendre_p(n,m,z)$ +2^m*(z1^2-1)^(-m/2)/gamma(1-m)*hgfred([1/2+n/2-m/2, -n/2-m/2], [1-m], 1-z1^2); +assoc_legendre_p(n,m,z1)$ /* * Formula 36: |