1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:

(******************************************************************************)
(******************************************************************************)

type

TComplex = record
 x,
 y : double;
end;

TQuat = record
 x,
 y,
 z,
 w : double;
end;

(******************************************************************************)
procedure Quat2Comp(const Q : TQuat; var C1, C2 : TComplex);
asm
  fld  TQuat(Q).X;// Delphi fld Q.X
  fsub TQuat(Q).w;
  fstp TComplex(C1).x;

  fld  TQuat(Q).y;
  fadd TQuat(Q).z;
  fstp TComplex(C1).y;

  fld  TQuat(Q).x;
  fadd TQuat(Q).w;
  fstp TComplex(C2).x;

  fld  TQuat(Q).y;
  fsub TQuat(Q).z;
  fstp TComplex(C2).y;
end;

Begin
end.					

1:

result.x := ((a.x*b.x)-(a.w*b.w)) / (sqr(b.x)-sqr(b.w)) + ((a.y*b.y)-(a.z*b.z)) / (sqr(b.y)-sqr(b.z));					

1:
2:
3:
4:
5:
6:
7:

const
  b = 640;
  h = 640;
  _x1 = -1.02;// Bei -1.0 kommt DivisionbyZeroError
  _x2 =  1.0;
  _y1 = -7.0;
  _y2 =  7.0;					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:

function rec2polC(const C : TComplex) : TComplex;
begin
 result := C;

 result.x := sqrt(sqr(result.x) + sqr(result.y));

 if result.x = 0 then
  begin
//Was denn nun?
    result.y :=  PIhalbe;
    result.y := -PIhalbe;
  end
 else
  begin
    result.y := arctan(result.y / result.x);
    if (result.x < 0)then result.y := result.y + pi;
  end;
end;					

1:
2:

result.x := hypot(c.x, c.y);
result.y := arctan2(c.y, c.x);					

1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:

//Zeichenprozedur
procedure jetztzeichnen;
var
    T1,T0: TDateTime;
    maxwx,diffx,diffy:real;

    i,j,b,h,plotfarbe:integer;
    red,green,blue,maxf,maxw:integer;
    k1:string;

begin
    //Funktionsterm
    k1:=edit0.text;
    termkorrektur(k1);

    b:=bitmap.width;
    h:=bitmap.height;
    diffx:=(_x2-_x1)/b;
    diffy:=(_y2-_y1)/h;
    //Zeichnen
    T0 := Time;
    for i:=0 to b do
    begin
      kx.re:=_x1+i*diffx;
      for j:=0 to h do
      begin
        kx.im:=_y1+j*diffy;
        //komplexer Funktionswert
        fehlerkomplex:=0;
        ky:=funktionswertkomplex(k1,kx);

        //Phase ermitteln
        maxwx:=arctan2(ky.im,ky.re);
        if maxwx<0 then maxwx:=maxwx+2*pi;

        //Zeichenfarbe auswählen
        if fehlerkomplex<>0 then plotfarbe:=clblack
        else
        begin
              maxw := trunc(3/pi*maxwx);
              maxf := round(3*255/pi*(maxwx-maxw*pi/3));
              case maxw of
                0 : begin Red := 255; Green := 0; Blue := maxf end;
                1 : begin Red := 255-maxf; Green := 0; Blue := 255 end;
                2 : begin Red := 0; Green := maxf; Blue := 255 end;
                3 : begin Red := 0; Green := 255; Blue := 255-maxf end;
                4 : begin Red := maxf; Green := 255; Blue := 0 end;
               else begin Red := 255; Green := 255-maxf; Blue := 0 end;
              end;
            plotfarbe:=rgb(red,green,blue);
        end;
        //bitmap.Canvas.pixels[i,h-j]:=plotfarbe;
      end;
//      paintbox1.canvas.draw(0,0,bitmap);
//      application.processmessages;
    end;
    T1 := Time;
    Bitmap.Canvas.TextOut(10,10,FormatDateTime('HH:NN:SS.ZZZ ',T1-T0));
end;