Abstract
A variation perturbation approach is presented for the Dirac equation of linear molecules. The method is tested on the lowest S1/2 state of H and σ1/2 state of H2+.