A is  an 2N by 2N Matrix
•For the purposes of this presentation a short FORTRN 90 program 
•to fill up the matrix elements of the A matrix was written.
•
•
   Once one has the eigenvalue matrix equation AY=λY. An IMSL, NAG, etc.
    subroutine can be called to solve for the eigenvector Y that
   contains the wave functions  u(x) and w(x) and the egenvalue λ form
   which the Binding Energy of the deuteron can be obtained.
    The wave functions are normalized  by requiring
A typical plot of the u(x) and w(x) wave functions is shown
 the next slide.