``zriddrp480''

double zriddrp480(double
	      (*func)(double,int,double[],double[],double[],int),
	      double x1, double x2, double xacc,
	      int kmatch, double Veff[], double r[], double dr[], int N)

Once you have ``rtbisp480'' working, make the same modifications to ``zriddr'', making certain to printout the final number of iterations needed and final function value.

Also, change your main program to print out the results for E both from ``rtbisp480'' and ``zriddrp480'' as well...

            .
            .
            .
E=rtbisp480(func_Schrodinger,E1,E2,1e-15,
    kmatch,Veff,r,dr,N);
printf("N=%6d E=%20.15f\n",N,E);

E=zriddrp480(func_Schrodinger,E1,E2,1e-15,
    kmatch,Veff,r,dr,N);
printf("N=%6d E=%20.15f\n",N,E);
            .
            .
            .