# My projects library #4! Gallery Megathread

#61 THE MANDELBROT SET!!!

It took me a lot of tries to complete. I didn’t want to watch the tutorials, since I wouldn’t understand the process. So after watching tons of youtube videos explaining the mandelbrot set, I had finally done it!

Thank you to Crystalize on discord, for showing me a neat way to create complex numbers in processing.
I made all the basic functions I needed with them and commented out the rest. So right now, there are:
addition, subtraction, multiplication ( C * C or C * R ), division ( C / C or C / R), squaring and z^2+c. If you want to use them, feel free!

Code
``````complex c = new complex(0, 0);
float scale = 150;
int iter = 1000;
float cx = 300, cy = 300;
color clr[] = new color[iter];
color bg = #FF0000;
void setup() {
size(600, 600);
colorMode(HSB);
noStroke();

float m = 0.95;
float v = 1;

for (int i = 0; i < iter; i++) { //can be replaced with: 255 * pow(m,i) and shorten to a single row
clr[i] = color((v)*255, 255, 255);
v*=m;
}
clr[iter-1] = #000000; //mandelbrot set is black in the center
}

void draw() {
for (float i = 0; i < width; i ++) for (float j = 0; j < height; j++) {
if (dist(i, j, width*0.5, height*0.5)<scale*2) {
complex c = new complex( (i-300.0)/scale, (j-300.0)/scale); //og
complex d = new complex( (i-300.0)/scale, (j-300.0)/scale); //changes

int it;
for (it = 0; it < iter-1 && inBounds(d) == true; it++) {
complex nw = d.z2pc(d, c);
d.setValueC( nw );
}
pixels[(int)i+(int)j*width] = clr[it];
} else pixels[(int)i+(int)j*width] = bg;
}
updatePixels();
//if(mousePressed) {
//  println(((float)mouseX-300.0)/scale, ((float)mouseY-300.0)/scale);
//}
noLoop();
}

class complex {
float re, im;
complex(float re, float im) {
this.re = re;
this.im = im;
}
void setValue(float re, float im) {
this.re = re;
this.im = im;
}
void setValueC(complex c) {
re = c.re;
im = c.im;
}
return new complex(re+z.re, im+z.im);
}
//complex sub(complex z) {
//  return new complex(re-z.re, im-z.im);
//}
complex multC(complex z) {
return new complex(re*z.re - im*z.im, re*z.im + z.re * im);
}
//complex multR(float z) {
//  return new complex(re*z, im*z);
//}
complex devC(complex z) {
return new complex( (re*z.re + im*z.im)/(z.re*z.re+z.im*z.im), (im * z.re - re * z.im)/(z.re*z.re+z.im*z.im) );
}
complex squ() {
return multC(this);
}
complex z2pc(complex z, complex c) { //f(z) = z^2 + c