- #1

- 5

- 0

My question is that if I want to solve the Schrodinger equation for an infinite wall potential at r = R, I would want my wave function to die at r = R.

Here is the code that I know for a square well potential;

NDSolve[ f''[r] + 2*mu*Energy - 2*mu*VOO[r]f[r] == 0, f[0] == 1.,

f'[0] == 0.0]

However, I do not know how to define specially boundary conditions for the finite wall problem.

If you have any suggestion, I will appreciate it.

Thanks