%  Ideal gas constant 
function Rg = Rgas(M)        % val=J/kg/K
  R =8.3144621;              % universal gas constant
  Rg=1000*R/M;               % gas constant
end
