New routine: ``func_SchrodingerNodes''

double func_SchrodingerNodes(double E,
                              int kmatch, double Veff[],
                              double r[], double dr[], int N)

Return value:

• (number of nodes at energy E) - (kmatch)

Input:

• E: proposed electron energy
• kmatch: targeted number of nodes in solution
• Veff[]: effective potential on grid
• r[], dr[], N: standard grid information

Write a routine ``func_SchrodingerNodes'', which is very much like ``func_Schrodinger'', except that it returns the number of nodes reported by schint minus the value of the argument kmatch. In C, for instance, your return statement will appear as something like:

  return ((double) (schint(&Psi,&Psip,NULL,0,N,Veff,E,r,dr,N)-kmatch));

Note: Although the return value is always an integer, we return it as a double so that we can use ``rtbisp480'' routine without modification.