Здравствуйте!
Мне необходимо найти минимум функции с помощью метода половинного деления на заданном отрезке (с точность Е)
Вот то, что получилось у меня, но работает не верно. Возможно Вы увидите ошибку в коде.
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.
Две поправки:
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.