% Testfunktion für numerov()
% Grundzustandsenergie des Valenzelektrons eines Natrium Atoms
Z = 11; %Natrium Ladungszahl
y(1) = 0;
%y(2) = h;
k1y1 = 2* Z;
anf = 0;
aend = 3.9405;
n = 6400;
h = (aend - anf)/n;
y(2) = h;
%eigw suche
eigw = -0.75;
d_eigw = 0.1;
f_alt = -1;
i = 1;
while (abs(f_alt) > 0.0000000001 & i < 100)
    i = i+1;
    y = numerov(anf,aend,y(1),y(2),k1y1,eigw,n,'natrium_elektron_kfunc');
    f = y(n+1) - aend * (y(n+2)-y(n))/(2*h);
    if (f*f_alt) < 0   % bei f-vorzeichenwechsel
        d_eigw = - d_eigw / 10;
    end
    eigw = eigw + d_eigw;
    f_alt = f;
    disp(num2str(eigw,10))
end
%plot(i,y,'r.-');
