Skriv in diffekv
> 
> de:=diff(y(t),t,t)+16*y(t)=6*sin(w*t);


                     / 2      \
                     |d       |
               de := |--- y(t)| + 16 y(t) = 6 sin(w t)
                     |  2     |
                     \dt      /

Lös ekv med dsolve (se?dsolve)
> svar:=dsolve({de,y(0)=0,D(y)(0)=0},y(t));

                                 w sin(4 t)   6 sin(w t)
              svar := y(t) = 3/2 ---------- - ----------
                                         2            2
                                  -16 + w      -16 + w

Sen med w=4
> w:=4:
> svar4:=dsolve({de,y(0)=0,D(y)(0)=0},y(t));

           svar4 := y(t) = - 3/4 cos(4 t) t + 3/16 sin(4 t)

Ta bort värdet 4
> w:='w':
plotta (höger led (rhs) av) svaret med olika värden på w.
> plot(subs(w=3.5,rhs(svar)),t=0..40,y=-10..10);
> plot(rhs(svar4),t=0..40);                              
> 
Vill man se hur losningen andras med w kan man anvanda kommandot animate (se ?animate). Då måsta man först ladda in fler grafikrutiner
> with(plots):
> animate(rhs(svar),t=0..40,w=0..6,frames=30,numpoints=300);
>