From: Dieter K. <cra...@us...> - 2009-09-23 21:53:23
|
Update of /cvsroot/maxima/maxima/tests In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv16696/tests Modified Files: rtest_hypgeo.mac Log Message: Adding a test for the Laplace transform of the square of bessel_i(1/2,t). Related bug report: Bug ID: 2864588 - specint(exp(-s*t)*bessel_i(1/2,t)^2,t) not correct Index: rtest_hypgeo.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtest_hypgeo.mac,v retrieving revision 1.10 retrieving revision 1.11 diff -u -d -r1.10 -r1.11 --- rtest_hypgeo.mac 13 Sep 2009 19:52:30 -0000 1.10 +++ rtest_hypgeo.mac 23 Sep 2009 21:53:11 -0000 1.11 @@ -1908,9 +1908,8 @@ radcan(specint(%e^(-s*t)*bessel_k(0,a*t),t)); specint(%e^(-s*t)*bessel_k(0,a*t),t); -/******************************************************************************/ - -/* +/******************************************************************************* + * * Bug ID: 2858045 - Laplace transform of asin or atan gives wrong noun form * Check the correct noun for the asin and atan function. */ @@ -1929,6 +1928,20 @@ specint(%e^(-s*t)*atan(1/sqrt(t)),t); sqrt(%pi)*%f[3,1]([1/2,1,-1/2],[3/2],s)/sqrt(s); +/******************************************************************************* + * + * Bug ID: 2864588 - specint(exp(-s*t)*bessel_i(1/2,t)^2,t) not correct + * + * The Laplace transform of the square of the bessel_i(1/2,t) function is equal + * 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; + -log(1-4/s^2)/(2*%pi); + +res2:specint(exp(-s*t)*bessel_i(1/2,t)^2,t),besselexpand:true; + -log(1-4/s^2)/(2*%pi); + /******************************************************************************/ kill(all); |