Orbital calculation FORTRAN program