How to stop ode integration?
조회 수: 1 (최근 30일)
이전 댓글 표시
I'm trying to solve a differential equations system with ode45 solver. But sometimes integration takes too much time because one of three values in the system becomes too small. As I understood I should use an event in ode function.
function [value,isterminal,direction] = events(t,y)
value = y(1) - 0.0000001;
isterminal = 1;
direction = 0;
How to write the event not just only for one y value, but for all three of them?
댓글 수: 0
채택된 답변
Walter Roberson
2014년 2월 6일
function [value,isterminal,direction] = events(t,y)
miny = [0.0000001, 0.000000025, -83.9]; %for example
value = y < miny;
isterminal = 1;
direction = 0;
댓글 수: 2
Walter Roberson
2014년 2월 6일
Ah, it should probably be
isterminal = ones(size(y));
direction = zeros(size(y));
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Ordinary Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!