#include <stdio.h>
#include <limits.h>
main() {
unsigned int resx = 4096, resy = 4096;  //Resolution;
unsigned char i = 0;                    //Iteration counter;
float xmin = -2, xmax = +2;             //Borders of the area
float ymin = -2, ymax = +2;             //to be plotted;
unsigned int xpix = 0, ypix = 0;        //Column and row counters;
float x, y;                             //Coordinates;
long double Rez, Imz, Rez_new, Imz_new; //Real and imaginary parts.
FILE *f;

f = fopen("tricorn.pgm", "w");
fprintf(f, "P2\n%d %d\n%d\n", resx, resy, UCHAR_MAX); //Portable Greymap file header
for (ypix = 0; ypix < resy; ypix++) {
	for (xpix = 0; xpix < resx; xpix++) {
		x = xmin + xpix * (xmax - xmin)/resx;
		y = ymin + ypix * (ymax - ymin)/resy;
		Rez = x; Imz = y;
		for (i = 0; Rez*Rez + Imz*Imz <= 4; i++) { //z = conj(z^2) + c
			Rez_new = Rez*Rez - Imz*Imz + x;       //(It'd be z = z^2 +c for Mandelbrot)
			Imz_new = -2*Rez*Imz + y; //remove minus sign here to obtain the Mandelbrot set
			Rez = Rez_new;
			Imz = Imz_new;
			if (i == UCHAR_MAX) break; }
		fprintf(f, "%d ", i); }
	fprintf(f, "\n"); }
}