MapleDESeries5.mws

A Series Solution Comparison

for various degree approximations

Truncated series solutions are graphed below.

>    ode:=(v+5)*diff(y(v),v,v)-3*diff(y(v),v)-3*(y(v))=0;

ode := (v+5)*diff(y(v),`$`(v,2))-3*diff(y(v),v)-3*y(v) = 0

>    Order:=7;

Order := 7

>    dsolve({ode,y(0)=0,D(y)(0)=2},y(v),type=series);

y(v) = series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+O(v^7),v,7)

>    rhs(%);

series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+O(v^7),v,7)

>    poly:=convert(%,polynom);

poly := 2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6

>    with(plots):

Warning, the name changecoords has been redefined

>    Order7soln:=plot(poly,v=0..5,color=red):

>    display(Order7soln);

[Maple Plot]

>    Order:=12;

Order := 12

>    dsolve({ode,y(0)=0,D(y)(0)=2},y(v),type=series);

y(v) = series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11+O(v^12),v,12)

>    rhs(%);

series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11+O(v^12),v,12)

>    poly:=convert(%,polynom);

poly := 2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11

>    SeriesSoln:=plot(poly,v=0..5,color=blue):

>    display(SeriesSoln,Order7soln);

[Maple Plot]

>    Order:=18;

Order := 18

>    dsolve({ode,y(0)=0,D(y)(0)=2},y(v),type=series);

y(v) = series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/35035...
y(v) = series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/35035...

>    rhs(%);

series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/350350000000...
series(2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/350350000000...

>    poly:=convert(%,polynom);

poly := 2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/35035000000...
poly := 2*v+3/5*v^2+7/25*v^3+11/250*v^4+21/2500*v^5+3/5000*v^6+3/35000*v^7+1/1400000*v^9-1/14000000*v^10+9/770000000*v^11-13/7700000000*v^12+23/91000000000*v^13-381/10010000000000*v^14+2031/35035000000...

>    SeriesSoln:=plot(poly,v=0..5,color=blue):

>    display(SeriesSoln,Order7soln);

[Maple Plot]

>