/*  旋转的立体椭球 */
　　#include graphics.h
　　#include math.h
　　#include conio.h
　　#define PI 3.14159
　　#define T PI/180
　　#define NN 36
　　#define DT PI/NN
　　#define DIST 0.8
　　#define R 190
　　static int n=1;
　　void trans(t,s)
　　float t[3];
　　float s[3][3];
　　{
　　 s[0][0]=cos(t[1])*cos(t[2]); s[0][1]=cos(t[1])*sin(t[2]);
　　 s[0][2]=-sin(t[1]);
　　 s[1][0]=sin(t[0])*sin(t[1])*cos(t[2])-cos(t[0])*sin(t[2]);
　　 s[1][1]=sin(t[0])*sin(t[1])*sin(t[2])+cos(t[0])*cos(t[2]);
　　 s[1][2]=sin(t[0])*cos(t[1]);
　　 s[2][0]=cos(t[0])*sin(t[1])*cos(t[2])+sin(t[0])*sin(t[2]);
　　 s[2][1]=cos(t[0])*sin(t[0])*sin(t[2])-sin(t[0])*cos(t[2]);
　　 s[2][2]=cos(t[0])*cos(t[1]);
　　}
　　void draw(m,da,db)
　　float m[3][3],da,db;
　　{   
　　 float f[3],f0 , x0,y0,x1,y1,x2,y2;
　　    f[0]=R*sin(da)*cos(db); f[1]=R*sin(da)*sin(db);
　　 f[2]=R*cos(da);
　　    f0=f[0]*m[0][2]+f[1]*m[1][2]+f[2]*m[2][2];
　　 if(f0=0) n=1;
　　    else
　　    {
　　  x0=300.0;y0=165.0;
　　  x2=(m[0][0]*f[0]+m[1][0]*f[1]+m[2][0]*f[2])+x0;
　　  y2=(m[0][1]*f[0]+m[1][1]*f[1]+m[2][1]*f[2])*DIST+y0;
　　  if(n==1) {n=2;x1=x2;y1=y2;}
　　  else
　　  {
　　   line(x1,y1,x2,y2); 
　　   x1=x2;y1=y2;
　　  }
　　 }
　　}main()
　　{
　　 float d[3],r[3][3],ta,tb;
　　 char k;
　　 int gdriver=VGA,gmode=VGAMED, i,p=1;
　　    initgraph(&gdriver,&gmode,"c:c");
　　 setbkcolor(BLACK);setcolor(GREEN);
　　    d[1]=30*T; d[2]=10*T;
　　do{
　　 for(i=0;i=361;i+=1)
　　 {
　　  k=kbhit();
　　  if(k!=0)break;
　　  setactivepage(p);
　　  d[0]=i*T;
　　  cleardevice();
　　        trans(d,r);
　　   for(tb=0.0;tbPI;tb+=DT)
　　   {
　　    n=1;
　　          for(ta=0.0;ta2.1*PI;ta+=DT) draw(r,ta,tb);}
　　          for(ta=0.0;taPI;ta+=DT)
　　    {   
　　    n=1;
　　    for(tb=0.0;tb2.1*PI;tb+=DT) draw(r,ta,tb);
　　    }
　　         setvisualpage(p);delay(150);p=1-p; 
　　  }
　　 }
　　    
　　  while(k==0);
　　  getch();closegraph();
　　}TC2.0调试通过，希望大家能喜欢！