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