Dr. Kevin G. TeBeest
Assoc. Prof. of Applied Mathematics
Kettering University

Example 7:    Pure Resonant Motion

> de := diff( x(t),t$2 ) + 4 * x(t) = 24 * sin(2*t) ;

Solve the above ODE, subject to zero initial position and zero initial velocity. Store the result in X:

> X := rhs( dsolve( { de, x(0) = 0, D(x)(0) = 0 }, x(t) ) );

Simplify the solution, and store the result in X:

> X := combine( X, trig ) ;

Plot the solution:

> plot( X, t = 0 .. 15, x = -100 .. 100, thickness = 6, color = blue, title = `Pure Resonance`, ytickmarks = 3, size = [560,400] ) ;