MapleDEinitSolve.mws

DE Example 2

Putting in initial conditions and solving the resulting system

>    f:=A+B*x+C*exp(-2*x)+D*exp(-x);

f := A+B*x+C*exp(-2*x)+D*exp(-x)

>    f0:=eval(f,x=0);

f0 := A+C+D

>    fx:=diff(f,x);

fx := B-2*C*exp(-2*x)-D*exp(-x)

>    fx0:=eval(fx,x=0);

fx0 := B-2*C-D

>    fxx:=diff(fx,x);

fxx := 4*C*exp(-2*x)+D*exp(-x)

>    fxx0:=eval(fxx,x=0);

fxx0 := 4*C+D

>    fxxx:=diff(fxx,x);

fxxx := -8*C*exp(-2*x)-D*exp(-x)

>    fxxx0:=eval(fxxx,x=0);

fxxx0 := -8*C-D

>    solve({f0=0,fx0=4,fxx0=-6,fxxx0=14},{A,B,C,D});

{A = 0, B = 2, C = -2, D = 2}

>