Aizawa chaotic system

Revision:
0:18b2474ec07b
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/main.cpp	Wed Mar 09 21:31:16 2011 +0000
@@ -0,0 +1,50 @@
+#include "TFT_4DGL.h"
+
+TFT_4DGL lcd(p9,p10,p11);
+
+int main() {
+
+     int   xp,yp;
+     
+     int n = WHITE;
+   
+     float x,y,z;
+   
+     float a = 0.95;
+     float b = 0.7;
+     float c = 0.6;
+     float d = 3.5;
+     float e = 0.25;
+     float f = 0.1;
+     
+     float dt = 0.01;
+ 
+     float x1 = 0.1;
+     float y1 = 0.1;
+     float z1 = 0.1;
+
+     
+     x = x1;
+     y = y1;
+     z = z1;
+
+    while(1)
+    
+    {
+    
+        x = x1+((z1-b)*x1-d*y1)*dt;
+        y = y1+(d*x1+(z1-b)*y1)*dt;
+        z = z1+(c+a*z1-(pow(z1,3)/3)-(pow(x1,2)+pow(y1,2))*(1+e*z1)+f*z1*pow(x1,3))*dt;
+ 
+        xp = (int)floor(x*95);
+        yp = (int)floor(z*85);
+
+        lcd.pixel(65+yp,160+xp,n);
+     
+        x1 = x;
+        y1 = y;
+        z1 = z;
+
+    }
+ 
+}
\ No newline at end of file