Моделирование атома (3D-модель)

#include <math.h>
#include <windows.h>
 
int main(){
    HWND hWnd = ::GetForegroundWindow();
    HDC  hDC  = GetDC(hWnd);
    RECT pRect= {0};
    GetClientRect(hWnd, &pRect);
    HBRUSH hBack = CreateSolidBrush(RGB(255, 255, 255));
    HBRUSH hProt = CreateSolidBrush(RGB(0, 0, 255));
    HBRUSH hElec = CreateSolidBrush(RGB(255, 0, 0));
 
    float x = 0;
    float y = 0;
    float r = 200;
    for( int i = 0; i < 512; i++ )
    {
        x = i * r / 50;
        y = sqrt(r*r - x*x);
        hBack = (HBRUSH)SelectObject(hDC, hBack);
        Rectangle(hDC, 0, 0, pRect.right, pRect.bottom);
        hBack = (HBRUSH)SelectObject(hDC, hBack);
 
        hProt = (HBRUSH)SelectObject(hDC, hProt);
        Ellipse(hDC, pRect.right / 2 - 50, pRect.bottom / 2 + 50, pRect.right / 2 + 50, pRect.bottom / 2 - 50);
        hProt = (HBRUSH)SelectObject(hDC, hProt);
 
        hElec = (HBRUSH)SelectObject(hDC, hElec);
        Ellipse(hDC, pRect.right / 2 - 5 + x, pRect.bottom / 2 + 5 + y, pRect.right / 2 + 5 + x, pRect.bottom / 2 - 5 + y);
        hElec = (HBRUSH)SelectObject(hDC, hElec);
 
        Sleep(100);
    }
    return 0;
}