function y = x2dot(x)
g=10; % gravitasi
l=0.28; % panjang pendulum
x4dot=x(1);
x3=x(2);
y = (g*tan(x3)) - ((l*x4dot)/cos(x3));