For each point kon the complex plane do:
    let x=0.
    repeat infinite times:
        x=x^2+k.
    end repeat
    if x goes to infinity,
        k is not
in the set. Color is white.
    else
        k is in
the set. Color is black.

Понятно, что бесконечных циклов быть не должно. Поэтому возьмем некоторое большое число I и проитерируем I раз. Чем большее I мы взяли, тем, понятнее, точнее ответ мы получим. Из практики, число 4000 дает довольно хороший результат. Да, но 4000 раз “крутить“ цикл для каждого пиксела изображения, это многовато. К счастью, мы можем воспользоваться результатами многолетней работы математиков в этой области. Оказывается, если в любой конкретный момент вычислений, для k расстояние от z_i(k) начала координат больше 2, то мы можем принять, что данная {Z_n(k)} уйдет в бесконечность (При сравнении: расстояние < 2, поэтому его квадрат меньше 4 и корень извлекать не нужно). Итак, теперь наш алгоритм выглядит так:

For each point k in the complex plane do:
    let x=0.
    repeat 4000 times
        let x=x^2+k
        if x^2 > 4 then Color it white
Break.
    end repeat
    if we reached 4000 then
        Color
it black.

Этот метод дает нам черно-белое изображение множества Мандельброта. Теперь надо подумать о том, как сделать его разноцветным.

For each point k in the complex plane do:
    let x:=0.
    for i:=0 to 4000
        let x=x^2+k
        if ( |x|^2 > 4) then Color point k color i
Break;
        end if
    end for
    if (i=4000)
        Color
point k black.
    end if
	

Конечно, просто рисовать точку цветом i мы не можем. Считая, что у нас есть только 256 градаций серого, а i меняется до 4000. Нам надо как-то отображать i на доступный нам диапазон цветов. Эту проблему мы оставляем вам. После того, как мы получили приличное изображение в градациях серого, очень легко чуть изменить алгоритм для получения цветного изображения. Например, в изображении в градациях серого, если точка вышла из области на n-ой, вы можете рисовать ее цветом (n, n, n). Можете попробовать и что-нибудь поинтереснее типа (n, 255 – n, 50 mod n * 3). Оставляем простор для вашей фантазии. И последнее: обычно, все множество Мандельброта расположено от -2 до 0.5 по действительной оси и от –1.25 до 1.25 по мнимой оси. Ваша программа не должна тестировать точки далеко за пределами этой области.

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <dos.h>
#include <conio.h>

#define COLOR 100
#define MAS 0.9

typedef struct complex Complex;

void Sqr(Complex *z)
{
  Complex Fool=*z;

  z->x=Fool.x*Fool.x-Fool.y*Fool.y;
  z->y=2*Fool.x*Fool.y;
}

char GetColor(Complex zInit)
{
  Complex z=zInit;
  int Color=COLOR;

  while(z.x*z.x+z.y*z.y <= 4 && Color)
  {
    Sqr(&z);

    z.x+=zInit.x;
    z.y+=zInit.y;

    Color--;
  }

  return Color;
}

void DrawMandelSet(double xMin,double xMax,double yMin,double yMax)
{
  double xInc,yInc;
  Complex zInit;

  int y,x;

  char far *Screen=(char far *)MK_FP(0xa000,0);

  zInit.y=yMin;

  xInc=(xMax-xMin)/320;
  yInc=(yMax-yMin)/200;

  for(y=0;y<200;y++,zInit.y+=yInc)
  {
    zInit.x=xMin;

    for(x=0;x<320;x++,zInit.x+=xInc,Screen++)
    *Screen=GetColor(zInit);
  }
}

void main(void)
{
  _AX=0x13;geninterrupt(0x10);

  DrawMandelSet(-2*MAS,1*MAS,-1*MAS,1*MAS);

  getch();

  _AX=0x03;geninterrupt(0x10);}

Copyright 2000-2002 © Ilia Kantor, при поддержке проекта MANUAL.RU