int schint(double *Psi, double *Psip, double Psiout[], int k1, int k2, double Veff[], double E, double r[], double dr[], int N);Return value:
- number of nodes in solution between r[k1] and r[k2]
Output:
- Psi, Psip: values of
and
- Psiout[]: not used (yet!)
Input:
- Psi, Psip: values of
and
- k1, k2: indices of initial and final integration points (r[k1]
r[k2])
- Veff[]: effective potential on grid
- E: proposed electron energy
- r[], dr[], N: standard grid information
Modify your ``schint'' routine to conform to the above prototype. The modified routine (now a function subroutine if you are working in Fortran) should now return an integer equal to the number of ``nodes'' (zero crossings of the solution Psi) in the range from k1 to k2. Note that the routine also accepts an array Psiout, which we shall use in a later problem set.