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

public static void Main(string[] args){
double rmax=10,dr=0.2;
foreach(var arg in args){
	var words = arg.Split(':');
	if(words[0]=="-rmax")rmax=double.Parse(words[1]);
	if(words[0]=="-dr")dr=double.Parse(words[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];

jacobi.cyclic(H,e,V:V);
WriteLine($"# lowest energies");
for(int k=0;k<5;k++) WriteLine($"{e[k]} {-0.5/(k+1)/(k+1)}");
WriteLine("\n");

Func<double,double> f0 = z => 2*z*Exp(-z);
WriteLine("# ground state wave-function");
WriteLine($"{0} {0} {0}");
for(int i=0;i<n;i++){
	WriteLine($"{r[i]} {V[i,0]/Sqrt(dr)} {f0(r[i])}");
	}
Error.WriteLine($"{rmax} {dr} {e[0]} {e[1]} {e[2]}");
}//Main

}