Numerical methods in FORTRAN