Brusselator chemical oscillator

Dependencies:   mbed

Revision:
0:7db0029f15d6
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/main.cpp	Fri Dec 31 16:54:57 2010 +0000
@@ -0,0 +1,41 @@
+#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.4;
+    float A = 0.05;
+    float b = 1.2;
+    float c = 0.9;
+
+    x1 = 0.3;
+    y1 = 2;
+    z1 = 0;
+    
+    dt = 0.02;
+
+    while(1)
+    
+    {
+    
+        x = x1+(pow(x1,2)*y1-(b+1)*x1+a+A*sin(c*z1))*dt;
+        y = y1+(-pow(x1,2)*y1+b*x1)*dt;
+        z = z1+(c)*dt;
+        t = t+dt;
+
+        xp = -70+(int)floor(395*x);
+        yp = -190+floor(125*y);
+
+        lcd.pixel(yp,xp,WHITE);
+      
+        x1 = x;
+        y1 = y;
+        z1 = z;
+
+    }
+ 
+}
\ No newline at end of file