function [Mv] = oun_Mv(gamma, Mm)
%oun_TvonTm
%
Mv=sqrt((2+(gamma-1)*Mm.^2)./(2*gamma*Mm.^2-(gamma-1)));
end