using System;
public class hydrogen{

public static double Fe(double e, double r){
	double rmin = 1e-3;
	if(r<rmin) return r-r*r;

	Func<double,vector,vector> shrodinger = (x,y)=>
	{ /* -(1/2)f'' - (1/r)f = e f */
		return new vector(new double[] { y[1], 2*(-1/x-e) *y[0] });
	};

	vector ymin = new vector(new double[]{rmin-rmin*rmin, 1-2*rmin});
	vector ymax = ode.rk23(shrodinger,rmin,ymin,r);
	return ymax[0];
}

}//class