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 );
>