Welcome to Fractal Forums

	r1=sqrt (sqr(sy) + sqr(sz));  // sx, sy, and sz are the starting x, y, and z components of the iteration
	r3=r1^(-n);
			
	victor=complex(sx,r1);     // complex (variable_1, variable_2) converts 2 reals into a complex number 
	bravo=complex(sy,sz);     //  z= variable_1 + [i]i[/i] * variable_2
	victor=victor^n;  // you could write a script, so if n is 4, do sqr(sqr(victor)), if n=6 do sqr(victor)*sqr(sqr(victor)) ... etc.  
	bravo=bravo^n;
			
	nx=real (victor);                            // real grabs the real component of a complex number 
	ny=imag (victor)*real (bravo) * r3;   //  imag grabs the imaginary component of a complex number
	nz=imag (victor)*real (bravo) * r3;   //  nx, ny, and nz are the intermediate x, y, and z components of the iteration
			
	if (juliaMode) {
		sx=nx+cr;
		sy=ny+ci;
		sz=nz+cj;
				
	} else {
		sx=nx+(pixelr);
		sy=ny+abs(pixeli);    //  because this is the "Christmas tree fractal" (Name from Bugman's website) 
		sz=nz+(pixelj);         //  and it's really distorted for -y values, I take the absolute value of y... reflects it over the axis as well
	}
	z=quaternion(sx,sy,sz,0);   // for color assignment, nothing to do with the calculation of the fractal
	bail=abs(sx)+abs(sy)+abs(sz);    // bailout is absolute value, avoid squaring numbers...