Recent content by serelha

You are right, but L=0 here. So the boundary conditions are correct. If L>0, then the w.f vanishes at the origin.

Here is the Schrödinger equation I mentioned Manipulate[ m = 6.; R = 5.; Vs2 = 4./(2*m*R^2); Vs = -10./(2*m*R^2) + Vs2; VCC[r_] = Vs*UnitStep[R - r] + Vs2*UnitStep[r - R]; L = 0; system = {RC''[r] + 2/r*RC'[r] + (-L*(L + 1)/r^2 - 2*mu*(VCC[r] - Energy))*RC[r] == 0...

Thanks for your respond.. Regarding the boundary conditions; my original differential equation is Bessel differential equations, so as you guessed I did substitution R[r]=RC[r]/r, however, the boundary conditions are correct for L = 0 (s-wave). The conditions that you suggested are correct...

Thanks for your help. I want to solve the whole system numerically, and I have already tried the method that you suggested. In this problem since I want to calculate the phase shift, I can vary the energy as much as I want. Here is a sample code mu = 6.; width = 5.; R = width; Vb =...

I am trying to use NDSolve to solve a Schrödinger equation. My question is that if I want to solve the Schrödinger 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[...