Фрактал Ньютона(Pascal.)

[CopperNickel]

Задача.
Реализовать программу, строящую фрактал с заданными границами расчета x0<=x<=xn, y0<=y<=yn.
Где x0=-1.0, xn=0.7, y0=-0.7, yn=1.0.
Нашел на форуме тему 3 летней давности , те же данные все тоже, мне бы радоваться. Одно но, там ошибка в одном из корней. Сам решить не могу. Почитал про комплексные числа, про бассейн ньютона, комплексную функцию. Ошибка 207 что-то про плавающую запятую. Что я пробывал, я менял знаки, удалял часть переменных из строки, удалял саму строку, все чего добился половина фрактала и половина экрана желтая. Помогите плиз.

program N3;
uses Graph, Crt;
type Complex = Record
                x : Real;
                y : Real;
               end;
const iter = 50;
      max  = 1e+6;
      min  = 1e-6;
var  z, t, d : Complex;
     y, x, p : Real;
   xc, yc, n : Integer;
     Cancel  : Boolean;
     gd, gm  : Integer;
     mx, my  : real;
 x0,xn,y0,yn : real;
begin
Cancel := False;
Randomize;
gd := Detect;
InitGraph(gd,gm,'c:\bp\bgi');
{дипазон по х,у}
x0:=-1;xn:=0.7;
y0:=-0.7;yn:=1.0;
{начало экранных координат, смещено влево и вверх}
xc := round(GetMaxX/(xn-x0));
yc := round(GetMaxY*(-y0)/(yn-y0));
{масштаб для перевода координат в экранные}
mx:=xc;
my:=yc/-y0;
y:=y0;{цикл по у}
repeat
 x:=x0; {цикл по х}
  repeat
   n := 0;
   z.x := x;
   z.y := y;
   d := z;
   while (sqr(z.x)+sqr(z.y) < max) and (sqr(d.x)+sqr(d.y) > min)
   and (n < iter) do
    begin
     t := z;
     {z^3 - 1}
     p := sqr(sqr(t.x)+sqr(t.y));
     z.x := 2/3*t.x + (sqr(t.x)-sqr(t.y))/(3*p);{ в этой строке ошибка}
     z.y := 2/3*t.y*(1-t.x/p);{}
     d.x := abs(t.x - z.x);
     d.y := abs(t.y - z.y);
     Inc(n);
     if keypressed then Cancel := true;
    end;
   PutPixel(xc + round(x*mx),yc + round(y*my),16 - (n mod 16));
   if cancel then exit;
   x:=x+0.005;{если 0.001 то картинка лучше, но рисует долго}
  until x>xn;
 y:=y+0.005;
until y>yn;
Readkey;
end.

[Serge_Bliznykov]

так, бреда ради, попробуйте записать так:

    z.x := 2.0/3.0*t.x + (sqr(t.x)-sqr(t.y))/(3.0*p);{ в этой строке ошибка}
    z.y := 2.0/3.0*t.y*(1.0-t.x/p);{}

[CopperNickel]

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

так, бреда ради, попробуйте записать так:

    z.x := 2.0/3.0*t.x + (sqr(t.x)-sqr(t.y))/(3.0*p);{ в этой строке ошибка}
    z.y := 2.0/3.0*t.y*(1.0-t.x/p);{}

Попробовал, ничего не изменилось. Я удалял эту строку z.x. Дальше ошибок нет, удается запустить и получить полфрактала. Ошибка именно в этой строчке.

[ViktorR]

У меня Free Pascal.
Скопировал и запустил. Все сработало нормально, без ошибки.
Закоментировал строку, которая у Вас выдает ошибку.
Наблюдаю только левую половину фрактала.




Так что, ...

[CopperNickel]

У меня выдает только правую. Спасибо, попробую фри паскаль.

[Serge_Bliznykov]

странно. у меня ваш код из пост #1 сработал без ошибок.
Turbo Pascal 7.1.

[CopperNickel]

Спасибо всем я нашел ошибку.

[Serge_Bliznykov]

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

Спасибо всем я нашел ошибку.

интересно, и в чём же она заключалась?
Почему она проявлялась только у Вас?

[CopperNickel]

Там запятая вместо точки.

[Serge_Bliznykov]

простите... не понял...
где запятая вместо точки?! в вашем код из пост #1 ?