function [rhovonrhom] = oun_rhovonrhom(gamma,Mm )
%oun_rhovonrhom
%
rhovonrhom=(gamma+1)*Mm.^2./(2+(gamma-1)*Mm.^2);
end