pro7.c

来自「TC 下的物理程序」· C语言 代码 · 共 46 行

/*洛伦兹吸引子*/
#include<d:\ccbp\ccbp.h>
 float dt = 0.001;
 
 int main()
 {
    float a,b,c,x,y,z,x1,y1,z1;
    int i,color;
    int gdriver = DETECT,gmode;
    initgraph(&gdriver,&gmode,"");
    setbkcolor(1);
 
     setbkcolor(5);
     a = 10; b=2.666;
     x = y = z=2;
     for(c = 0;c < 40;c += 1){
        outhzxy3(400,360,"洛仑兹引力",2,0,15);
 
     line(10,10,10,250);line(10,250,300,250);
     line(310,250,310,10);line(310,250,600,250);
     line(50,450,50,280);line(50,450,400,450);
     
     outtextxy(420,410,"press any key to continue ...");
     printf("the parameter c is %f\r",c);
     outtextxy(20,20,"x");outtextxy(290,260,"y");
     outtextxy(320,10,"x");outtextxy(600,260,"z");
     outtextxy(60,280,"y");outtextxy(400,460,"z");
     
     for(i=1;i<=22000;i++){
        x1=-a*(x-y);
        y1=-y-x*z+c*x;
        z1=-b*z+x*y;
        x=x+x1*dt;
        y=y+y1*dt;
        z=z+z1*dt;
        putpixel(170+5*x,150+2.5*y,(int)(i/15));
        putpixel(470+5*x,20+2.5*z,(int)(i/14));
        putpixel(210+5*y,300+2.5*z,(int)(i/13));
     }
     getch();
     cleardevice();
   }
   getch();
   return 0;
}