function rhs = odef1_rhs(z,y)

  global C K R

  rhs = [y(2);   K^2*(C-K*z)*y(1) - 2*K*y(2)];
