Фрактал Мандельброта

#include <fstream>
#include <iostream>
 
using namespace std;
 
 
 
 
int findMandelbrot(double cr, double ci, double max_iterations) {
    int i=0;
    double zr = 0.0, zi = 0.0;
    while (i < max_iterations && zr * zr + zi *zi < 4.0)
 
    {
        double temp = zr * zr - zi * zi + cr;
        zi = 2.0 * zr * zi  + ci;
        zr = temp;
        i++;
    }
    return i;
 
}
 
double mapToReal(int x, int imageWidth, double minR, double maxR)
{
    double range = maxR - minR;
    return x * ( range / imageWidth ) + minR;
 
}
double mapToImaginary(int y, int imageHeight, double minI, double maxI)
{
    double range = maxI - minI;
    return y * (range / imageHeight ) + minI;
 
}
 
 
int main() {
    ifstream fin ("input.txt"); // берёт данные из файла
    int imageWidth, imageHeight, maxN;
    double minR, maxR, minI, maxI;
    if (!fin)
    {
        cout << " Could not open file " << endl;
        cin.ignore();
        return 0;
    }
    fin >> imageWidth >> imageHeight >> maxN;
    fin >> minR >> maxR >> minI >> maxI;
    fin.close();
 
    ofstream fout ("output_image.ppm"); //открывает файл
    fout << "P3" << endl;
    fout << imageWidth << " " << imageHeight << endl; //размеры
    fout << "256" << endl;  //максимально число пикслелей в ргб
 
    for (int y = 0; y < imageHeight; y++) //строки
    {   
        for (int x = 0; x < imageWidth; x++) //пиксели в строке
        {
 
        double cr = mapToReal (x , imageWidth , minR ,maxR);
        double ci = mapToImaginary (y, imageHeight, minI, maxI);
 
        int n = findMandelbrot ( cr , ci , maxN);
 
        int r = ((int)(n* sinf(n)) % 256);
        int g = ((n * 10) % 256);
        int b = (n % 256);
 
        fout << r << "  " << g << "  " << b << "  ";
 
        }
        fout << endl;
 
        
 
    }
    fout.close();
    cout << "Finished " << endl;
    return 0;
 
 
 
}

512 512 256 -1.5 0.7 -1.0 1.0
 
image Width, image Height
maximum value of N in z_n (mandelbrot formula)
minReal, maxReal
minImaginary, maxImaginary.