Welcome to Fractal Forums

s1=sqrt(2.0)/2.0;

 while ((norm<8)&&(iter<imax)) {
		
   x=fabs(x);
   y=fabs(y);
   z=fabs(z);
   x2=x;y2=y;z2=z;
   
   if (x>2) {norm=1000;}
   if (y>2) {norm=1000;}
   if (z>2) {norm=1000;}
   
   if ((x2<y2)&&(x2<=z2)) {d=1;}
   if ((y2<z2)&&(y2<=x2)) {d=2;}
   if ((z2<x2)&&(z2<=y2)) {d=3;}
   
   if (d==1) {
     if (x2<thresh){iter=imax+1;maxnorm=0.3;} else
     { if (y>0) { y=(y-1); }
        if (z>0) { z=z-1; }
        y2 = s1*y-s1*z;
        z2 = s1*z+s1*y;
        y=y2;z=z2;
     }
   }
   if (d==2) {
     if (y2<thresh){iter=imax+1;maxnorm=0.6;} else
     { if (x>0) { x=(x-1); }
        if (z>0) { z=z-1; }
        x2 = s1*x-s1*z;
        z2 = s1*z+s1*x;
        x=x2;z=z2;
     }	
   }
   if (d==3) {
     if (z2<thresh){iter=imax+1;maxnorm=0.9;} else
     { if (x>0) { x=(x-1); }
        if (y>0) { y=y-1; }
        x2 = s1*x-s1*y;
        y2 = s1*y+s1*x;
        x=x2;y=y2;
     }	
   }
   (*nfunc)++;
   iter++;
   r = x*x+y*y+z*z;
   x = x*2.5;y=y*2.5;z=z*2.5;
  } 
  if (iter==imax){iter=1;}
  
  if (norm>10) {iter=1;}