double arshud(double g, double d, int *ierr);
Compute pressure factor based upon density factor .
ENTRY: Specific heat constant .
ENTRY: Density fraction .
EXIT: (optional) Return status code.
Specific heat ratio .
Given out of range.
Return Value
Pressure fraction or NAN on error.