simulate_points.m

来自「一个matlab的将军模型」· M 代码 · 共 12 行

function Pf = simulate_points(P0,sys_eq,ode_param,tf);

Pf = vertices;
for k = 1:length(P0)
  % syntax: ODE45('F',TSPAN,Y0,OPTIONS,P1,P2,...)
  [T,X] = ode45(sys_eq,[0 tf],P0(k),[],ode_param);
  last = length(T);
  Pf = Pf | X(last,:)';
end

return