runSim.m

function results = runSim(model, sim_in)
% Runs a Simulink model. The `model` takes as input
% a workspace variable `sim_in` which is a matrix where
% the first column is a time stamp and the remaining
% columns are individual inputs. The model outputs to
% a workspace variable `sim_out` which is a timestamp
% object.
% The function parses the `sim_out` and returns `result`
% in the same form as `sim_in`.
stop_t = sim_in(end,1);
steps = diff(sim_in(:, 1));
minStep = min(steps(:));
set_param(model, 'StopTime', num2str(stop_t), 'MaxStep', num2str(minStep));
% set_param(model, 'StopTime', num2str(stop_t));
sim(model);
% results = [sim_out.Time sim_out.Data];
[y, t] = resample(sim_out.Data, sim_out.Time, 1 / minStep);
results = [t, y];
end