function xdot=fcnspring(t,x,flag,m,c,k)
xdot=[x(2);-(c/m)*x(2)-(k/m)*x(1)];