Метод половинного деления, находим минимум унимодальной ф-ции Опции

[Sensitive]

Здравствуйте!
Мне необходимо найти минимум функции с помощью метода половинного деления на заданном отрезке (с точность Е)
Вот то, что получилось у меня, но работает не верно. Возможно Вы увидите ошибку в коде.

uses crt;

function f(x:real):real;
begin
f:=sqr(x-4)*sqr(x-4); {сама функция}
end;

var x,y:array [1..50] of real;
    eps,min,a,b:real;
    i,n:byte;

begin
clrscr;
write('vvedite a = '); readln(a); {начало отрезка}
write('vvedite b = '); readln(b); {конец отрезка}
write('vvedite eps = '); readln(eps); {точность}
min:=f(x[1]);
x[1]:=a;
x[5]:=b;
repeat

x[3]:=(x[5]-x[1])/2;
x[2]:=(x[3]-x[1])/2;
x[4]:=(x[5]-x[3])/2;

for i:=1 to 5 do
begin
y[i]:=f(x[i]);
if y[i]<min then
begin
min:=y[i]; n:=i;
end;
end;

x[1]:=x[n-1];
x[5]:=x[n+1]

until x[5]-x[1]<eps;

writeln('x = ',x[n]:0:5);
writeln('min = ',min:0:5);

readln;
end.

[Lapp]

Две поправки:

uses crt;

function f(x:real):real;
begin
f:=sqr(x-4)*sqr(x-4); {сама функция}
end;

var x,y:array [1..50] of real;
    eps,min,a,b:real;
    i,n:byte;

begin
clrscr;
write('vvedite a = '); readln(a); {начало отрезка}
write('vvedite b = '); readln(b); {конец отрезка}
write('vvedite eps = '); readln(eps); {точность}
min:=f(a);  { ---------------------------------------------------------------- here }
x[1]:=a;
x[5]:=b;
repeat

x[3]:=(x[5]-x[1])/2;
x[2]:=(x[3]-x[1])/2;
x[4]:=(x[5]-x[3])/2;

for i:=1 to 5 do begin
  y[i]:=f(x[i]);
  if y[i]<min then begin
    min:=y[i]; n:=i;
  end;
end;

if n=1 then n:=2 else if n=5 then n:=4;  { ------------------------------- here }

x[1]:=x[n-1];
x[5]:=x[n+1]

until x[5]-x[1]<eps;

writeln('x = ',x[n]:0:5);
writeln('min = ',min:0:5);

readln;
end.

Вроде так, но не проверял.