  We present a multigrid scheme for the solution of finite-element Hartree-Fock
equations for diatomic molecules. It is shown to be fast and accurate, the time
effort depending linearly on the number of variables. Results are given for the
molecules LiH, BH, N_2 and for the Be atom in our molecular grid which agrees
very well with accurate values from an atomic code. Highest accuracies were
obtained by applying an extrapolation scheme; we compare with other numerical
methods. For N_2 we get an accuracy below 1 nHartree.
