function [M2] = machshock(M1,gamma);

M2 = sqrt(((gamma-1)*M1^2 + 2)/(2*gamma*M1^2 - (gamma - 1)));

end