//consider a NACA-0012 airfoil
Cpomin=-0.43;//minimum pressure coefficient on the surface of airfoil at low speed from figure of Cp vs x/c given in question.
M=linspace(0.4,0.9,6);//Mach number over which we have to calculate Cp critical.
y=1.4;//specific heat ratio for air.