function [value,isterminal,direction] = evente(t,x) % theta=x(1); % omega=x(2); global L yc theta0 delta=real(L*sin(x(1))-L*sin(theta0)); yc0=1.9*yc; value =delta-yc0; isterminal = 1; direction = 1;