Доработка программы. Вывести количество итераций

[vitaly38]

Добрый день, уважаемые пользователи.
Суть программы
решение нелинейного уравнения методом деления отрезка пополам.
Сама программа готова, но не могу сделать вывод количества итераций (для каждого корня)
вот сам код

unit Unit1;

interface

uses
  Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
  Dialogs, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, Buttons, Series,math;

type
  TForm1 = class(TForm)
    Edit1: TEdit;
    Edit2: TEdit;
    Edit3: TEdit;
    Edit4: TEdit;
    Edit5: TEdit;
    Chart1: TChart;
    Label2: TLabel;
    Label3: TLabel;
    Label4: TLabel;
    Label5: TLabel;
    Memo1: TMemo;
    Button1: TButton;
    BitBtn1: TBitBtn;
    Series1: TLineSeries;
    procedure Button1Click(Sender: TObject);
  private
    { Private declarations }
  public
    { Public declarations }
  end;

var
  Form1: TForm1;

implementation

{$R *.dfm}

type fun=function(x:extended):extended;

function f(x:extended):extended;
begin
f:=x-7*power(sin(x),2)
end;

function delpop(a,b,eps:extended; f:fun; var i:integer):extended;
var x1,x2,y1,y2,xm:extended;

begin
 i:=0;
repeat

x1:=(a+b-eps)/2; y1:=f(x1);
x2:=(a+b-eps)/2; y2:=f(x2);
if y1>y2
  then a:=x1
  else b:=x2;

   until abs(b-a)<=eps;
       i:=i+1;
 xm:=(a+b)/2;

result:=xm;

end;

procedure TForm1.Button1Click(Sender: TObject);
var a,b,h,eps,x:extended;
i,n:integer;
begin

series1.Clear;
a:=strtofloat(edit1.Text);
b:=strtofloat(edit2.Text);
h:=strtofloat(edit3.Text);
eps:=strtofloat(edit4.Text);
x:=a;
repeat
series1.AddXY(x,f(x),'',clmaroon);
x:=x+strtofloat(edit5.Text);
until x>b+0.00000001;

x:=a; n:=0; memo1.Clear;
repeat
if (f(x)<f(x-h))and(f(x)<f(x+h))then
begin inc(n);
memo1.Lines.Add(inttostr(n)+'-й минимум= '+
floattostrf(delpop(x-h,x+h,eps,f,i),fffixed,9,4)+'    '+inttostr(i));
end;
x:=x+h;
until x>b+0.0000000001;

end;

end.

количество итераций- параметр i, но почему то программа выдает значение 1 всегда для каждого корня. В чем проблема? почему счетчик не работает? заранее благодарю.

[Grag]

Так ты параметр i увеличиваешь на 1 за пределами цикла !!! Вот он и равен всегда 1 (i=0+1)...

[vovk]

 i:=0;
repeat
x1:=(a+b-eps)/2; y1:=f(x1);
x2:=(a+b-eps)/2; y2:=f(x2);
if y1>y2
  then a:=x1
  else b:=x2;
   until abs(b-a)<=eps;// тут цикл кончается
      i:=i+1; // это уже вне цикла
 xm:=(a+b)/2;

ну и сами понимаете
для i выполняется при каждом вызове delpop

i:=0;
i:=i+1;

и сколько получится?

//Опоздал

[vitaly38]

я пробовал и перед антилом, и везде вставлял, все равно 1 равно.

[Grag]

В тело цикла вставлять надо! В тело, а не за его пределами!

[vitaly38]

Grag
пробовал!!!

[vitaly38]

не понимаю .. почему счетчик не работает.. уже везде его пробовал ставить..