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