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