function [value,isterminal,direction] = events(t,y)
% Locate the time when height passes through zero in a decreasing direction
% and stop integration.  
value = y(1);     % detect height = 0
isterminal = 1;   % stop the integration
direction = -1;   % negative direc