function M = odef3_mass(eta)

  global C K R

if (eta == 1/R)
  M = eye(2);
else
  M = [1  0;  0  eta*(1-R^2*eta^2)];
end
