and 
SUBROUTINE AR21AB(G,M,CP,RESULT,IERR) DOUBLE PRECISION G, M, CP, RESULT(3) INTEGER IERR
Compute the Laitone pressure correction coefficient relations given the mach number
and the the Laitone pressure
correction coefficient
.
ENTRY: Specific heat constant
.
ENTRY: Mach number
.
ENTRY: Pressure correction coefficient
.
EXIT: Array with result properties as described in 19.
EXIT: (optional) Return status code.
Specific heat ratio
.
Given
out of range.