global m l g;
m = 1;
l = 1;

g = -9.81;

y0 = [	pi/4;
		pi/4;
		0;
		0; ];

dt = 5e-4;
T = 10;

Y = euler(@delta, y0, dt, T);
plot(0:dt:T, Y(1:2,:))