метод итераций по образцу с форума, выдает ошибку Опции

[Ven]

нужно решить уравнение методом простых итераций. Задачу нашла на форуме в FAQ


сделала по образцу, в итоге ошибка 205: Floating Point Overlow

program ITERAT;
uses crt;

const max_iter=100; {maksimal'noe chislo iteraciy}

var

i: integer;
x,x0,eps,M: real;

function F(x:real):real; {funkciya}

begin

F:= x*x + x + 0.09;
end;

begin {osnovnaya programma}
clrscr;
write('Vvedite priblijennoe znachenie x='); readln(x);
write('Vvedite tochnost vichisleniy eps='); readln(eps);

i:=0;
repeat
x0:=x;
x:=x0+F(x0);
inc(i);
writeln('---Iteraciya  ' , i:3,'  x=',x);
writeln('F(x)=' , F(x),' tochnost=' , abs(x-x0));
until (abs(x-x0)<=eps) or (i>max_iter);

if (abs(x-x0)<=eps) then writeln ('Otvet: X=' , x)
else writeln('OTVET NE NAYDEN!!! za ',max_iter:0,' shagov iteraciya ne soshlas');

end.

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

Добавлено через 15 мин.
кстати,когда выдает ошибку - возращает меня на строку

F:= x*x + x + 0.09;

Сообщение отредактировано: Ven - 9.10.2011 22:21

[IUnknown]

Ven, я не понимаю, в чем проблема? Ты задаешь точность в одну сотую, то есть хочешь, чтобы в ответе были как минимум 2 цифры после запятой правильны. Так? Получаешь
Otvet : X=-9.01...E-0001
, то есть, если перейти от экспоненциального представления к более привычному: - 0.901... Все, как заказывала, 2 цифры после запятой правильные...

Когда точность задана в одну тысячную - значит, надо чтоб как минимум три цифры после запятой были точными. Что имеем?
Otvet : X=-9.00075...E-0001
, или -0.900075... Опять все правильно. Три цифры точны, дальше - как повезёт. Тебе повезло, точны даже 4. Где подвох? Что не так?

Печатай результат не в экспоненциальной форме, а в обычной:

if (abs(x-x0)<=eps) then writeln ('Otvet: X=', x:15 :10)

, не надо будет переводить в привычный вид, сразу будет печататься в нём.

Сообщение отредактировано: IUnknown - 16.01.2012 21:16