function [value,isterminal,direction] = eventr(t,x)

global L theta0 deltar 

delta=L*sin(x(1))-L*sin(theta0);

value = delta-deltar;

isterminal = 1;

direction = 0;