func_SchrodingerNodes()

double func_SchrodingerNodes(double E,
                             int match, double V[],
                             double r[], double dr[], int N)

Return value:

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

Input:

• E: proposed electron energy
• match: desired number of nodes in solution
• V[]: potential on grid
• r[], dr[], N: standard grid information

Write a routine func_SchrodingerNodes() in ``physics.c'' (don't forget to put the prototype in ``physics.h'' - last reminder!), which is very much like func_Schrodinger(), except that it returns the number of nodes reported by schint minus the value of the argument ``match''. Your return statement, for instance, will appear something like:

  return ((double) (schint(&Psi,&Psip,Psiout,k1=0,k2=N,V,E,r,dr,N)-match));

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