Numerical methods for Fortran programmers