A two-dimensional, fully numerical approach to the four-component first-ord
er Dirac-differential equation utilizing the Finite-Element-Method (FEM) is
employed for H-2(+) and Th-2(179+) Using elliptic-hyperbolic coordinates a
nd further one-dimensional singular transformations, scaling transformation
s and extrapolation techniques (geometrical over iteration steps and logari
thmic over grid points) ve achieve for the molecules H-2(+) and Th-2(179+)
relative accuracies better than 10(-12) for 1(112)g energies.