Ueda chaotic oscillator

Dependencies:   mbed

Committer:
JLS
Date:
Sat Jan 01 11:34:23 2011 +0000
Revision:
0:46b738bd46bb

        

Who changed what in which revision?

UserRevisionLine numberNew contents of line
JLS 0:46b738bd46bb 1 #include "TFT_4DGL.h"
JLS 0:46b738bd46bb 2
JLS 0:46b738bd46bb 3 TFT_4DGL lcd(p9,p10,p11);
JLS 0:46b738bd46bb 4
JLS 0:46b738bd46bb 5 int main() {
JLS 0:46b738bd46bb 6
JLS 0:46b738bd46bb 7 double x1,y1,z1,x,y,z,t,dt;
JLS 0:46b738bd46bb 8 int xp,yp;
JLS 0:46b738bd46bb 9
JLS 0:46b738bd46bb 10 float a = 7.5;
JLS 0:46b738bd46bb 11 float b = 0.05;
JLS 0:46b738bd46bb 12 float c = 1;
JLS 0:46b738bd46bb 13
JLS 0:46b738bd46bb 14 x1 = 2.5;
JLS 0:46b738bd46bb 15 y1 = 0;
JLS 0:46b738bd46bb 16 z1 = 0;
JLS 0:46b738bd46bb 17
JLS 0:46b738bd46bb 18 dt = 0.005;
JLS 0:46b738bd46bb 19
JLS 0:46b738bd46bb 20 while(1)
JLS 0:46b738bd46bb 21
JLS 0:46b738bd46bb 22 {
JLS 0:46b738bd46bb 23
JLS 0:46b738bd46bb 24 x = x1+(y1)*dt;
JLS 0:46b738bd46bb 25 y = y1+(-pow(x1,3)-b*y1+a*sin(c*z1))*dt;
JLS 0:46b738bd46bb 26 z = z1+(c)*dt;
JLS 0:46b738bd46bb 27 t = t+dt;
JLS 0:46b738bd46bb 28
JLS 0:46b738bd46bb 29 xp = 165+(int)floor(40*x);
JLS 0:46b738bd46bb 30 yp = 120+(int)floor(15*y);
JLS 0:46b738bd46bb 31
JLS 0:46b738bd46bb 32 lcd.pixel(yp,xp,WHITE);
JLS 0:46b738bd46bb 33
JLS 0:46b738bd46bb 34 x1 = x;
JLS 0:46b738bd46bb 35 y1 = y;
JLS 0:46b738bd46bb 36 z1 = z;
JLS 0:46b738bd46bb 37
JLS 0:46b738bd46bb 38 }
JLS 0:46b738bd46bb 39
JLS 0:46b738bd46bb 40 }