Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project!

## [Maxima-commits] CVS: maxima/tests rtest_hypgeo.mac,1.12,1.13

 [Maxima-commits] CVS: maxima/tests rtest_hypgeo.mac,1.12,1.13 From: Dieter Kaiser - 2009-09-26 14:53:12 Update of /cvsroot/maxima/maxima/tests In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv31296/tests Modified Files: rtest_hypgeo.mac Log Message: Adding two examples for the Laplace transform of the square of bessel_j(v,t) where v=-1/2 and v=1/2. Index: rtest_hypgeo.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtest_hypgeo.mac,v retrieving revision 1.12 retrieving revision 1.13 diff -u -d -r1.12 -r1.13 --- rtest_hypgeo.mac 25 Sep 2009 18:37:18 -0000 1.12 +++ rtest_hypgeo.mac 26 Sep 2009 12:25:11 -0000 1.13 @@ -1936,12 +1936,36 @@ * to the Laplace transform of the expanded form 2*sinh(t)^2/(%pi*t). */ -res1:specint(exp(-s*t)*bessel_i(1/2,t)^2,t),besselexpand:false; +specint(exp(-s*t)*bessel_i(1/2,t)^2,t),besselexpand:false; -log(1-4/s^2)/(2*%pi); -res2:specint(exp(-s*t)*bessel_i(1/2,t)^2,t),besselexpand:true; +specint(exp(-s*t)*bessel_i(1/2,t)^2,t),besselexpand:true; -log(1-4/s^2)/(2*%pi); +/******************************************************************************* + * + * Bug ID: 2866802 - specint(exp(-s*t)*t^(5/2)*bessel_j(-1/2,sqrt(t))^2,t) wrong + * + * We add two tests to show that the square of bessel_j(v,t) works for + * v=-1/2 (the orignally bug) and v=1/2 as expected. + */ + +res:factor(ratsimp(specint(exp(-s*t)*t*bessel_j(-1/2,t)^2,t))), + besselexpand:false; +2*(s^2+2)/(%pi*s*(s^2+4)); + +expand(res-ratsimp(specint(exp(-s*t)*t*bessel_j(-1/2,t)^2,t))), + besselexpand:true; +0; + +res:factor(ratsimp(specint(exp(-s*t)*t*bessel_j(1/2,t)^2,t))), + besselexpand:false; +4/(%pi*s*(s^2+4)); + +expand(res-ratsimp(specint(exp(-s*t)*t*bessel_j(1/2,t)^2,t))), + besselexpand:true; +0; + /******************************************************************************/ kill(all);