Abstract
The determination procedure for shielded diatomic orbitals which are eigenfunctions of the two-center problem is generalized for any polynomial form of the separated parametric shielding potentials. The matrix technique using a Newton-Raphson scheme for simultaneous resolution of the separated equations, and the matrix elements of the general shielding potential are presented.