Welcome to Fractal Forums

() 
{
enum {SIZE=512,PRECALC=100, PARTNUM=10};
static tab[SIZE][SIZE];
static extabr[10000],extabg[10000],extabb[10000];
static expo,HSIZE,QSIZE,i1,j1;

if (numframes==0){
   HSIZE=SIZE/2;
   QSIZE=SIZE/4;
   expo=1.0;expor=0.5;expog=0.5;expob=0.025;
   for (i=0; i<1000;i++){
      extabr[i]=256*(exp(-expor*i/256));
      extabg[i]=256*(exp(-expog*i/256));
      extabb[i]=256*(exp(-expob*i/256));
   }
   i1=0;j1=0;
}

al=alpha(); be=beta();

for(i=0;i<SIZE;i++)
   for(j=0;j<SIZE;j++) tab[i][j]=0;

//cls(255,255,255);
cls(255,255,255);

x=100; y=100; INFRM=-1/1000;
   for (i=0;i<PRECALC;i++)
   {
      a=x-al;
      b=y-be;
      x=sqrt(a*a+b*b);
      y=sqrt((x-a)*0.5);
      x=sqrt((x+a)*0.5);
      if (b<0) y=-y;
      a=rdm();b=rdm();
      x*=a; y*=b;//we should do y*=a instead :)
   }
   do
   {
      for (i=0;i<PARTNUM;i++)
      {
         a=x-al;
         b=y-be;
         x=sqrt(a*a+b*b);
         y=sqrt((x-a)*0.5);
         x=sqrt((x+a)*0.5);
         if (b<0) y=-y;
         a=rdm();b=rdm();
         x*=a; y*=b;//we should do y*=a instead :)
         
         //putpixel(x,y);
         i1=x*QSIZE+HSIZE;
         j1=y*QSIZE+HSIZE;
         tab[i1][j1]+=1;
         c=(1-exp(INFRM*tab[i1][j1]))*1000;
         setcol(extabr[c],extabg[c],extabb[c]);
         setpix(i1,j1);
      };
   }while((al==alpha())&&(be==beta()));
}

rdm()
{
   a=RND;
   if (a>0.5) return 1;
   return -1;
}

alpha()
{
   return 4*(mousx-0.5*xres)/xres;
}

beta()
{
   return 4*(mousy-0.5*yres)/xres;
}

putpixel(x,y)
{
   setpix(x*0.25*xres+0.5*xres,y*0.25*xres+0.5*yres);
}