Struggling at 'Events' of ODE45

1 Ansicht (letzte 30 Tage)
Felix Lauwaert
Felix Lauwaert am 30 Nov. 2015
Kommentiert: Felix Lauwaert am 1 Dez. 2015
Hello,
I have a ODE and need to compute the "exact" time and values of y at y(2)=0 and I'm using ode45 and 'Events'.
When I run my function, I get the error message: 'error using events. Not enough input variables.'
The part of my code:
options=['RelTol',1e2*AbsTol,'AbsTol',AbsTol*ones(1,length(condIni))...
,'Refine',8,'Events',seccio];
[T,Y,Te,Ye,Ie] = ode45(funRTBP,incT,condIni,options );
And my events function:
function [value,isterminal,direction] = seccio(t,y)
value = y(2);
direction = 0;
isterminal = 1;
end
AbsTol and condIni are correctly defined.
Thanks for your help.

Akzeptierte Antwort

Walter Roberson
Walter Roberson am 30 Nov. 2015
options = odeset('RelTol', 1e2*AbsTol, 'AbsTol', AbsTol*ones(1,length(condIni)), ...
'Refine', 8, 'Events', @seccio);
  3 Kommentare
Walter Roberson
Walter Roberson am 1 Dez. 2015
I think I will need to see more of your code.
Did you notice that I changed from [] around the options to odeset() around the options?
Felix Lauwaert
Felix Lauwaert am 1 Dez. 2015
I hadn't seen it, now it works! Thank you vey much.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Produkte

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by