In the present work, a mass equation obtained as the solution of an inhomogeneous partial difference equation is used to predict masses of unknown neutron-rich and proton-rich nuclei. The inhomogeneous source terms contain shell-dependent symmetry energy expressions (quadratic in isospin), and include, as well, an independently derived shell-model Coulomb energy equation which describes all known Coulomb displacement energies with a standard deviation of [sigma]c = 41 keV. Perturbations of higher order in isospin, previously recognized as a cause of systematic effects in long-range mass extrapolations, are also incorporated. The most general solutions of the inhomogeneous difference equation have been deduced from a [chi]2-minimization proc...