Radial schrodinger and shooting method (1 Viewer)

Users Who Are Viewing This Thread (Users: 0, Guests: 1)


I am to write a program using "shooting method"
to solve schrodinger radial equation with Coulomb potential.
I've already written the code but it doesn't work well.
I will be very thankful if someone would be so kind as to take
a look at my code(I can send it by e-mail) or give me some
advices what can be wrong.

To make everything clear :

I assume h=1, m=1

program finds function K(r) and then
R(r)(which is solution of schr. radial eq.) from the formula:
R(r) = K(r)/r

boundary conditions for K(r):
- close to zero K(r) behaves as r^(l+1)
- when values run to infinity the aproximation is : exp(-r*sqrt(-2*E))
and E is energy

I've implemented numerov algorithm and zero search method -
I've tried out this implementation for 1D schr. equation with some
potential (example taken from book "An introduction to
computational physics" Tao pang) and it seems to work well.

best regards

The Physics Forums Way

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving