function xdot=myfunc5(t,x)
 xdot=zeros(2,1); % —ÌˆæŠm•Û
 F=0.01*sign(sin(0.967*t));
 g=1;
 a=0.5;
 m=1;
 l=1;
 xdot(1)=(-x(1)*a + (-sin(x(2))*g + F)/l )/m;
 xdot(2)=x(1);
% end of function
