Численные методы. Решение нелинейного уравнения методом половинного деления.

[gree]

Привет помогите плиз доразобраться, Короче говоря задание нужно написать программу вычисляющую корень уравнения на заданном промежутке, с заданой точностью.. по данной блок схеме..

уравнение и его корни можно увидеть из графика.

Среда программирования Delphi.
вот что у меня получилось:

procedure TForm1.Button1Click(Sender: TObject);
var a,b,c,eps,x:real;
begin
a:=strtofloat(edit1.Text);
b:=strtofloat(edit2.Text);
eps:=strtofloat(edit3.Text);
 if (sin(0.5*a+0.2)-sqr(a)*cos(0.5*a+0.2))*(sin(0.5*b+0.2)-sqr(b)*cos(0.5*b+0.2))<0 then
begin
edit4.Text:='числа заданы не верно';
exit;
end;
repeat
c:=(a+b)/2;
if (sin(0.5*a+0.2)-sqr(a)*cos(0.5*a+0.2))*(sin(0.5*c+0.2)-sqr(c)*cos(0.5*c+0.2))<0 then
b:=c
else a:=c
until abs(b-a)<0;
x:=(a+b)/2;
edit4.Text:=floattostr(x);
end;
end.

[Ti_pain)]

until abs(b-a)<0;
помойму ето не правельно модуль любого чилса не меньше 0))