Symbolic polynomial approximation of an LDE solution

Copyright © 2004-2012  DigiArea, Inc. . All rights reserved.

Description:

This worksheet illustrates ©LdeApprox ©package capability of doing symbolic polynomial approximation of an LDE solution. First of all we load LdeApprox package and define an IVP. Then ©we use ApproxSol procedure to find 3-rd degree polynomial approximation for the IVP solution on interval x = [-1,1]. After that we find exact solution by Maple™ procedure dsolve . Finally we compare exact and approximate results using Maple™ procedure plot3d .

This loads the package.

restart:
with(LdeApprox):

Example

Initial value problem

ivp:={diff(y(x),x,x,x)+lambda*x*diff(y(x),x)+lambda^2*y(x) = 0,y(0)=0,D(y)(0)=0,`@@`(D,2)(y)(0)=lambda};

ivp := {diff(y(x),`$`(x,3))+lambda*x*diff(y(x),x)+lambda^2*y(x) = 0, D(y)(0) = 0, `@@`(D,2)(y)(0) = lambda, y(0) = 0}

Finding polynomial approximation for solution of the IVP.

apr:=ApproxSol(ivp,y(x),x=-1..1, 3);

apr := y(x) = -2304*lambda^3*(2+lambda)/(41472*lambda^3+23040*lambda^4+11*lambda^6+6*lambda^7+lambda^8-46080*lambda^2+6*lambda^5+84934656)+48*lambda^2*(9216+4608*lambda+6*lambda^3+5*lambda^4+lambda^5)/...
apr := y(x) = -2304*lambda^3*(2+lambda)/(41472*lambda^3+23040*lambda^4+11*lambda^6+6*lambda^7+lambda^8-46080*lambda^2+6*lambda^5+84934656)+48*lambda^2*(9216+4608*lambda+6*lambda^3+5*lambda^4+lambda^5)/...
apr := y(x) = -2304*lambda^3*(2+lambda)/(41472*lambda^3+23040*lambda^4+11*lambda^6+6*lambda^7+lambda^8-46080*lambda^2+6*lambda^5+84934656)+48*lambda^2*(9216+4608*lambda+6*lambda^3+5*lambda^4+lambda^5)/...
apr := y(x) = -2304*lambda^3*(2+lambda)/(41472*lambda^3+23040*lambda^4+11*lambda^6+6*lambda^7+lambda^8-46080*lambda^2+6*lambda^5+84934656)+48*lambda^2*(9216+4608*lambda+6*lambda^3+5*lambda^4+lambda^5)/...

Finding exact solution of the IVP using Maple™ ©function dsolve .

sol:=dsolve(ivp,y(x));

sol := y(x) = 1/2*lambda*x^2*hypergeom([2/3+1/3*lambda],[4/3, 5/3],-1/9*lambda*x^3)

Comparing exact and approximate results.

plot3d(subs(sol,y(x))-subs(apr,y(x)),lambda=-1/8..1/8,x=-1..1,axes=boxed,orientation=[45,45]);

[Maple Plot]

Note

The method applied in the package is numerically - analytical one. It means that you can use symbolic expressions ©as boundary conditions, interval of approximation etc. However these kind of examples leads to huge output so its not for Web. This reason force us to introduce simple example with one parameter only. You can try more complex examples in your computer.