FORTRAN subroutines for minimization by quasi-Newton methods