From: David B. <bil...@us...> - 2007-02-27 13:47:21
|
Update of /cvsroot/maxima/maxima/share/contrib/diffequations/tests In directory sc8-pr-cvs7.sourceforge.net:/tmp/cvs-serv1116 Modified Files: rtestode_murphy_1_2.mac rtestode_murphy_1_4.mac Log Message: Remove special checks for Bessel function solutions. There is now code in ode_check() to do this. Index: rtestode_murphy_1_2.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/diffequations/tests/rtestode_murphy_1_2.mac,v retrieving revision 1.6 retrieving revision 1.7 diff -u -d -r1.6 -r1.7 --- rtestode_murphy_1_2.mac 24 Feb 2007 13:49:00 -0000 1.6 +++ rtestode_murphy_1_2.mac 27 Feb 2007 13:47:13 -0000 1.7 @@ -549,7 +549,7 @@ /* 173 Too general */ -/* 174 FIXME - 2007-02-04 - Solution correct. Add check. */ +/* 174 */ assume(a>0); [a>0]; (pn_(174),ans:contrib_ode(eqn:x*'diff(y,x)+a+x*y^2=0,y,x)); @@ -557,8 +557,8 @@ +bessel_j(2,2*sqrt(a)*sqrt(x))-bessel_j(0,2*sqrt(a)*sqrt(x)))*sqrt(a)*sqrt(x) -bessel_y(1,2*sqrt(a)*sqrt(x))*%c-bessel_j(1,2*sqrt(a)*sqrt(x))) /((2*bessel_y(1,2*sqrt(a)*sqrt(x))*%c+2*bessel_j(1,2*sqrt(a)*sqrt(x)))*x)]; -method; -riccati; +[method,ode_check(eqn,ans[1])]; +[riccati,0]; forget(a>0); [a>0]; @@ -626,8 +626,8 @@ (pn_(183),ans:contrib_ode(eqn:x*'diff(y,x)+a*x^2*y^2+2*y=b,y,x)); [y=(bessel_y(1,-sqrt(a)*sqrt(-b)*x)*%c+bessel_j(1,-sqrt(a)*sqrt(-b)*x))*sqrt(-b) /((bessel_y(0,-sqrt(a)*sqrt(-b)*x)*%c+bessel_j(0,-sqrt(a)*sqrt(-b)*x))*sqrt(a)*x)]; -method; -riccati; +[method,ode_check(eqn,ans[1])]; +[riccati,0]; /* 184 Riccati */ assume(not equal(m+n,0)); Index: rtestode_murphy_1_4.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/diffequations/tests/rtestode_murphy_1_4.mac,v retrieving revision 1.4 retrieving revision 1.5 diff -u -d -r1.4 -r1.5 --- rtestode_murphy_1_4.mac 24 Feb 2007 13:51:42 -0000 1.4 +++ rtestode_murphy_1_4.mac 27 Feb 2007 13:47:13 -0000 1.5 @@ -136,8 +136,8 @@ -bessel_j(n-1,2*sqrt(-x))) +%c*bessel_y(n-1,2*sqrt(-x))+bessel_j(n-1,2*sqrt(-x))) /(2*%c*bessel_y(n-1,2*sqrt(-x))+2*bessel_j(n-1,2*sqrt(-x)))]; -method; -riccati; +[method,ode_check(eqn,ans[1])]; +[riccati,0]; /* 387 Riccati */ (pn_(387),ans:contrib_ode(eqn:x^n*'diff(y,x)+x^(2*n-2)+y^2+(1-n)*x^(n-1)*y=0,y,x)); |