Метод Хорд(Pascal)

[k4k7u3]

Подскажите как написать программу для нахождение точек пересечения графика функции, с осью X, с точностью до 0,0001 методом хорда?
x1(-10;0)
x2(0;10)
x3(10;20)
x4(20;30)
Функция задана уравнением:f(x)=a0*x^4+a1*x^3+a2*x^ 2+a3*x^+a4
x=-9,-9+2,5..29
a0=1
a1=-45
a2=593
a3=-1731
a4=-6210

[k4k7u3]

Подскажите, если найдете ошибку(


program lab_1;
var a,b,c,eps:real;i:integer;
function func(x:real):real;
begin
Func:=(x*x*x*x)-(45*x*x*x)+(593*x*x)-(1731*x)-6210;
end;
begin
read ();
i:=0;
eps:=0.0001;
for i:=0 to 3 do
begin
a:=-10+i*10;
b:=i*10;
writeln('a=',a,'b=',b);
if func(a)*func(b)<0 then
while abs(b-a)<=eps do
begin
c:=a-func(a)*(b-a)/(func(b)-func(a));
if func(a)*func(c)<0 then
b:=c else a:=c;

end;
writeln('x',i+1,'=',a:2:7);
end;
readln();
end.