MapleDESystemIC.mws

Entering Initial Conditions, Solving The System

For The PowerPoint Presentation

A = C[1], B = C[2], C = C[3], D = C[4]

>    f:=A*exp(3*x)+B*exp(-2*x)+C*sin(2*x)+D*cos(2*x)+(1/65)*x*exp(3*x)-(15/104)*x*sin(2*x)+(3/104)*x*cos(2*x);

f := A*exp(3*x)+B*exp(-2*x)+C*sin(2*x)+D*cos(2*x)+1/65*x*exp(3*x)-15/104*x*sin(2*x)+3/104*x*cos(2*x)

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

f0 := A+B+D

>    fx:=diff(f,x);

fx := 3*A*exp(3*x)-2*B*exp(-2*x)+2*C*cos(2*x)-2*D*sin(2*x)+1/65*exp(3*x)+3/65*x*exp(3*x)-15/104*sin(2*x)-15/52*x*cos(2*x)+3/104*cos(2*x)-3/52*x*sin(2*x)

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

fx0 := 3*A-2*B+2*C+23/520

>    fxx:=diff(fx,x);

fxx := 9*A*exp(3*x)+4*B*exp(-2*x)-4*C*sin(2*x)-4*D*cos(2*x)+6/65*exp(3*x)+9/65*x*exp(3*x)-15/26*cos(2*x)+15/26*x*sin(2*x)-3/26*sin(2*x)-3/26*x*cos(2*x)

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

fxx0 := 9*A+4*B-63/130-4*D

>    fxxx:=diff(fxx,x);

fxxx := 27*A*exp(3*x)-8*B*exp(-2*x)-8*C*cos(2*x)+8*D*sin(2*x)+27/65*exp(3*x)+27/65*x*exp(3*x)+45/26*sin(2*x)+15/13*x*cos(2*x)-9/26*cos(2*x)+3/13*x*sin(2*x)

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

fxxx0 := 27*A-8*B-8*C+9/130

>    solve({f0=0,fx0=1,fxx0=1,fxxx0=1},{A,B,C,D});

{A = 502/4225, C = 395/1352, B = -3/400, D = -301/2704}

>