using System;
using static System.Console;
using static System.Math;
class main{

public static void Main(string[] args){
double rmax=10,dr=0.2;
if(args.Length>0)rmax=double.Parse(args[0]);
if(args.Length>1)dr=double.Parse(args[1]);
int n=(int)(rmax/dr-1);
vector r = new vector(n);
vector e = new vector(n);
for(int i=0;i<n;i++)r[i]=dr*(i+1);
matrix H = new matrix(n,n);
matrix V = new matrix(n,n);
for(int i=0;i<n-1;i++){
  matrix.set(H,i,i,-2);
  matrix.set(H,i,i+1,1);
  matrix.set(H,i+1,i,1);
  }
matrix.set(H,n-1,n-1,-2);
H*=-0.5/dr/dr;
for(int i=0;i<n;i++)H[i,i]+=-1/r[i];
//var ev = new EVD(H);
jacobi.cyclic(H,e,V);
//for(int k=0;k<5;k++) WriteLine($"{ev.e[k]} {-0.5/(k+1)/(k+1)}");
for(int k=0;k<5;k++) WriteLine($"{e[k]} {-0.5/(k+1)/(k+1)}");
WriteLine("\n");
int ifac=5;
//double factor = ev.V[ifac,0]/(r[ifac]*Exp(-r[ifac]));
double factor = V[ifac,0]/(r[ifac]*Exp(-r[ifac]));
Func<double,double> f0 = z => z*Exp(-z)*factor;
WriteLine($"{0} {0} {0}");
for(int i=0;i<n;i++){
	//WriteLine($"{r[i]} {ev.V[i,0]} {f0(r[i])}");
	WriteLine($"{r[i]} {V[i,0]} {f0(r[i])}");
	}
Error.WriteLine($"{rmax} {dr} {e[0]} {e[1]} {e[2]}");
}//Main

}