vdpol_events.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% vdpol_events.m
% 响应事件函数
mu=2;
options=odeset('Events',@vdpolevents);
[t,y,te,ye]=ode45(@vdpol,tspan,yo,options,mu);
plot(t,y,te,ye(:,2),'o');
title('Van der Pol resutls for |y(2)=1.5|')
legend('x','dx/dt');