Что не так в программе? программа для вычисления вектора X по заданной формуле. Паскаль.

[motivator]

Вот задача.
Вот решение.

program zad_3;
function tg(x:real):real;
begin
tg:=0;
if cos(x)=0 then writeLn(' ошибка при расчете tg x=',x)
else
tg:=sin(x)/cos(x)
end;
var
j, m: integer;
x: array[1..10] of real;
begin
readln(m);
readln(x[1], x[2], x[3]);
for j:=4 to m do
begin
x[j]:=power(tg(x[j-1]),1/3)-ln(abs(x[j-1]));
writeln(x[j]:3:3);
end;
end.

Что не так, выводит наны.

[Вадим Мошев]

Рекомендую посмотреть на значение каждой переменной, которая участвует в расчётах.
NaN Возникает в том случае, когда:
1. Делится нуль на нуль
2. Из отрицательного числа извлекается корень чётной степени
3. Аргумент логарифма - нуль или отрицателен.

И вообще, в тех случаях, когда математическое выражение нельзя вычислить в силу недопустимости аргумента (за исключением деления на нуль НЕ нуля)

[motivator]

Даже если не везде нан выводит
6
8
2
4
-0.336
NaN
NaN
????

[ura_111]

А ну, проанализируй фотку:

12.jpg



Что ты думаешь по этому поводу?

[Serge_Bliznykov]

попробуйте такой код:

program zad_3;

function tg(x: real): real;
begin
  tg := 0;
  if cos(x) = 0 then writeLn('#error# tg(',x,') undefined.')
  else
    tg := sin(x) / cos(x)
end;

function sqrt3(x:real):real;
var absx : real;
begin
   if x<0 then sqrt3 := -1*Power(abs(x),1/3)
   else sqrt3 := Power(x,1/3)
end;

var
  k, m: integer;
  x: array[1..10] of real;

begin
  readln(m);
  readln(x[1]);
  for k := 2 to m do
  begin
    x[k] := sqrt3(tg(x[k - 1])) - ln(abs(x[k - 1]));
    writeln(x[k]);
  end;
end.

[motivator]

Большое спасибо, Сергей.

[Serge_Bliznykov]

Не за что! Успехов!