Numerical recipes in Fortran 77 See original record