From: Dieter K. <cra...@us...> - 2009-10-12 21:26:08
|
Update of /cvsroot/maxima/maxima/tests In directory 23jxhf1.ch3.sourceforge.com:/tmp/cvs-serv20371/tests Modified Files: rtest_hypgeo.mac Log Message: Adjusting one example because of revision 1.68 of hypgeo.lisp. We have now an algorithm for bessel_k(0,a*t). Index: rtest_hypgeo.mac =================================================================== RCS file: /cvsroot/maxima/maxima/tests/rtest_hypgeo.mac,v retrieving revision 1.16 retrieving revision 1.17 diff -u -d -r1.16 -r1.17 --- rtest_hypgeo.mac 1 Oct 2009 18:04:14 -0000 1.16 +++ rtest_hypgeo.mac 12 Oct 2009 21:25:55 -0000 1.17 @@ -1899,14 +1899,12 @@ specint(%e^(-s*t)*bessel_y(0,a*t),t); -log((sqrt(s^2+a^2)+s)/(sqrt(s^2+a^2)-s))/(%pi*sqrt(s^2+a^2)); -/* Problem 247: No. 18 (OK noun form New) - We have no algorithm for bessel_k when the order is an integer or zero. - Tabulated result: - (log(s+sqrt(s^2-a^2))-log(a))/sqrt(s^2-a^2); -*/ +/* Problem 247: No. 18 OK New (10/2009) + * A special algorithm for bessel_k(0,a*t) has been added. + */ -radcan(specint(%e^(-s*t)*bessel_k(0,a*t),t)); specint(%e^(-s*t)*bessel_k(0,a*t),t); +acosh(s/a)/sqrt(s^2-a^2); /******************************************************************************* * |