Struggling at 'Events' of ODE45

1 visualizzazione (ultimi 30 giorni)
Felix Lauwaert
Felix Lauwaert il 30 Nov 2015
Commentato: Felix Lauwaert il 1 Dic 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.

Risposta accettata

Walter Roberson
Walter Roberson il 30 Nov 2015
options = odeset('RelTol', 1e2*AbsTol, 'AbsTol', AbsTol*ones(1,length(condIni)), ...
'Refine', 8, 'Events', @seccio);
  3 Commenti
Walter Roberson
Walter Roberson il 1 Dic 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 il 1 Dic 2015
I hadn't seen it, now it works! Thank you vey much.

Accedi per commentare.

Più risposte (0)

Prodotti

Community Treasure Hunt

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

Start Hunting!

Translated by