Ошибка в методе простых итераций

program chm;
var
me,kol,i,j:integer;
k,q,n,b,c,y0,x0,x1,y1,r1,r2,r3,r4,nach,kon,shag,toch,x:real;
a:array[1..100,1..100] of real;
it:array[1..3,1..2] of real;
function fun(x:real):real;
begin
fun:=x*x-2*x;
end;
function fun1(x:real):real;
begin
fun1:=2*x-2;
end;
function iter(var k,x:real):real;
begin
iter:=x-((x*x*x-5*x*x+3*x+7)/k);
 end;
function funrk(var x0,y0:real):real;
begin
funrk:=1-sin(1.25*x0+y0)+0.5*y0/(x0+2);
end;
procedure dihotomia(var n,k:real; var t:real);
var ser,r:real;
  begin
  while k-n>=2*t do
         begin
         ser:=(k+n)/2;
         if fun(n)*fun(ser)<=0 then
         k:=ser
         else n:=ser;
         end;
  ser:=(k+n)/2;
  writeln('koren=', ser:3:3,' s tochnostyu ',t:2:4);
  end;
    procedure hord(var n,k:real; var t:real);
        begin
        if n>=0 then
            begin
            while k-n>=t do begin n:=n-(fun(n)*(k-n)/(fun(k)-fun(n))); writeln(n,' ',k,' '); end;
            writeln('koren =',n:3:3, ' s tochnostyu ',t:3:3);
            end
        else
        begin
        while k-n>=t do begin   k:=k-(fun(k)*(k-n)/(fun(k)-fun(n))); end;
        writeln('koren =',k:3:3, ' s tochnostyu ',t:3:3);
        end;
    end;
    procedure kasat(var n,k:real; var t:real);
        begin
        if n>=0 then
            begin
            while k-n>=t do begin k:=k-(fun(k)/fun1(k)); end;
            writeln('koren =',n:3:3, ' s tochnostyu ',t:3:3);
            end
        else
        begin
        while k-n>=t do begin   n:=n-(fun(n)/fun1(n)); end;
        writeln('koren =',k:3:3, ' s tochnostyu ',t:3:3);
        end;
    end;
 
 
 
begin
writeln('1-3:  x*x-2*x');
writeln('4:  0.5*(x+(1/x))');
writeln('6:  S((lg(4+2*x*x)/(x-1))');
writeln('1- Metod dihotomii');
writeln('2- Metod hord');
writeln('3- Metod kasatelnih');
writeln('4- Metod prostih iteraciy');
writeln('5- Formula Lagranga');
writeln('6- Formula Simpsona');
writeln('7- Metod Runge-Kutta');
writeln('viberete metod');
readln(me);
if (me>=1) and (me<=3) then
    begin
    writeln('Vvedite: nachalo otrezka, konec otrezka, shag, tochnost');
    readln(nach,kon,shag,toch);
    n:=nach;
    k:=n+shag;
    kol:=0;
    i:=1;
    while k<=kon do
        begin
        if fun(n)*fun(k)<=0 then
              begin
              kol:=kol+1;
              a[i,1]:=n;
              a[i,2]:=k;
              i:=i+1;
              end;
        n:=n+shag;
        k:=k+shag;
        end;
    for i:=1 to kol do
    for j:=1 to 2 do write(a[i,j]:3:4,' ');
    end
else if (me=4) then
    begin
    writeln('vvedite tochnost');
    readln(toch);
    end
     else if(me=6) then
        begin
        writeln('vvedite: nachalo, konets, kollichestvo shagov');
        readln(nach,kon,kol);
        end
          else if (me=7) then
    begin
    writeln('vvedite: nachalo otrezka, konets otrezka, Y0, kollichestvo shagov');
    readln(nach,kon,y0,kol);
    end;
case me of
    1:begin
       for i:=1 to kol do
         begin
         n:=a[i,1];
         k:=a[i,2];
         dihotomia(n,k,toch);
         end;
      end;
    2:begin
      for i:=1 to kol do
            begin
            if ((a[i,1]<0) and (a[i,2]>0)) then
                    begin
                    kol:=kol+1;
                    a[kol+1,1]:=0;
                    a[kol+1,2]:=a[i,2];
                    a[i,2]:=0;
                    end;
            end;
      for i:=1 to kol do hord(a[i,1],a[i,2],toch);
      end;
3:begin
      for i:=1 to kol do
            begin
            if ((a[i,1]<0) and (a[i,2]>0)) then
                    begin
                    kol:=kol+1;
                    a[kol+1,1]:=0;
                    a[kol+1,2]:=a[i,2];
                    a[i,2]:=0;
                    end;
            end;
      for i:=1 to kol do kasat(a[i,1],a[i,2],toch);
      end;
4:      begin
    it[1,1]:=10;
        it[2,1]:=-3;
        it[3,1]:=10;
        it[1,2]:=1/3;
        it[2,2]:=3;
        it[3,2]:=4;
 
        for i:=1 to 3 do
            begin
            n:=iter(it[i,1],it[i,2]); writeln(n:4:4,'     ');
            k:=iter(it[i,1],n); write(k:4:4,'    ');
   repeat
            n:=k; writeln(n:4:4,'     ');
            k:=iter(it[i,1],n); write(k:4:4,'    ');
            until (abs(k-n)<toch);
            writeln('x= ',k:4:4);
            end;
        end;
5:
        begin
        writeln('vvedite kollichestvo uzlov');
        readln(kol);
        for i:=1 to kol do
            begin
            writeln('vvedite ',i,' uzel i znachenie funkcii v nem');
            readln(a[i,1],a[i,2]);
            end;
        writeln('vvedite X ');
        readln(x);
        nach:=0;
        toch:=1;
        for i:=1 to kol do
            begin
            for j:=1 to kol do  if i<>j then
            toch:=toch*(x-a[j,1])/(a[i,1]-a[j,1]);
            nach:=nach+toch*a[i,2];
            toch:=1;
            end;
        writeln('Ln(x)=',nach:4:4);
        end;
6:
        begin
        x:=0;
        q:=0;
        shag:=(kon-nach)/kol;
        for i:=1 to (kol+1) do
            begin
            a[i,1]:=(((ln(1+nach*nach))/(ln(10)))/(2*nach-1));
            writeln(a[i,1]:4:4);
            nach:=nach+shag;
            end;
        k:=a[1,1]+a[kol+1,1];
        for i:=2 to kol do if (i mod 2)=0 then x:=x+a[i,1]
        else q:=q+a[i,1];
        writeln(x:4:4,q:4:4,k:4:4);
        k:=(k+4*x+2*q)*shag/3;
        writeln('znachenie integrala = ',k:4:4);
        end;
7:
        begin
        shag:=(kon-nach)/kol;
        x0:=nach;
 
        writeln('x0= ',x1:4:4,'     y0= ',y1:4:4);
        for i:=1 to kol do
            begin
            q:=x0+shag/2;
            x:=y0+r1/2;
            toch:=y0+r2/2;
            b:=x0+shag;
            c:=y0+r3;
            r1:=shag*funrk(x0,y0);
            r2:=shag*funrk(q,x);
            r3:=shag*funrk(q,toch);
            r4:=shag*funrk(b,c);
            writeln(r1:3:3,r2:3:3,r3:3:3,r4:3:3);
            y1:=y0+((r1+2*r2+2*r3+r4)/6);
            x1:=x0+shag;
            writeln('x',i,'= ',x1:4:4,'     y',i,'= ',y1:4:4);
            y0:=y0+((r1+2*r2+2*r3+r4)/6);
            x0:=x0+shag;
            end;
        end;
end;
end.