Lorenz-Mod2 chaotic attractor

Dependencies:   mbed

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,dt;
00008     int   xp,yp;
00009    
00010     float a = 0.9;
00011     float b = 5; 
00012     float c = 9.9;
00013     float d = 1;
00014 
00015     x1 = 0;
00016     y1 = 1;
00017     z1 = 0;
00018     
00019     dt = 0.005;
00020 
00021     while(1)
00022     
00023     {
00024     
00025         x = x1+(-a*x1-pow(y1,2)-pow(z1,2)+a*c)*dt;
00026         y = y1+(x1*(y1-b*z1)+d)*dt;
00027         z = z1+(-z1+x1*(b*y1+z1))*dt;
00028 
00029         xp = 100+(int)floor(48*x);
00030         yp = 160+(int)floor(38*y);
00031 
00032         lcd.pixel(xp,yp,WHITE);
00033       
00034         x1 = x;
00035         y1 = y;
00036         z1 = z;
00037 
00038     }
00039  
00040 }