Hadley chaotic circulation

Embed: (wiki syntax)

« Back to documentation index

Show/hide line numbers main.cpp Source File

main.cpp

00001 #include "TFT_4DGL.h"
00002 
00003 TFT_4DGL lcd(p9,p10,p11);
00004 
00005 int main() {
00006 
00007     double x1,y1,z1,x,y,z,t,dt;
00008     int   xp,yp;
00009    
00010     float a = 0.2;
00011     float b = 4;
00012     float f = 8;
00013     float g = 1;
00014 
00015     x1 = 0;
00016     y1 = 0;
00017     z1 = 1.3;
00018     
00019     dt = 0.005;
00020 
00021     while(1)
00022     
00023     {
00024     
00025         x = x1+(-pow(y1,2)-pow(z1,2)-a*x1+a*f)*dt;
00026         y = y1+(x1*y1-b*x1*z1-y1+g)*dt;
00027         z = z1+(b*x1*y1+x1*z1-z1)*dt;
00028         t = t+dt;
00029 
00030         xp = 90+(int)floor(90*x);
00031         yp = 115+(int)floor(50*y);
00032 
00033         lcd.pixel(yp,xp,WHITE);
00034       
00035         x1 = x;
00036         y1 = y;
00037         z1 = z;
00038 
00039     }
00040  
00041 }