using Laitone pressure approximationdouble armcrla(double g, double cp0, int *ierr);
Compute critical mach number given
Laitone pressure approximation.
ENTRY: Specific heat constant
.
ENTRY: Critical pressure coefficent
.
EXIT: (optional) Return status code.
Specific heat ratio
.
Iterative solver failed.
Return Value
Critical mach number
or NAN on error.