Welcome to Fractal Forums

static n=2;
(x,y) {
   static cx=-1,cy=0.5;
   static bo=2;
   static frm=-1;
   if(numframes>frm){
      frm=numframes;
      cx=(mousx-0.5*xres)*0.005;
      cy=(0.5*yres-mousy)*0.005;
      bo=newton(1,sqrt(cx*cx+cy*cy))^2;
   }
   1-julia(x,y,cx,cy,32,bo)
}
julia(x,y,cx,cy,it,bo){
   r2=x*x+y*y;
   for(i=0;i<it && r2<bo;i++){
      r=sqrt(r2);
      t=atan2(y,x);
      r=r^n;t=n*t;
      x=r*cos(t)+cx;
      y=r*sin(t)+cy;
      /*x1=x*x-y*y+cx;
      y=2*x*y+cy;
      x=x1;*/
      r2=x*x+y*y;
   }
   i/it
}
newton(x,c){//Solves for the radius
   i=0;
   do{
      x1=x-(x^n-x-c)/(n*x^(n-1)-1);
      if(abs(x-x1)/x<2^-53) break;
      x=x1;
      i++;
   }while(i<5);
   x
}

static n=2;
(x,y) {
   cx=x;cy=y;
   bo=newton(1,sqrt(cx*cx+cy*cy))^2;
   1-julia(x,y,cx,cy,32,bo)//in fact this is the mandelbrot set :) replace "bo" by 2 and see the difference
}
julia(x,y,cx,cy,it,bo){
   r2=x*x+y*y;
   for(i=0;i<it && r2<bo;i++){
      r=sqrt(r2);
      t=atan2(y,x);
      r=r^n;t=n*t;
      x=r*cos(t)+cx;
      y=r*sin(t)+cy;
      /*x1=x*x-y*y+cx;
      y=2*x*y+cy;
      x=x1;*/
      r2=x*x+y*y;
   }
   i/it
}
newton(x,c){
   i=0;
   do{
      x1=x-(x^n-x-c)/(n*x^(n-1)-1);
      if(abs(x-x1)/x<2^-53) break;
      x=x1;
      i++;
   }while(i<5);
   x
}