formula for the f2c fortran-to-C tool from Netlib