Нахождение максимумов функции с использованием ODE45

Я пытаюсь определить местонахождение одного из уравнений в системе дифференциальных уравнений в MATLAB. Я пытаюсь использовать свойства событий odeset. Как мне выбрать конкретное уравнение в моей функции?

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

В частности, я пытаюсь записывать всякий раз, когда X (3,1) = 0.

Спасибо

Ответы на вопрос(2)

Ваш ответ на вопрос