Посмотрите где может быть ошибка в коде

[student-99]

Реализация метода золотого сечения на Делфи.
a=-10; b=3; e=0,001
Вот код:

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, XPMan;

type
  TForm1 = class(TForm)
    Label1: TLabel;
    Label2: TLabel;
    Label3: TLabel;
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Edit4: TEdit;
    Button1: TButton;
    Label4: TLabel;
    XPManifest1: TXPManifest;
    Edit5: TEdit;
    Label5: TLabel;
    Button2: TButton;
    procedure Button1Click(Sender: TObject);
    procedure Button2Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var xm:real;
  Form1: TForm1;

implementation
function f(x:real):real;
begin
f:=3*sqr(x)*x-28*sqr(x)+2
end;

{$R *.dfm}

procedure TForm1.Button1Click(Sender: TObject);
var  a,b,e,opt:real;
     y,z,fy,fz:real;
     count:Integer;
begin
  e:=strtofloat(edit3.Text);
  a:=strtofloat(edit1.Text);
  b:=strtofloat(edit2.Text);
  count:=0;
  y:=0.618*a+0.382*b;
  z:=0.382*a+0.618*b;
  fy:=f(y);
  fz:=f(z);
    while  (b-a)>e do
    begin
    if f(y)<f(z) then
    begin
    b:=z;
    z:=y;
    fz:=fy;
    y:=0.618*a+0.382*b;
    fy:=f(y);
    end
    else
    a:=y;
    y:=z;
    fy:=fz;
    z:=0.382*a+0.618*b;
    fz:=f(z);
		Inc(count);
    end;
    xm:=(b+a)/2;
  opt:=Round(xm*10000)/10000;
  edit4.Text:=FloatToStr(opt);
  edit5.Text:=IntToStr(count)
end;

procedure TForm1.Button2Click(Sender: TObject);
begin
Close;
end;

end.

http://vk.com/doc50596190_163304908?...404363d7f323c3
Ссылка на архив.
Должно получится в конце около -0,2635
Спасибо.

[Sciv]

Что-то мне подсказывает, что в этом куске кода одного бегина-энда не хватает, не?

    else
    a:=y;
    y:=z;
    fy:=fz;
    z:=0.382*a+0.618*b;
    fz:=f(z);
		Inc(count);
    end;

[student-99]

Я поставил после else бегина-энд, но программа все равно подсчитывает оптимум не верно. Видимо ошибка где-то в алгоритме, но где не могу найти...((((

[student-99]

Вот блок схема по которой я создал программу:

[Sciv]

Ну во-первых, покажите, какой теперь код получился (после установки бегинэнда)
А во-вторых, где у Вас на блок-схеме вот это?

opt:=Round(xm*10000)/10000;

[student-99]

opt:=Round(xm*10000)/10000; - это округление к 4 знаков после комы...
А бегин я поставил после else
else
begin
a:=y;
y:=z;
fy:=fz;
z:=0.382*a+0.618*b;
fz:=f(z);
end;
Inc(count);
end;