Halvorsen 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 = 1.4;
00011 
00012     x1 = -5;
00013     y1 = 0;
00014     z1 = 0;
00015     
00016     dt = 0.005;
00017 
00018     while(1)
00019     
00020     {
00021     
00022         x = x1+(-a*x1-4*y1-4*z1-pow(y1,2))*dt;
00023         y = y1+(-a*y1-4*z1-4*x1-pow(z1,2))*dt;
00024         z = z1+(-a*z1-4*x1-4*y1-pow(x1,2))*dt;
00025         t = t+dt;
00026 
00027         xp = 210+(int)floor(13*x);
00028         yp = 155+(int)floor(9*y);
00029 
00030         lcd.pixel(yp,xp,WHITE);
00031       
00032         x1 = x;
00033         y1 = y;
00034         z1 = z;
00035 
00036     }
00037  
00038 }