function dx=ode_model(t,x,FLAG,C,K)

%below is the system of first order ODEs for spring-mass-dashpot system
dx(1)=x(2);
dx(2)=-K*x(1)-C*x(2);
dx=dx';