double armcrc(double g, double cp, int *ierr);
Compute given the critical mach number given .
ENTRY: Specific heat constant .
ENTRY: Critical pressure coefficent .
EXIT: (optional) Return status code.
Specific heat ratio .
Given out of range.
Iterative solver failed.
Return Value
Critical mach number or NAN on error.