function M = odef1_mass(z)
  
  global C K R
  
  M = [1  0;  0  (C-K*z)*(1-R^2*(C-K*z)^2)];
  