Часть моего кода от решателя ODE:
function rdot = f(t, r)
rdot(1) = r(2);
calculation expressions
...
...
rdot(2) = some expression;
rdot(3) = r(4);
calculation expressions;
...
...
rdot(4) = some expression;
rdot = rdot';
Я хочу, чтобы значения rdot (1) и rdot (4) сохранялись для каждого входного значения t. Я нуждаюсь в них для некоторых других вычислений для двух других переменных f1 и f2 позже. Итак, я думаю, мне нужно сделать массив из двух новых переменных для хранения rdot (3) и rdot (4). Но как мне это реализовать?
Или я должен сделать
function rdot = f(t, r, f1, f2)
и производите ли вычисления для f1 и f2 в одном и том же разделе этого кода функции? Но это может запутаться с решателем оды?