%
%  spring_rhs
%

  function dzc = spring_rhs(t,z,params)

  C = params(1);
  K = params(2);
  A = [0 1;-K -C];

  dzc = A*z;