[Maxima-commits] CVS: maxima/tests rtestode.mac,NONE,1.1 Makefile.am,1.9,1.10 tests.lisp,1.12,1.13

[David B. <bil...@us...>]

Update of /cvsroot/maxima/maxima/tests
In directory usw-pr-cvs1:/tmp/cvs-serv10162

Modified Files:
	Makefile.am tests.lisp 
Added Files:
	rtestode.mac 
Log Message:
Add ODE tests in rtestode.mac

--- NEW FILE: rtestode.mac ---
/* ODE tests */

kill(all);
done;

/* Examples from "The Maxima Book" */

ode2(x^2*'diff(y,x)+3*x*y=sin(x)/x, y, x);
y = (%C-cos(x))/x^3;
ic1(%, x=1, y=1);
y = -(cos(x)-cos(1)-1)/x^3;
method;
linear;

ode2('diff(y,x,2) + y = 4*x, y, x);
y = %K1*sin(x) + %K2*cos(x) + 4*x;
method;
variationofparameters;

ode2((3*x^2+4*x+2)=(2*y-1)*'diff(y,x), y, x);
y^2-y = x^3+2*x^2+2*x+%C;
method;
separable;

/* ode2(x^2*cos(x*y)*'diff(y,x) + (sin(x*y)+x*y*(cos(x*y)))=0, y, x);
x*sin(x*y)=%C;
method;
exact; */

/* ode2( (2*x*y-exp(-2*y))*'diff(y,x)=0, y, x);
x*exp(2*y) - log(y) = %C;
method;
exact; 
intfactor;
exp(2*y)/y; */

/* ode2( 'diff(y,x)=(y/x)^2+2*(y/x), y, x);
-(x*y+x^2)/y = %C;
method;
exact;*/ /*homogeneous;*/

ode2( 'diff(y,x)+(2/x)*y=(1/x^2)*y^3, y, x);
y = 1/(sqrt( 2/(5*x^5) + %C)*x^2);
method;
bernoulli;
odeindex;
3;

ode2( 'diff(y,x,2)-3*'diff(y,x)+2*y=0, y, x);
y = %K1*exp(2*x) + %K2*exp(x);
method;
constcoeff;

soln9: ode2( 'diff(y,x,2)-4*'diff(y,x)+4*y=0, y, x);
y = (%K2*x + %K1)*exp(2*x);
method;
constcoeff;

ode2(x^2*'diff(y,x,2)+x*'diff(y,x)-y=0, y, x);
y=%K2*x-%K1/(2*x);
method;
exact;

ode2( x^2*'diff(y,x,2)+4*x*'diff(y,x)+2*y=0, y, x);
y=%K1/x+%K2/x^2;
method;
exact; /*euler*/

ode2( x^2*'diff(y,x,2)+5*x*'diff(y,x)+4*y=0, y, x);
y=(%K2*log(x)+%K1)/x^2;
method;
euler;

ode2( x^2*'diff(y,x,2)+x*'diff(y,x)+(x^2-1/4)*y=0, y, x);
y=(%K1*sin(x)+%K2*cos(x))/sqrt(x);
method;
bessel;

ode2( x^2*'diff(y,x,2)+x*'diff(y,x)+(x^2-4)*y=0, y, x);
y=%K1*%J[2](x)+%K2*%Y[2](x);
method;
bessel;

ode2( (x-1)^2*'diff(y,x,2)+(x-1)*'diff(y,x)+((x-1)^2-4)*y=0, y, x);
y=%K1*%J[2](x-1)+%K2*%Y[2](x-1);
method;
bessel;

ode2( 'diff(y,x,2)+2*'diff(y,x)+y=exp(x), y, x);
y=exp(x)/4+(%K2*x+%K1)*exp(-x);
method;
variationofparameters;
/* Particular solution */
yp;
exp(x)/4;

/*ode2( x*'diff(y,x,2)+('diff(y,x))^2=0, y, x);
y='integrate(1/(log(x)+%K1),x)+%K2;
method;
freeofy; */

/* ode2( y*'diff(y,x,2)+('diff(y,x))^2=0, y, x);
y^2/(2*%K1)=x+%K2;
method;
freeofx; */


Index: Makefile.am
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/Makefile.am,v
retrieving revision 1.9
retrieving revision 1.10
diff -u -d -r1.9 -r1.10
--- Makefile.am	30 Aug 2002 13:36:58 -0000	1.9
+++ Makefile.am	20 Sep 2002 04:40:34 -0000	1.10
@@ -53,6 +53,7 @@
 rtest8.mac \
 rtest9a.mac \
 rtest9.mac \
+rtestode.mac \
 tests.lisp
 
 EXTRA_DIST = $(genericdirDATA) summarize-log known-bugs

Index: tests.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/tests/tests.lisp,v
retrieving revision 1.12
retrieving revision 1.13
diff -u -d -r1.12 -r1.13
--- tests.lisp	3 Jul 2002 02:13:30 -0000	1.12
+++ tests.lisp	20 Sep 2002 04:40:34 -0000	1.13
@@ -29,7 +29,7 @@
 	  "rtest8"
 	  "rtest9"
 	  "rtest9a" "rtest10" "rtest11" "rtest12" "rtest13" "rtest13s"
-	  "rtest14" "rtest15" "rexamples"
+	  "rtest14" "rtest15" "rexamples" "rtestode"
 	  )
 	do
 	(format t "~%Testing ~a.mac" testv)