function dy=fPendul(t,y,flag,g,L)
dy=[y(2);-g/L*sin(y(1))];