events.m

来自「matlab宝典,电子工业出版社,包含该书的源码2-4章」· M 代码 · 共 7 行

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