2D-Circular - Circular limit cycle - 2D nonlinear attractor

Users » JLS » Code » 2D-Circular

Circular limit cycle - 2D nonlinear attractor

main.cpp@0:ba94dffa2c21, 2011-02-08 (annotated)

User	Revision	Line number	New contents of line
JLS	0:ba94dffa2c21	1	#include "TFT_4DGL.h"
JLS	0:ba94dffa2c21	2
JLS	0:ba94dffa2c21	3	TFT_4DGL lcd(p9,p10,p11);
JLS	0:ba94dffa2c21	4
JLS	0:ba94dffa2c21	5	int main() {
JLS	0:ba94dffa2c21	6
JLS	0:ba94dffa2c21	7	int xp,yp;
JLS	0:ba94dffa2c21	8
JLS	0:ba94dffa2c21	9	int n = WHITE;
JLS	0:ba94dffa2c21	10
JLS	0:ba94dffa2c21	11	float x,y,xn,yn;
JLS	0:ba94dffa2c21	12
JLS	0:ba94dffa2c21	13	float dt = 0.814;
JLS	0:ba94dffa2c21	14
JLS	0:ba94dffa2c21	15	float x1 = 0.1;
JLS	0:ba94dffa2c21	16	float y1 = 0;
JLS	0:ba94dffa2c21	17
JLS	0:ba94dffa2c21	18	x = x1;
JLS	0:ba94dffa2c21	19	y = y1;
JLS	0:ba94dffa2c21	20
JLS	0:ba94dffa2c21	21	while(1)
JLS	0:ba94dffa2c21	22
JLS	0:ba94dffa2c21	23	{
JLS	0:ba94dffa2c21	24
JLS	0:ba94dffa2c21	25	xn = y1;
JLS	0:ba94dffa2c21	26	yn = (1-(x1x1)-(y1y1))*y1-x1;
JLS	0:ba94dffa2c21	27
JLS	0:ba94dffa2c21	28	x = x+xn*dt;
JLS	0:ba94dffa2c21	29	y = y+yn*dt;
JLS	0:ba94dffa2c21	30
JLS	0:ba94dffa2c21	31	xp = (int)floor(x*65);
JLS	0:ba94dffa2c21	32	yp = (int)floor(y*110);
JLS	0:ba94dffa2c21	33
JLS	0:ba94dffa2c21	34	lcd.pixel(120+xp,160+yp,n);
JLS	0:ba94dffa2c21	35
JLS	0:ba94dffa2c21	36	x1 = x;
JLS	0:ba94dffa2c21	37	y1 = y;
JLS	0:ba94dffa2c21	38
JLS	0:ba94dffa2c21	39	}
JLS	0:ba94dffa2c21	40
JLS	0:ba94dffa2c21	41	}