У меня проблема с решением ODE с использованием события, и мне нужна помощь.
ODE прост: (y1 '(t), y2' (t)) = (-exp (-t), exp (-t)).
и я определил функцию следующим образом:
function [Y1 Y2]= myftn(y1 y2 t)
Y1=-exp(-t);
Y2=exp(-t);
end
Исходные данные: y1 (0) = 1, y2 (0) =-1, поэтому z: = y1 ^ 2 + y2 ^ 2 сходится к 0.
Я хочу остановить численное интегрирование, когда | z (t_n) -z (t_ {n-1}) | <0.01где </p>
{t_n} - последовательность времени.
Заранее спасибо.