Реализация алгоритма Мандельброта

#include "windows.h"
#include <conio.h>
 
int depth=2000;
 
int main(int argc, char* argv[])
{
    float cr, ci, zr, zi, tmp;
    HWND wnd=GetConsoleWindow();
    HDC dc=GetDC(wnd);
    int k;
 
    for(int n=1; n<10000; n+=100)
    {
        for(int i=-320;i<320;i++)
        {
            ci=((float)i)/320.0;
            for(int j=-380;j<160;j++)
            {
                cr=((float)j)/240.0;
                zi=zr=0.0;
                for(k=0;k<depth;k++)    
                {
                    tmp=zr*zr-zi*zi;
                    zi=2*zr*zi+ci;
                    zr=tmp+cr;
                    if (zr*zr+zi*zi>1.0E16) break;
                }
                if (k<depth)
                    SetPixel(dc, i+320, j+380, k*255*n);
                else SetPixel(dc, i+320, j+380, 10);
            }
        }
        if(_kbhit()) break;
    }
 
    while(!_kbhit()) ;
    return 0;
}

GetConsoleWindow

#define _WIN32_WINNT 0x0501
#include "windows.h"
#include <wincon.h>
#include <wingdi.h>
#include <conio.h>