function [M]=invPrandtlMeyerFunction(gamma,nu)

M_min=1;
M_max=10;

M_list=linspace(M_min,M_max,1000);
nu_list=PrandtlMeyerFunction(gamma,M_list);
[~,closestIndex] = min(abs(nu_list-nu));
M_min=M_list(closestIndex-1);
M_max=M_list(closestIndex+1);
M_list=linspace(M_min,M_max,1000);
nu_list=PrandtlMeyerFunction(gamma,M_list);
[~,closestIndex] = min(abs(nu_list-nu));
M=M_list(closestIndex);

end