From: David B. <bil...@us...> - 2007-02-27 13:49:20
|
Update of /cvsroot/maxima/maxima/share/contrib/diffequations/tests In directory sc8-pr-cvs7.sourceforge.net:/tmp/cvs-serv2860 Modified Files: rtestode_kamke_1_2.mac Log Message: Remove special checks for Bessel function solutions. There is now code in ode_check() to do this. Index: rtestode_kamke_1_2.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/diffequations/tests/rtestode_kamke_1_2.mac,v retrieving revision 1.13 retrieving revision 1.14 diff -u -d -r1.13 -r1.14 --- rtestode_kamke_1_2.mac 26 Feb 2007 09:50:47 -0000 1.13 +++ rtestode_kamke_1_2.mac 27 Feb 2007 13:49:14 -0000 1.14 @@ -15,20 +15,6 @@ (pn_(n_):=print("Kamke ODE 1.",n_),true); true; -/* Another way to simplify Bessel function solutions */ -(besjysimp(e,n,x):=block([expr], - expr:subst(2*(n-1)*bessel_j(n-1,x)/x-bessel_j(n-2,x),bessel_j(n,x),e), - expr:subst(2*(n-1)*bessel_y(n-1,x)/x-bessel_y(n-2,x),bessel_y(n,x),expr), - ratsimp(expr)), 0); -0; -/* Apply besjysimp for order n and n-1 */ -(besjysimp2(e,n,x):=block([expr:besjysimp(e,n,x)],besjysimp(expr,n-1,x)), 0); -0; -/* Apply besjysimp for order n, n-1 and n-2 */ -(besjysimp3(e,n,x):=block([expr:besjysimp(e,n,x)], - expr:besjysimp(expr,n-1,x),besjysimp(expr,n-2,x)), 0); -0; - /* 101 */ (pn_(101),ans:contrib_ode(eqn:x*'diff(y,x)+x*y^2-y,y,x)); [(x^2*y-2*x)/(2*y)=%c]; @@ -475,7 +461,7 @@ +((2*bessel_y(-3/2,2*a/sqrt(-x))-2*bessel_y(1/2,2*a/sqrt(-x)))*%c -2*bessel_j(1/2,2*a/sqrt(-x))+2*bessel_j(-3/2,2*a/sqrt(-x)))*a*sqrt(-x)) /(4*bessel_y(-1/2,2*a/sqrt(-x))*%c+4*bessel_j(-1/2,2*a/sqrt(-x)))]; -[method,besjysimp3(ode_check(eqn,ans[1]),3/2,2*a/sqrt(-x))]; +[method,ode_check(eqn,ans[1])]; [riccati,0]; /* 165 */ |