function zdot=fcnspring(t,z,flag,m,c,k)

%Forcing function
f=0;

%Compute the right hand side of the differential equations
zdot=[z(2);-(c/m)*z(2)-(k/m)*z(1)+f/m];