Encontrar los máximos de una función usando ODE45

Estoy tratando de ubicar las ubicaciones de una de las ecuaciones en un sistema de ecuaciones diferenciales en MATLAB. Estoy tratando de usar la propiedad de Eventos de odeset. ¿Cómo selecciono la ecuación particular en mi función?

options = odeset('Events',@event);
[t x tm xm ie] = ode45(@Lorenz,[0 200],I,options);


function X = Lorenz(t,x)
r = 15;
sigma = 10;
b = 8/3;
X(1,1) = sigma*(x(2,1)-x(1,1));
X(2,1) = r*(x(1,1)) - x(2,1) -x(1,1)*x(3,1);
X(3,1) = x(1,1)*x(2,1) - b*x(3,1);
end

function [value,isterminal,direction] = event(t,x)
value  = Lorenz(t,x); %I can't specify X(3,1) here
isterminal = 0;
direction = -1;
end

En particular, estoy tratando de grabar cada vez que X (3,1) = 0.

Gracias

Respuestas a la pregunta(2)

Su respuesta a la pregunta