I was interested in making what I thought would be a simple simulation of an electron encountering a positron by numerically solving the Schrodinger equation over several time steps, but I've run into 1 problem: the potential used in most textbooks (to solve the hydrogen atom for example) are still of a classical point particle, rather than a wavefunction. How do I generate the potential field from two real, moving and interacting wavefunctions?