double 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.