%Diese Funktion gibt ein Bild der Konvergenzgeschwindigkeit aus
function draw_nits(xrange,yrange,tol)

ImageMatrix=zeros(length(xrange),length(yrange));
ns=zeros(length(xrange)*length(yrange),2);

for i=1:length(xrange)
    for j=1:length(yrange)
        %[x,y] soll nun immer ein Punkt im von xrange und yrange
        %definierten Gitter sein
        
        x=xrange(i);
        y=yrange(j);
        
        %%%%%%%%%%%%%%%%%%%%%%%%%%%
        %AUFRUF IHRER NEWTON FKT  % 
        %kann z.B wie folgt heissen%
        %%%%%%%%%%%%%%%%%%%%%%%%%%%
        [x_null,n_it]=newton2d([x;y],tol);
        
        %merke Nullstelle
        ns((i-1)*length(yrange)+j,:)=x_null';
        
        %Schreibe bis zur Konvergenz benoetigte Anzahl in Matrix
        ImageMatrix(length(yrange)-j+1,i)=n_it;
    end
end

clf;

%Setze Farbtabelle und plotte
colormap(hot);
imagesc(ImageMatrix);

axis off;

%plotte noch die Nullstellen
figure(2);
hold on;
for i=1:length(xrange)*length(yrange)
   plot(ns(i,1),ns(i,2),'X','MarkerSize', 24);
end

end