double rtbisp480(double (*func)(double,int,double [],double [],double[],int), double x1, double x2, double xacc, int kmatch, double Veff[], double r[], double dr[], int N)
Modify the Numerical Recipes bisection routine ``rtbis'' with the above prototype so that it (1) accepts and then (2) passes to all calls to (*func) the additional arguments kmatch, Veff[], r[], dr[] and N which ``func_Schrodinger'' requires. Also, make sure to remember to convert ``rtbis'' to double precision, for which you should also double the value of JMAX.
Finally, so that we can study the working of our code, modify the routine so that it prints out the final number of iterations and the final function value just before it returns rtb, with a print statement of the form
printf(" (%d iterations in rtbisp480, func=%e)\n",j,fmid);