Welcome to Fractal Forums

	float r=length(p);
	float r2=r*r;

	float r2=dot(p, p);
	float r=sqrt(r2);

kl = new V2D(1.95, 0.03);
			b_a = kl.y / kl.x;
			box = new V2D(1.0);
			extraIter = 5;
			balls = true;
			fourGen = false;
			clamp = new V2D(0.3, 2.0);
			doInv = false;
			invCenter = new V3D(1, 0.96, 0);
			invR = 1;
			recenter = new V3D(0);

void DE() {
			V3D p = previous.ergebnis;
			V3D q = p.copy();
			V3D lz = q.getAdd(V3D.V1);
			V3D llz = q.getSub(V3D.V1);
			DF = 1.0;
			double DE = 1e10;
			double f = Math.signum(kl.y);
			for(int i = 0 ; i < maxIterM ; i++){
				q.x = q.x + b_a * q.y;
				if(fourGen){
					q.x = wrap(q.x, 2.0 * box.x, -box.x);
					q.x = wrap(q.y, kl.x, 0);
					q.x = wrap(q.z, 2.0 * box.y, -box.y);
				}else{
					q.x = wrap(q.x, 2.0 * box.x, -box.x);
					q.x = wrap(q.z, 2.0 * box.y, -box.y);
				}
				q.x = q.x - b_a * q.y;
				if(q.y >= kl.x * (0.5 + f * (2.0 * kl.x - 1.95) * 0.25 * Math.signum(q.x + kl.y * 0.5) * (1.0 - Math
							.exp(-(7.2 - (1.95 - kl.x) * 15.0) * Math.abs(q.x + kl.y * 0.5))))){
					q.x = -kl.y - q.x;
					q.y = kl.x - q.y;
					q.z = -q.z;
				}
				q = TransA(q);
				V3D tmp = q.getSub(llz);
				if(tmp.dot(tmp) < 1e-5){
					break;
				}
				llz = lz;
				lz = q;
			}
			double y = balls ? Math.min(q.y, kl.x - q.y) : q.y;
			DE = Math.min(DE, Math.min(y, clamp.x) / Math.max(DF, clamp.y));
			dist = DE;
			hitMat = mat;
		}

		private V3D TransA(V3D q) {
			double r = 1.0 / q.dot(q);
			q.mult(-r);
			q.x = -kl.y - q.x;
			q.y = kl.x + q.y;
			DF *= r;
			return q;
		}

		private double wrap(double x, double a, double s) {
			x -= s;
			return x - a * Math.floor(x / a) + s;
		}

	saveOTColor(lzy/(i+1)); // good enuff :)

public override bool GetBool(double x, double y, double z)
{
  double KleinR=1.94;
  double KleinI=0.033;
  for (int i = 0;i < 44;  i++)
  {
    x = x-Math.Floor(x+0.5);
    z = z-Math.Floor(z+0.5);
    if ( y >= KleinR*(0.5+0.1*Math.Sin(-2.0*Math.PI*x)) )
    {
      x = KleinI-x;
      y = KleinR-y;
      z = -z;
    }
    double r=x*x+y*y+z*z;
    if(r>0)
    {
      x /= r;
      y /= r;
      z /= r;
      if(r<0.2) return false;
    }
    x = KleinI + x;
    y = KleinR - y;
    z = -z;
  }
  return true;
}

public override bool GetBool(double x, double y, double z)
{
  double KleinR=1.94;
  double KleinI=0.033;
  for (int i = 0;i < 44;  i++)
  {
    x = x-Math.Floor(x+0.5);
    z = z-Math.Floor(z+0.5);
    double r=x*x+y*y+z*z;
    if(r>0)
    {
      x /= r;
      y /= r;
      z /= r;
      if(r<0.2) return false;
    }
    x = KleinI + x;
    y = KleinR - y;
    z = -z;
  }
  return true;
}

public override bool GetBool(double x, double y, double z)
{
  double KleinR=1.94;
  double KleinI=0.033;
  for (int i = 0;i < 44;  i++)
  {
    x=x%2;
    z=z%2;
    while(Math.Abs(x)>0.5)
    {
      if(x>0.5) x = 1-x;
      if(x<-0.5) x = -1-x;
    }
    while(Math.Abs(z)>0.5)
    {
      if(z>0.5) z = 1-z;
      if(z<-0.5) z = -1-z;
    }
    if ( y >= KleinR*(0.5+0.1*Math.Sin(-2.0*Math.PI*x)) )
    {
      x = KleinI-x;
      y = KleinR-y;
      z = -z;
    }
    double r=x*x+y*y+z*z;
    if(r>0)
    {
      x /= r;
      y /= r;
      z /= r;
      if(r<0.2) return false;
    }
    x = KleinI + x;
    y = KleinR - y;
    z = -z;
  }
  return true;
}

public override bool GetBool(double x, double y, double z)
{
  double KleinR=1.94;
  double KleinI=0.033;
  for (int i = 0;i < 44;  i++)
  {
    x=x%2;
    z=z%2;
    while(Math.Abs(x)>0.5)
    {
      if(x>0.5) x = 1-x;
      if(x<-0.5) x = -1-x;
    }
    while(Math.Abs(z)>0.5)
    {
      if(z>0.5) z = 1-z;
      if(z<-0.5) z = -1-z;
    }
    double r=x*x+y*y+z*z;
    if(r>0)
    {
      x /= r;
      y /= r;
      z /= r;
      if(r<0.2) return false;
    }
    x = KleinI + x;
    y = KleinR - y;
    z = -z;
  }
  return true;
}

float  JosKleinian(vec3 z)
{...}

dist = z.getLength() - r;