Не работает программа в Pascal

[seriy_volk 228]

Возникла проблема при подготовке лабы в Pascal
Ошибку не выдаёт, лишь "NaN" в окне вывода

var s,x,a,b,c:real;y,z:integer;
  begin
    z:= 3;
    y:=1;
    b:=0;
    s:=0;
    x:=0.95;
    while z<=4 do
    begin
      while z<=5 do
        begin
        a:=(abs(sqr(sin(pi/y)))-exp(ln(cos((2*pi-x)/z))*3))/(1+(x*y)/(sqr(z-1)+sqr(z-1)));
        b:=b+a;
        z:=z+1
        end;
    c:=(ln(-(z+x)/8)+arctan((2*x)/z))/(1+exp(ln(cos((pi+x)/(2*z)))*3))-b;
    s:=s+c;
    y:=y+1
    end;
  writeln(s)
  end.

Прошу помочь, 2 дня до защиты

[p51x]

Зачем вам внешний цикл, если он выполняется один раз?
ln(-(z+x)/8) = ln(-(5+0)/8)=ln(-5/8). Чему равен логарифм отрицательного числа?

[digitalis]

Хорошо бы опубликовать исходное задание - понять, накосячено ли при переносе в код или само задание "с придурью" (бывает и такое)

[seriy_volk 228]

Большое спасибо, я совершенно забыл про ОДЗ
Если кому было интересно, прикладываю задания
Что касается внешнего цикла, то он выполняется 2 раза

[Аватар]

Сообщение от seriy_volk 228

Что касается внешнего цикла, то он выполняется 2 раза

Нет. 1 и тут z=6 и все, если раньше не сдохнет

c:=(ln(-(z+x)/8)+arctan((2*x)/z))/(1+exp(ln(cos((pi+x)/(2*z)))*3))-b;

[p51x]

А логарифм то зачем?

[digitalis]

А кубический корень чем извлекать будем ? Но там логарифм квадрата = положительного числа, так что с логарифмом затыка не будет.

[p51x]

Я про первый. Хотя и со вторым беда, в коде не квадрата, не корня.

[digitalis]

Я вот думаю: кому моча мысль в бошку стукнула - закручивать такие мастодонто-образные бессмысленные формулы. Попробую, думаю, хватит ли терпения. На Делфи.

procedure TForm1.Button1Click(Sender: TObject);
var z: integer ; x,s : Extended ;
//
function cmpt(z:integer;x:Extended) : Extended ;
var tmp1,tmp2,tmp3,tmp4: Extended ; y: integer ;
begin
  tmp1 := Exp((z+x)/-8)+ ArcTan (2*x/z) ;
  tmp3 := Cos((Pi+x)/(z+z)) ;
  tmp2 := 1.0 + Exp(Ln(tmp3*tmp3)/3) ;
//  Memo1.Lines.Add(Format('%9.6f %9.6f',[tmp3,tmp2])) ;
  tmp1 := tmp1/tmp2 ;
  tmp2 := 0 ;
  for y := 1 to 5 do
    begin
       tmp3 := Sin(Pi/y) ;
       tmp4 := Cos((2*Pi-x)/z) ;
       tmp3 := Abs(tmp3*tmp3-tmp4*tmp4*tmp4) ;
//       Memo1.Lines.Add(Format('%3d %9.6f %9.6f',[y,tmp3,tmp4])) ;
       tmp2 := tmp2 + tmp3/(1 + x*y/(z*z-z-z+y*y-y-y+2))
// в эту смешную скобку я раскрутил 2 квадратных трехчлена
    end ;
    cmpt := tmp1 - tmp2;
end ;
//
begin
  x := 0.95 ; s := cmpt(3,x) + cmpt(4,x) ;
  Memo1.Lines.Add(Format('= %9.6f',[s]))
end ;

Вроде правильно. И даже что-то выдает. Цены на кондомы в Антарктиде.

А он мне говорит: "Изобразите-ка квадратный трехчлен". А я не то что изобразить - вообразить такое не могу.