OpenGL 地球，月亮，太阳 相对运动

/*
file:earch_moon_sun.c
每个月，月亮绕地球一圈，每年，地球绕太阳一圈。

为简单起见，假设每个月30天，一年共30x12 = 360 天。
假设天体都是球形，处于同一平面，太阳不动。

         单位 KM  放大倍数        
太阳半径 = 696000*100= 69600000
地球半径 = 6378*2500 = 15945000
月球半径 = 1738*5000 = 4345000
地月距离 = 380000*100
地日距离 = 150000000*1

把太阳放在原点。

观察点=(0,-200000000,200000000)

透视，最近可视距离为1，最远可视距离为20000000*2，
可视角75度，宽高比为1

设今天为第day天，则地球转过的角度为day/一年的天数*360度 = day

对于月球：
glRotatef(月球旋转的角度,0,0,-1);
glTranslatef(月球轨道半径,0,0);
glutSolidSphere(月球半径,20,20);

月球旋转地角度 = day/一个月的天数*360度 - 地球转过的角度

注记：若按上面的数据填入，则太阳好像破裂了；下面的长度数据
都除以1000,使得太阳正常显示。另外，月球时快时慢 没有处理好
相对运动。

*/
#include <GL/glut.h>
static int day = 200;

void display()
{
   glEnable(GL_DEPTH_TEST);
   glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
   glMatrixMode(GL_PROJECTION);
   glLoadIdentity();
   gluPerspective(75,1,1,400000);

   glMatrixMode(GL_MODELVIEW);
   glLoadIdentity();
   gluLookAt(0,-200000,200000, 0,0,0, 0,0,1);

   glColor3f(1,0,0);//sun
   glutSolidSphere(69600,20,20);

   glColor3f(0,0,1);
   glRotatef(day,0,0,-1);
   glTranslatef(150000,0,0);//earth
   glutSolidSphere(15945,20,20);

   glColor3f(1,1,0);
   glRotatef(day/30.0*360 - day, 0,0,-1);
   glTranslatef(38000,0,0);//moon
   glutSolidSphere(4345,20,20);

   glutSwapBuffers();

}

void timer(int p)
{
   day ++;
   if(day >360)day = 0;
   glutTimerFunc(50,timer,0);
   glutPostRedisplay();
}

int main(int argc,char **argv)
{
   glutInit(&argc,argv);
   glutInitDisplayMode(GLUT_RGB | GLUT_DOUBLE);
   glutCreateWindow("earth,moon,sun");
   glutInitWindowSize(400,400);

   glutDisplayFunc(display);
   glutTimerFunc(50,timer,0);
   glutMainLoop();
   return 0;
}