|
ru.algorithms- RU.ALGORITHMS ---------------------------------------------------------------- From : Michael Ryazanov 2:5030/1072.2 21 May 2001 21:40:00 To : Slava Arharov Subject : Re: Фрактал -------------------------------------------------------------------------------- 20.05 22:29 Slava Arharov, 2:5090/96.15 -> All SA> Может уже избито, но не подскажет ли мне олл _HОРМАЛЪHЫЙ_ и _РАБОТАЮЩИЙ_ SA> алгоритм сабжа Манделъброта ? Какие-то сабжи я уже строитъ научился, но SA> этот... Hикак не получится, хотъ ты тресни...Можно исходник, толъко SA> работающий, но так, чтоб алгоритм явно просматривался. А то у меня естъ SA> куча сурсов, где 97% занимает процесс записи в банки VESA и прочая чушъ. SA> Можно просто записатъ в массив цвета точек, а на экран сам выведу. Заранее SA> благодарствую. ___/~~~mandelbr.ps~~~\___ __ _ . /iter { % на входе -- координаты точки c /y exch def /x exch def /x0 x def /y0 y def /n 0 def 20 -1 0 { % 20 -- количество итераций /re x dup mul def /im y dup mul def re im add 4 ge { % |z| >= 2 -- выход за пределы круга => /n exch def % => последовательность расходится exit % n -- номер итерации } if pop /y 2 x y mul mul y0 add def /x re im sub x0 add def } for n 20 div % нормализация } def 220 170 translate 100 100 scale /dx 0.02 def % шаг по x /dy 0.02 def % y -1.2 dx 0 { -2 dy 0.6 { dup 2 index iter setgray 2 copy 1 index dx dy rectfill neg dx dy rectfill % симметрично } for pop } for showpage ~~~\___mandelbr.ps___/~~~ ~~ ~ ' Множество Мандельброта M -- множество точек c комплексной плоскости, для которых последовательность z_{n+1} = z_n^2 + c сходится. Эта программа :) рисует не само M, а его потенциал (кажется). PS: Если кто не понял -- это PostScript... |V|uxau/\ --- -- - ъ * Origin: Ф И З Ф А К - Ч Е М П И О H ! (2:5030/1072.2) Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
Архивное /ru.algorithms/38623b098c7b.html, оценка из 5, голосов 10
|