function xdot=myfunc4(t,x)
 xdot=zeros(2,1); % —ÌˆæŠm•Û
 a=0.5;
 xdot(1)=-a*x(1)-sin(x(2));
 xdot(2)=x(1);
% end of function
