A dipole fit to electromagnetic form factors is used to determine a qu
ark density rho in the nucleon. A radial tenser potential is used to b
ind the quarks into states of good J, J(Z), and parity. The tenser pot
ential radial component is taken to satisfy the equation T' = T-0 rho,
where T-0 is a parameter of the model. This linear divergence equatio
n can be simultaneously solved with the Dirac equation for the bound q
uark wave functions. A self-consistent solution is possible where the
mass density used as the source for the binding potential is the same
as that determined from the solution for the quark wave functions.