Welcome to Fractal Forums

z=quaternion(sx,sy,sz,0);
if(orbitincreaser==true)
{
z=z+5;
}
}
bool bailout(void)
return(sx*sx+sy*sy+sz*sz<bailout);

...
z=quaternion(sx,sy,sz,0);
if(orbitincreaser==true)
{
z= z*1.618033988749894848204586834365;
}
}
bool bailout(void)
return( |z| <= bailout);

NumberSeeker {
// By Edgars Malinovskis 20.12.2011
// -Z seeker marks respective negative escaping z values with dots (zeros)
// around 2D mandelbrots and on 3D mandelbulb 
// Colours positive and negative quaternion hemispheres differently
// Colour method based on Exponenthial Smoothing developed by Ron Barnett
// and orbit colouring for quaternions. Zoom requires iteration increase.
// Good with mandel bulbs and brots but not so with sophisticated formulas.

double sum;
double sum2;
double cabsz;
complex zold;
parameter bool diverge;
parameter bool converge;
parameter bool modulusz;
parameter complex posneg;
parameter int iternum;
int cnt, iIter;

	void init(void)
	{
		sum= 0.0;
		sum2= 0.0;
		zold= (0,0);
		cabsz=0;
		
		cnt=0;
		iIter=iternum;
		if (iIter>=maxiter)
		{
			iIter=maxiter-1;
		}
	}
	void loop(void)
	{
	 
if (cnt <= iIter) //calculates colours just of first iterations
{
		if  (diverge)
		{
			sum= sum+ exp(-cabs(z+posneg));
		}
		if  (converge)
		{
			sum2= sum2+ exp(-1/cabs(zold-z+posneg));
		}
		zold= z;
		cnt=cnt+1;
				if (modulusz)  // if nothing works this gives result
				{
				cabsz= cabs(z+posneg);
				}
}
	}
	void final(void)
	{
		if  (|z- zold| < 0.5)  // bounded
		{
			if  (converge)
			{
				index= sum2;
			}
			else
			{
				index= cabsz;
			}
    
		}
		else
		{
			// to infinity
			if  (diverge)
			{
				index= sum;
			}
			else
			{

				 index= cabsz;
			}
		}
	}
	void description(void)
	{
		this.title = "NumberSeeker";
 
		diverge.caption = "Exp Smooth toInfinity";
		diverge.default =  true;
		diverge.hint = "If set, points which escape to infinity will be colored.";
 
		converge.caption = "Exp Smooth toBounded";
		converge.default =  true;
		converge.hint = "If set, points which collapse to set will be colored.";
		
		modulusz.caption = "Colouring by modulus of Z";
		modulusz.default =  true;
		modulusz.hint = "Use simple modulusZ method";
 
		posneg.caption = "-Z Seeker (r,i)";
		posneg.default =  (2.0, 0.0);
		posneg.hint = "Adding number diverses negative Z areas";
		
		iternum.caption="Use first N iterations";
		iternum.default=15;
		iternum.min=0;
		iternum.hint="Lets you specify the iteration cycle at which you want to apply coloring";
	}
}