HYDROGEN-RADIAL
!--------------------------------------------------------------- !Program name: HYDROGEN-RADIAL !Radial equation of the hydrogen-like atom !Objective: plot the radial wave function for each n,l pair ! The equation is atomic units !--------------------------------------------------------------- !The Coulomb potential energy equation is: ! V(r) = -2Z/r ! The ceentrifugal potential equation is: ! VL=l(l+1)/r^2 ! The energy is: ! E=-2*Z/n^2 !--------------------------------------------------------------- ! Use of Numerov method to find the wave function ! Function fn from Numerov equation is ! Fn = 1 + K2*h^2/12 , where K2=-2(Veff - E) ! In the code, h is changed into dr !--------------------------------------------------------------- Program hydrogen Implicit none integer :: mesh, n, l, i, j, zeta, m, maxiter=100 Real*8 rmin, dr, zmesh, y_out_m double precision :: e, fh12, norm, eps=1.0E-6 double precision, allocatable :: r(:), r2