Изменить программу

@Стас45rus · Регистрация: 12.04.2010

Помогите изменить программу, чтобы получить фрактал на прикреплённом рисунке. Формула рисунка - ch(z) -гиперболический косинус. Формула гиперболического косинуса ch(z)=(e^z+e^-z)/2. Я выделил действительную и мнимую часть, но получается совсем не то. Что я не правильно делаю?

Pascal

PROGRAM Zadanie_14a;
Uses Graph,Crt;
CONST
  iter=50;
  max=16;
TYPE
  TComplex=Record
            x:Real;
            y:Real;
           End;
VAR
  z,t,c:TComplex;
  x,y,n:Integer;
  Cancel:Boolean;
  gd,gm:Integer;
  mx,my:Integer;
BEGIN
 ClrScr;
 Cancel:=False;
 Randomize;
 gd:=Detect;
 InitGraph(gd,gm,'');
 mx:=GetMaxX div 2;
 my:=GetMaxY div 2;
 For y:=-my To my Do
  For x:=-mx To mx Do
    begin
      n:=0;
      z.x:=x*0.005;
      z.y:=y*0.005;
      c.x:=0.11;
      c.y:=-0.66;
      While (sqr(z.x)+sqr(z.y)<max) and (n<iter) Do
       begin
         t:=z;
         z.x:=(exp(t.x+t.y)+exp(-t.x-t.y))/2;
         z.y:=(exp(t.x+t.y)+exp(-(t.x+t.y)))/2;
         Inc(n);
         If KeyPressed Then Cancel:=True;
       end;
      If n<iter Then PutPixel(mx+x,my+y,8+(n mod 16));
      If Cancel Then Exit;
    end;
 ReadKey;
 CloseGraph;
END.

@Стас45rus · 11.06.2012, 14:49 **[ТС]**

Переделал программу. Выделил из ch z действительную и мнимую часть. Формула: ch z=ch x*cos y+i*sh x*sin y. Re: ch x*cos y; Im:sh x*sin y. ch x=(e^x+e^-x)/2; sh x=(e^x-e^-x)/2. Итоговая формула: Re(e^x+e^-x)/2)*cos y; Im(e^x-e^-x)/2)*sin y. Подставил в код но рисунок совсем другой. Уважаемые программисты помогите.

Pascal

PROGRAM Zadanie_14a;
Uses Graph,Crt;
CONST
  iter=50;
  max=16;
TYPE
  TComplex=Record
            x:Real;
            y:Real;
           End;
VAR
  z,t,c:TComplex;
  x,y,n:Integer;
  Cancel:Boolean;
  gd,gm:Integer;
  mx,my:Integer;
BEGIN
 ClrScr;
 Cancel:=False;
 Randomize;
 gd:=Detect;
 InitGraph(gd,gm,'');
 mx:=GetMaxX div 2;
 my:=GetMaxY div 2;
 For y:=-my To my Do
  For x:=-mx To mx Do
    begin
      n:=0;
      z.x:=x*0.05;
      z.y:=y*0.05;
      c.x:=0.11;
      c.y:=-0.66;
      While (sqr(z.x)+sqr(z.y)<max) and (n<iter) Do
       begin
         t:=z;
         z.x:=((exp(t.x)+exp(t.x))/2)*cos(t.y);
         z.y:=((exp(t.x)-exp(t.x))/2)*sin(t.y);
         Inc(n);
         If KeyPressed Then Cancel:=True;
       end;
      If n<iter Then PutPixel(mx+x,my+y,8+(n mod 16));
      If Cancel Then Exit;
    end;
 ReadKey;
 CloseGraph;
END.

Опции темы