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