Program numerov

  Program numerov
!-----------------------------------------------------------------
! April 2009
! 1D Schrodinger ODE y"(x)+2[E-V(x)]y(x)=0
! Numerical method: the Numerov's method for 2-nr oder ODE
! Root finding (e): the bisectional method
!-----------------------------------------------------------------
implicit none
integer, parameter :: n=101
double precision, dimension(-n:n) :: x, wl, wr
double precision e, a, b, emin, emax, energy, eps, fe, dx
double precision xmin, xmax
integer iflag, i
common/wf/ x, wl, wr
common/eqn/e, a, b, dx
external fe

open (unit=7,file="wfunc.dat")

 

https://ww2.odu.edu/~agodunov/computing/programs/f90/numerov.f90