43 lines
1.3 KiB
ObjectPascal
43 lines
1.3 KiB
ObjectPascal
|
// Демонстрация фрактальной графики. Папоротник
|
|||
|
|
// Для каждой точки комплексной плоскости z=(x,y) выполняем итерационный процесс z=z^2+c, c=(cx,cy)
|
|||
|
|
// Считаем количество итераций i до тех пор пока не выполнится условие |x|>max и |y|>max
|
|||
|
|
// После этого рисуем точку x,y с насыщенностью красного цвета, пропорциональной i
|
|||
|
|
uses GraphABCLinux;
|
|||
|
|
|
|||
|
|
const
|
|||
|
|
max = 10;
|
|||
|
|
cx = 0.251;
|
|||
|
|
cy = 0.95;
|
|||
|
|
coef1 = 0.5;
|
|||
|
|
coef2 = 0.88;
|
|||
|
|
|
|||
|
|
scalex = 0.001;
|
|||
|
|
scaley = 0.001;
|
|||
|
|
dx = 200;
|
|||
|
|
dy = 130;
|
|||
|
|
|
|||
|
|
begin
|
|||
|
|
Window.Title := 'Фракталы: папоротник';
|
|||
|
|
SetWindowSize(800,600);
|
|||
|
|
CenterWindow;
|
|||
|
|
for var ix:=0 to Window.Width-1 do
|
|||
|
|
for var iy:=0 to Window.Height-1 do
|
|||
|
|
begin
|
|||
|
|
var x := scalex*(ix-dx);
|
|||
|
|
var y := scaley*(iy-dy);
|
|||
|
|
var i := 1;
|
|||
|
|
while i<255 do
|
|||
|
|
begin
|
|||
|
|
var x1 := coef1*x*x-coef2*y*y+cx;
|
|||
|
|
var y1 := x*y+cy;
|
|||
|
|
x := x1;
|
|||
|
|
y := y1;
|
|||
|
|
if (abs(x)>max) and (abs(y)>max) then break;
|
|||
|
|
i += 1;
|
|||
|
|
end;
|
|||
|
|
if i>=255 then SetPixel(ix,iy,clRed)
|
|||
|
|
else SetPixel(ix,iy,RGB(255,255-i,255-i));
|
|||
|
|
end;
|
|||
|
|
writeln('Время расчета = ',Milliseconds/1000,' с');
|
|||
|
|
end.
|