function dy = ode_model(t,y,C,K)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = -K*y(1) - C*y(2);