[Maxima-commits] CVS: maxima/tests rtest14.mac,1.66,1.67

[Dieter K. <cra...@us...>]

Update of /cvsroot/maxima/maxima/tests
In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv22233/tests

Modified Files:
	rtest14.mac 
Log Message:
Adding examples to show that the functions hankel_1 and hankel_2 now works too.


Index: rtest14.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/rtest14.mac,v
retrieving revision 1.66
retrieving revision 1.67
diff -u -d -r1.66 -r1.67
--- rtest14.mac	26 Sep 2009 12:26:46 -0000	1.66
+++ rtest14.mac	26 Sep 2009 17:40:36 -0000	1.67
@@ -459,12 +459,18 @@
  *
  */
 
-specint( t*%h[2/3,1](t^(1/2))*%e^(-p*t),t);
+specint(t*%h[2/3,1](t^(1/2))*%e^(-p*t),t);
 -4*%i*gamma(1/3)*%m[-3/2,1/3](-1/(4*p))*%e^-(1/(8*p))/(3*(-1)^(5/6)*sqrt(3) 
  *gamma(2/3)*p^(3/2))+4*gamma(1/3)*%m[-3/2,1/3](-1/(4*p))*%e^-(1/(8*p))/(3*( 
  -1)^(5/6)*gamma(2/3)*p^(3/2))-8*%i*gamma(2/3)*%m[-3/2,-1/3](-1/(4*p))*%e^- 
  (1/(8*p))/(3*(-1)^(1/6)*sqrt(3)*gamma(1/3)*p^(3/2)) $
 
+specint(t*hankel_1(2/3,t^(1/2))*%e^(-p*t),t);
+-4*%i*gamma(1/3)*%m[-3/2,1/3](-1/(4*p))*%e^-(1/(8*p))/(3*(-1)^(5/6)*sqrt(3) 
+ *gamma(2/3)*p^(3/2))+4*gamma(1/3)*%m[-3/2,1/3](-1/(4*p))*%e^-(1/(8*p))/(3*( 
+ -1)^(5/6)*gamma(2/3)*p^(3/2))-8*%i*gamma(2/3)*%m[-3/2,-1/3](-1/(4*p))*%e^- 
+ (1/(8*p))/(3*(-1)^(1/6)*sqrt(3)*gamma(1/3)*p^(3/2)) $
+ 
 /*
  * H[3/4,2](t) = bessel_j(3/4,t)-%i*bessel_y(3/4,t)
  *
@@ -492,13 +498,20 @@
  *    =  r^(-3/2)*(-gamma(9/4)*assoc_legendre_p(1/2,-3/4,p/r)
  *                  -gamma(3/4)*sqrt(2)*assoc_legendre_p(1/2,3/4,p/r))
  */
-ratsimp(specint( t^(1/2)*%h[3/4,2](t)*%e^(-p*t),t));
+ratsimp(specint(t^(1/2)*%h[3/4,2](t)*%e^(-p*t),t));
 ''(ratsimp(5*gamma(1/4)/16/(p^2+1)^(3/4)*assoc_legendre_p(1/2,-3/4,p/sqrt(p^2+1))
 +sqrt(2)*%i*assoc_legendre_p(1/2,3/4,p/sqrt(p^2+1))*gamma(3/4)
        /(p^2+1)^(3/4)
        +5*%i*gamma(1/4)*assoc_legendre_p(1/2,-3/4,p/sqrt(p^2+1))
 	/(16*(p^2+1)^(3/4))))$
 
+ratsimp(specint(t^(1/2)*hankel_2(3/4,t)*%e^(-p*t),t));
+''(ratsimp(5*gamma(1/4)/16/(p^2+1)^(3/4)*assoc_legendre_p(1/2,-3/4,p/sqrt(p^2+1))
++sqrt(2)*%i*assoc_legendre_p(1/2,3/4,p/sqrt(p^2+1))*gamma(3/4)
+       /(p^2+1)^(3/4)
+       +5*%i*gamma(1/4)*assoc_legendre_p(1/2,-3/4,p/sqrt(p^2+1))
+   /(16*(p^2+1)^(3/4))))$
+
 /*
  * %h[1/2,1](t) = bessel_j(1/2,t)+%i*bessel_y(1/2,t)
  *
@@ -533,7 +546,10 @@
  *    -> -sqrt(2)*(p-1)*(p+1)/(sqrt(%pi)*(p^2+1)^2)
  */
 
-ratsimp(specint( t^(3/2)*%h[1/2,1](t)*%e^(-p*t),t));
+ratsimp(specint(t^(3/2)*%h[1/2,1](t)*%e^(-p*t),t));
+-((sqrt(%pi)*(sqrt(2)*%i*p^2-2*sqrt(2)*p-sqrt(2)*%i))/(%pi*p^4+2*%pi*p^2+%pi))$
+
+ratsimp(specint(t^(3/2)*hankel_1(1/2,t)*%e^(-p*t),t));
 -((sqrt(%pi)*(sqrt(2)*%i*p^2-2*sqrt(2)*p-sqrt(2)*%i))/(%pi*p^4+2*%pi*p^2+%pi))$
 
 /*