Метод итераций

[nik-kang]

Требуется написать программу для 2-х функций f(x) = x^2 и f(x) = sqrt(x)
Алгоритм для метода прилагается:


Пишу модуль:

unit Unit2;

interface
function F(x: real): real;
function IterF(a,b,e:real): real;
function G(x: real): real;
function IterG(a,b,e:real): real;

implementation

function f(x: real): real;
begin
result:=sqr(x);
end;

function IterF(a,b,e:real): real;
var h, x, y, m, c, x0, x1: real;
begin
h:=(b-a)/100;
x0:=(a+b)/2;
x:=a;
y:=a+h;
m:=Abs(F(y)-F(x))/h;
while y < b do
begin
x:=y;
y:=x+h;
c:=Abs(F(y)-F(x))/h;
if c > m then
m:=c;
end;
m:=m+1;
if f(a) < f(b) then
c := 1
else
c := -1;
x1:=x0-c*f(x0)/m;
while abs(x1-x0) > e do
begin
x0:=x1;
x1:=x0-c*F(x0)/m;
end;
result:=x1;
end;

function g(x: real): real;
begin
result:=sqrt(x);
end;

function IterG(a,b,e:real): real;
var h, x, y, m, c, x0, x1: real;
begin
h:=(b-a)/100;
x0:=(a+b)/2;
x:=a;
y:=a+h;
m:=Abs(G(y)-G(x))/h;
while y < b do
begin
x:=y;
y:=x+h;
c:=Abs(G(y)-G(x))/h;
if c > m then
m:=c;
end;
m:=m+1;
if G(a) < G(b) then
c := 1
else
c := -1;
x1:=x0-c*G(x0)/m;
while abs(x1-x0) > e do
begin
x0:=x1;
x1:=x0-c*G(x0)/m;
end;
result:=x1;
end;

end.

Но при вводе например данных в функцию f(x) = sqrt(x)
a: 0.1
b: 2
e: 0.001

выходит сообщение "Invalid floating point operation"
В чем ошибка?