Rabinovich 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,t,dt;
00008     int   xp,yp;
00009    
00010     float a = 0.84;
00011     float b = 1.1; 
00012 
00013     x1 = 1.1881;
00014     y1 = 0;
00015     z1 = 0.4;
00016     
00017     dt = 0.005;
00018 
00019     while(1)
00020     
00021     {
00022     
00023         x = x1+(y1*(z1-1+pow(x1,2))+a*x1)*dt;
00024         y = y1+(x1*(3*z1+1-pow(x1,2))+a*y1)*dt;
00025         z = z1+(-2*z1*(b+x1*y1))*dt;
00026         t = t+dt;
00027 
00028         xp = -25+(int)floor(135*x);
00029         yp = 270+(int)floor(105*y);
00030 
00031         lcd.pixel(xp,yp,WHITE);
00032       
00033         x1 = x;
00034         y1 = y;
00035         z1 = z;
00036 
00037     }
00038  
00039 }