Kamil Ondrousek / Mbed 2 deprecated Hadley

Hadley chaotic circulation

main.cpp

Committer:
JLS
Date:
2011-01-01
Revision:
0:1427f7906a08

File content as of revision 0:1427f7906a08:

#include "TFT_4DGL.h"

TFT_4DGL lcd(p9,p10,p11);

int main() {

    double x1,y1,z1,x,y,z,t,dt;
    int   xp,yp;
   
    float a = 0.2;
    float b = 4;
    float f = 8;
    float g = 1;

    x1 = 0;
    y1 = 0;
    z1 = 1.3;
    
    dt = 0.005;

    while(1)
    
    {
    
        x = x1+(-pow(y1,2)-pow(z1,2)-a*x1+a*f)*dt;
        y = y1+(x1*y1-b*x1*z1-y1+g)*dt;
        z = z1+(b*x1*y1+x1*z1-z1)*dt;
        t = t+dt;

        xp = 90+(int)floor(90*x);
        yp = 115+(int)floor(50*y);

        lcd.pixel(yp,xp,WHITE);
      
        x1 = x;
        y1 = y;
        z1 = z;

    }
 
}