Поправьте код.

// Моделирование сил притяжения
#include <cmath>
#include <GL/glut.h>
#pragma comment (lib, "OpenGL32.lib")
#pragma comment (lib, "Glu32.lib")
#pragma comment (lib, "Glaux.lib")
struct Particle
{ 
        float x;
        float y;
        float vx;
        float vy;
        float m;
};
 
const int N = 10;
 
Particle particles[N];
 
void display()
{
        glClear(GL_COLOR_BUFFER_BIT);
        glBegin(GL_POINTS);
        for (int i = 0; i < N; ++i)
                glVertex2f(particles[i].x, particles[i].y);
        glEnd();
 
        glutSwapBuffers();
}
void timer(int=0)
{
        for (int i = 0; i < N; ++i)//;
        {
                Particle &p0 = particles[i];
                for (int j = 0; j < N; ++j)//;
                {
                        if (j == i)
                                continue;
                        const Particle &p = particles[j];
                        float d = sqrt((p0.x - p.x) * (p0.x - p.x) + 
                                (p0.y - p.y) * (p0.y - p.y));
                        if (d > 2000)
                        {
                              p0.vx +=  0,1* p.m / d / d * (p.x - p0.x) / d;
                              p0.vy +=  0,1* p.m / d / d * (p.y - p0.y) / d; 
                        }
                }
                p0.x += p0.vx;
                p0.y += p0.vy;
        }
        display();
        glutTimerFunc(1, timer, 1);
}
int main(int argc, char **argv)
{
        particles[0].x = 100;
        particles[0].y = 100;
        particles[0].vx = 0;
        particles[0].vy = 0;
        particles[0].m = 100;
 
        particles[1].x = 105;
        particles[1].y = 100;
        particles[1].vx = 0,5;
        particles[1].vy = 0,;
        particles[1].m = 10;
 
        particles[2].x = 110;
        particles[2].y = 100;
        particles[2].vx = 0,2;
        particles[2].vy = 0;
        particles[2].m = 7;
 
        particles[3].x = 115;
        particles[3].y = 100;
        particles[3].vx = 0,15;
        particles[3].vy = 0;
        particles[3].m = 5;
 
        particles[4].x = 116;
        particles[4].y = 100;
        particles[4].vx = 0,15;
        particles[4].vy = 0;
        particles[4].m = 5;
 
        particles[5].x = 120;
        particles[5].y = 100;
        particles[5].vx = 0,15;
        particles[5].vy = 0;
        particles[5].m = 5;
 
        particles[6].x = 165;
        particles[6].y = 100;
        particles[6].vx = 0,15;
        particles[6].vy = 0;
        particles[6].m = 5;
 
        particles[7].x = 170;
        particles[7].y = 100;
        particles[7].vx = 0,15;
        particles[7].vy = 0;
        particles[7].m = 5;
 
        particles[8].x = 180;
        particles[8].y = 100;
        particles[8].vx = 0,15;
        particles[8].vy = 0;
        particles[8].m = 5;
 
        particles[9].x = 190;
        particles[9].y = 100;
        particles[9].vx = 0,15;
        particles[9].vy = 0;
        particles[9].m = 5;
 
 
        glutInit(&argc, argv);
        glutInitDisplayMode(GLUT_DOUBLE | GLUT_RGB);
        glutInitWindowSize(200, 200);
        glutInitWindowPosition(20, 810);
        glutCreateWindow("Gravitation");
        
        glMatrixMode(GL_PROJECTION);
        glLoadIdentity();
        glOrtho(0, 200, 200, 0, -1, 1);
        glutDisplayFunc(display);
        timer();
        glutMainLoop();
        return 0;
}