Kamil Ondrousek
/
Ueda
Ueda chaotic oscillator
main.cpp@0:46b738bd46bb, 2011-01-01 (annotated)
- Committer:
- JLS
- Date:
- Sat Jan 01 11:34:23 2011 +0000
- Revision:
- 0:46b738bd46bb
Who changed what in which revision?
User | Revision | Line number | New 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 | } |