From: David B. <bil...@us...> - 2004-02-14 11:17:33
|
Update of /cvsroot/maxima/maxima/share/contrib/rand In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv21236 Modified Files: mathieu0.mac Log Message: Removing control characters Index: mathieu0.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/contrib/rand/mathieu0.mac,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- mathieu0.mac 14 Feb 2004 11:10:48 -0000 1.1 +++ mathieu0.mac 14 Feb 2004 11:10:59 -0000 1.2 @@ -1,4 +1,4 @@ -/* Filename <name>.mac +/* Filename mathieu0.mac *************************************************************** * * @@ -12,7 +12,7 @@ * * *************************************************************** */ /* - (d4) This program computes the transition curve through +(d4) This program computes the transition curve through the origin (N = 0) in Mathieu's equation using a perturbation @@ -25,7 +25,7 @@ */ - mathieu0():=(input(),setup1(),setup2(), +mathieu0():=(input(),setup1(),setup2(), for i thru m do (step1(),step2a(),step3a()),output())$ input():=m:read("ENTER DEGREE OF TRUNCATION")$ setup1():=(w:0,for i thru m do w:w+k[i]*e^i,x:1, @@ -38,4 +38,3 @@ step2a():=(temp1:ode2(temp1,y[i](t),t),temp1:ev(temp1,%k1:0,%k2:0))$ step3a():=(f[i]:solve(coeff(expand(rhs(temp1)),t,2),k[i]),e[i]:ev(temp1,f[i]))$ output():=(print("delta=",ev(w,makelist([f[j]],j,1,m))),print(" "))$ - |