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

%system of first order ODE of spring-mass-dashpot system

dx=zeros(2,1);
dx(1)=x(2);
dx(2)=-K*x(1)-C*x(2);