готов лаби 4 курс 1 сим / цеом / lilu2

program Fractal;
uses windows;
type
BMP_HEADER=record
bfType:WORD; {BM}
FileSize:DWORD;
NoneUse:DWORD;
BytesToData:DWORD;
end;
BitMapInfo=record
StructSize:DWORD;
bmpWidth:DWORD;
bmpHeight:DWORD;
Planes:WORD;
BitPerPixel:WORD;
TypeCompression:DWORD;
DataSize:DWORD;
XPelsPerMeter:DWORD;
YPelsPerMeter:DWORD;
ColorUsed:DWORD;
ColorImportant:DWORD;
end;
PixelType=record
B,G,R:BYTE;
end;
Ekran =array [0..512*512-1] of PixelType;
PixEkr=array [1..512,1..512] of PixelType;
Complex=record
Re,Rez1,Rez2,Rez4,Rez3:single;
Im:single;
end;
const
Shapka:array [1..54] of byte=($42,$4D,$36,0,$C,0,0,0,0
,0,$36,0,0,0,$28,0,0,0,0,2,0,0,0,2,0,0,1,0,$18
,0,0,0,0,0,0,0,$C,0,$12,$B,0,0,$12,$B,0,0,0,0,
0,0,0,0,0,0);
var
Zagolovok:BMP_HEADER;
Opis :BitMapInfo;
Pixels :^PixEkr;
Palette :array [1..256] of PixelType;
f1:file;
T,x,y:DWORD;
Odin,C,Z,A,D:Complex;
minx,miny,maxx,maxy,tx:single;
n:DWORD;
Begin
GetMem(Pixels,3*512*512);
Odin.Re:= 1.0;
Odin.Im:=-1.0;
minx:=-1.5;
maxx:= 0.5;
miny:=-1.0;
maxy:= 1.0;
for n:=1 to 256 do
begin
Palette[n].R:=255-abs(127-n);
Palette[n].B:=256-n;
Palette[n].G:=n-1;
end;
T:=GetTickCount;
for x:=1 to 512 do {}
begin
for y:=1 to 512 do
begin
C.Re:=minx+(maxx-minx)*x/512.0;
C.Im:=miny+(maxy-miny)*y/512.0;
Z:=C;
n:=0;
repeat
begin
{A.Re:=C.Re-((sqr(Z.Re)-sqr(Z.Im))*Z.Re-2*Z.Re*sqr(Z.Im)-1);
A.Im:=C.Im-(Z.Im*(sqr(Z.Re)-sqr(Z.Im))+2*sqr(Z.Re)*Z.Im);
Z.Re:=A.Re;
Z.Im:=A.Im;
tx:=sqr(Z.Re)+sqr(Z.Im);}
asm

fld Z.Re {Z.Re;-;-;-;-;-;-;-}
fmul Z.Re {Z.Re*Z.Re;-;-;-;-;-;-;-}
fld Z.Im {Z.Im;Z.Re*Z.Re;-;-;-;-;-;-}
fmul Z.Im {Z.Im*Z.Im;Z.Re*Z.Re;-;-;-;-;-;-}
fsub {Z.Re*Z.Re-Z.Im*Z.Im;-;-;-;-;-;-;-}
fld Z.Im {Z.Im;Z.Re2-Z.Im2+C.Re;-;-;-;-;-;-}
fmul Z.Re {Z.Im*Z.Re;Z.Re2-Z.Im2+C.Re;-;-;-;-;-;-}
fadd st,st {2*Z.Im*Z.Re;Z.Re2-Z.Im2+C.Re;-;-;-;-;-}
fstp A.Im {Z.Re2-Z.Im2+C.Re;-;-;-;-;-;-;-}
fstp A.Re {Z.Re;-;-;-;-;-;-;-}


fld A.Re {A.Re;-;-;-;-;-;-;-}
fmul Z.Re {A.Re*B.Re;-;-;-;-;-;-;-}
fld A.Im {A.Im;A.Re*B.Re;-;-;-;-;-;-}
fmul Z.Im {A.Im*B.Im;A.Re*B.Re;-;-;-;-;-;-}
fsub {A.Re*B.Re-A.Im*B.Im;-;-;-;-;-;-;-}
fstp D.Re {-;-;-;-;-;-;-;-}
fld A.Im {A.Im;-;-;-;-;-;-;-}
fmul Z.Re {A.Im*B.Re;-;-;-;-;-;-;-}
fld A.Re {A.Re;A.Im*B.Re;-;-;-;-;-;-}
fmul Z.Im {A.Re*B.Im;A.Im*B.Re;-;-;-;-;-;-}
fadd {A.Im*B.Re+A.Re*B.Im;-;-;-;-;-;-;-}
fstp D.Im {-;-;-;-;-;-;-;-}


fld C.Re
fld D.Re
fsub
fld1
fadd

fld C.Im
fld D.Im
fsub

fstp Z.Im
fst Z.Re
fmul st,st {Z.Re2;-;-;-;-;-;-;-}
fld Z.Im {Z.Im;Z.Re2;-;-;-;-;-;-}
fmul st,st {Z.Im2;Z.Re2;-;-;-;-;-;-}
fadd {Z.Im2+Z.Re2;-;-;-;-;-;-;-}
fstp tx {-;-;-;-;-;-;-;-}

end;
end;

n:=n+1;
until( (tx>4.0) or (n=255) );{}
Pixels^[y,x]:=Palette[n];
end;
end;
T:=GetTickCount-T;
Writeln('Time calckulate is T=',T,' ms.');
assign(f1,'lilu2.bmp');
rewrite(f1,1);
blockwrite(f1,Shapka,54);
blockwrite(f1,Pixels^,3*512*512);
FreeMem(Pixels,3*512*512);
close(f1);
readln;
End.