Numerical recipes in FORTRAN 90 See original record