Задача в паскале. Работа с циклом while do.

[Йохан Шмидт]

Здравствуйте! Я недавно начал изучать pascal. Взял для изучения книгу "Алгоритмизация и программирование на Турбо Паскале".
В 5 главе столкнулся с небольшой проблемой: Паскаль выдаёт ошибку. К тому же не совсем понимаю задание.




До цикла мы вычисляем погрешность
eps:=exp(ln(10)*(-3));

Потом сравниваем
abs(y1-y) и eps

В цикле мы вычисляем примерно следующее:
while abs(y1-y)<=eps do
begin
y:=y1;
y1:=y1+(1/n)*(((x/(exp(ln(y1)*(n-1)))))-y1);
end;

При попытке запуска, паскаль выдаёт ошибку:
error 207 invalid floating point operation
Если я правильно понимаю, то значения y,x равны нулю?
Как исправить проблему, подскажите пожалуйста.

[Stilet]

Если мне не изменяет память y1 должен быть больше нуля. Логарифм от отрицательного не бывает.
И кстати,

eps:=exp(ln(10)*(-3));

Дюж расточительно, я бы это убрал, а писал так:

while abs(y1-y)<=0.001 do

[Йохан Шмидт]

Спасибо, точно...

var
z,y,x,eps,y1:real;
i,n:integer;
Begin
clrscr;
eps:=0.001;
y1:=1;
y:=1;
x:=PI/3;
n:=1;
while abs(y1-y)<=eps do
begin
y:=y1;
y1:=y1+(1/n)*(x/(exp(ln(y1)*(n-1))-y1));
end;
readkey;
End.

Сейчас паскаль ругается на строчку
y1:=y1+(1/n)*(x/(exp(ln(y1)*(n-1))-y1));
Пишет что деление на 0. Как это возможно исправить?

[Sanprof]

попробуйте в блоке

try 

except

end

обрабатывать ошибку DivByZero

[Йохан Шмидт]

error 3 unknown identifier
Нужно что-то прописывать чтобы использовать try..except?

Я так понимаю, что деление на 0 происходит тут:
(x/(exp(ln(y1)*(n-1))))
Можно как-то без try..except исправить деление на 0. Вот вопрос как.
У меня такое чувство, что скобка где-то не так стоит. А где найти не могу.
А можно сделать так?

while abs(y1-y)<=eps do
begin
y:=y1;
z:=x/(exp(ln(y1)*(n-1)));
if z<>0 then
y1:=y1+(1/n)*(z-y1);
WriteLn(y1);
end;

[Ol'ga_new]

Ноль получается в скобке: n-1 на первом шаге. Наверное нужно просто начинать с 2!

[Йохан Шмидт]

Ну допустим... А как теперь определить корень 2 и 4 степени из числа PI/3?

[Stilet]

Exp(2*ln(pi/3)) и Exp(4*ln(pi/3)) если не ошибаюсь

[Йохан Шмидт]

Сообщение от Stilet

Exp(2*ln(pi/3)) и Exp(4*ln(pi/3)) если не ошибаюсь

Скорее Exp(ln(pi/3)/2) и Exp(ln(pi/3)/4)
Но это всё не то.
Наверняка это надо сделать как-то через цикл, тоесть вместо 2 и 4 подставить n. Я просто не совсем понимаю задание корректно. Может кто-то изложит более ясно?

помогите пожалуйста...

[Stilet]

Я просто не совсем понимаю задание корректно

А мы его тем паче не понимает не видя глазами