Kamil Ondrousek
/
VanDerPol
main.cpp@0:03f336a2328f, 2010-12-30 (annotated)
- Committer:
- JLS
- Date:
- Thu Dec 30 21:32:54 2010 +0000
- Revision:
- 0:03f336a2328f
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
JLS | 0:03f336a2328f | 1 | #include "TFT_4DGL.h" |
JLS | 0:03f336a2328f | 2 | |
JLS | 0:03f336a2328f | 3 | TFT_4DGL lcd(p9,p10,p11); |
JLS | 0:03f336a2328f | 4 | |
JLS | 0:03f336a2328f | 5 | int main() { |
JLS | 0:03f336a2328f | 6 | |
JLS | 0:03f336a2328f | 7 | double x1,y1,z1,x,y,z,t,dt; |
JLS | 0:03f336a2328f | 8 | int xp,yp; |
JLS | 0:03f336a2328f | 9 | |
JLS | 0:03f336a2328f | 10 | float a = 5; |
JLS | 0:03f336a2328f | 11 | float b = 3; |
JLS | 0:03f336a2328f | 12 | float c = 1.788; |
JLS | 0:03f336a2328f | 13 | |
JLS | 0:03f336a2328f | 14 | x1 = -1.9; |
JLS | 0:03f336a2328f | 15 | y1 = 0; |
JLS | 0:03f336a2328f | 16 | z1 = 0; |
JLS | 0:03f336a2328f | 17 | |
JLS | 0:03f336a2328f | 18 | dt = 0.005; |
JLS | 0:03f336a2328f | 19 | |
JLS | 0:03f336a2328f | 20 | while(1) |
JLS | 0:03f336a2328f | 21 | |
JLS | 0:03f336a2328f | 22 | { |
JLS | 0:03f336a2328f | 23 | |
JLS | 0:03f336a2328f | 24 | x = x1+(y1)*dt; |
JLS | 0:03f336a2328f | 25 | y = y1+(-x1+b*(1-pow(x1,2))*y1+a*sin(c*z1))*dt; |
JLS | 0:03f336a2328f | 26 | z = z1+(c)*dt; |
JLS | 0:03f336a2328f | 27 | t = t+dt; |
JLS | 0:03f336a2328f | 28 | |
JLS | 0:03f336a2328f | 29 | xp = 160+(int)floor(70*x); |
JLS | 0:03f336a2328f | 30 | yp = 120+(int)floor(15*y); |
JLS | 0:03f336a2328f | 31 | |
JLS | 0:03f336a2328f | 32 | lcd.pixel(yp,xp,WHITE); |
JLS | 0:03f336a2328f | 33 | |
JLS | 0:03f336a2328f | 34 | x1 = x; |
JLS | 0:03f336a2328f | 35 | y1 = y; |
JLS | 0:03f336a2328f | 36 | z1 = z; |
JLS | 0:03f336a2328f | 37 | |
JLS | 0:03f336a2328f | 38 | } |
JLS | 0:03f336a2328f | 39 | |
JLS | 0:03f336a2328f | 40 | } |