function ex408(f,m,c,k)
    dydt=@(t,y) [y(2);(sin(2*pi*f*t)-c*y(2)-k*y(1))/m];
    [t,y]=ode45(dydt,[0,20],[0;0]);
    plot(t,y(:,1));grid('on');
    title(['f=',num2str(f),' m=',num2str(m),' c=',num2str(c),' k=',num2str(k)]);
    ylabel('Displacement[m]');