function f=pendulum(t,x)

L=1; W=1;g=32.2; I=L^2*W/(3*g);

f(1)=-W*(L/2)*sin(x(2));
f(2)=x(1)/I;

f=f';