Fortran 90 code for two-body 1D calculations.