[Maxima-commits] CVS: maxima/src hypgeo.lisp,1.10,1.11

[Raymond T. <rt...@us...>]

Update of /cvsroot/maxima/maxima/src
In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv30148/src

Modified Files:
	hypgeo.lisp 
Log Message:
Add missig function implementing formula 50, page 188 from Table of
Integral Transforms.


Index: hypgeo.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/src/hypgeo.lisp,v
retrieving revision 1.10
retrieving revision 1.11
diff -u -d -r1.10 -r1.11
--- hypgeo.lisp	9 Feb 2005 04:21:51 -0000	1.10
+++ hypgeo.lisp	19 Feb 2005 21:17:00 -0000	1.11
@@ -2639,6 +2639,26 @@
 	(return (f50p188-simp d m v a))))
      (return 'fail-in-f50cond))) 
 
+;; Table of Integral Transforms
+;;
+;; p. 188, formula 50:
+;;
+;; t^(u-1/2)*bessel_y(2*v,2*sqrt(a)*sqrt(t))
+;;    -> a^(-1/2)*p^(-u)*exp(-a/2/p)
+;;       * [tan((u-v)*%pi)*gamma(u+v+1/2)/gamma(2*v+1)*M[u,v](a/p)
+;;          -sec((u-v)*%pi)*W[u,v](a/p)]
+(defun f50p188-simp (d u v a)
+  (mul d
+       (power a (inv -2))
+       (power par (mul -1 u))
+       (power '$%e (div a (mul -2 par)))
+       (sub (mul (tan% (mul '$%pi (sub u v)))
+		 (gm (add u v (inv 2)))
+		 (inv (gm (add v v 1)))
+		 (mwhit (div a par) u v))
+	    (mul `((%sec) ,(mul '$%pi (sub u v)))
+		 (wwhit (div a par) u v)))))
+
 (defun f2p105v2cond-simp (m v a) 
   (mul -2.
        (power '$%pi -1.)