для моего кода
options=odeset('OutputFcn',@myfun)
[T,Y]=ode15s(@(T,X)sys(T,X,vin),[t0 .0005],X(:,1),options)
function Xdot = sys(T,X,vin)
Xdot=A*X+B*U
end
Я не знаю, как написать функцию myfun
для выполнения OutputFcn
для сохранения вектора решения на каждом шаге.