Step 4: Solution to Initial Value Problem (optional)

This step should, by now, be fairly routine. It's a little more difficult to get Maple to help us because the derivative of the solution must be computed in a form in which the second initial condition can be substituted into the derivative of the solution.

> IC1 := y(0)= ? ;

> IC2 := D(y)(0)= ? ;

>

> EQ1 := value( subs( t=0, IC1, SOLN ) ); # apply first IC to the solution

>

> dSOLN := diff( SOLN, t );

> EQ2 := value( subs( t=0, IC2, convert(dSOLN,D) ) ); # apply second IC to the derivative

>

> PARAMic := solve( { EQ1, EQ2 }, { _C1, _C2 } );

> SOLNivp := subs( PARAMic, SOLN );

>