1function varargout = interpfast(t_in, dt, t_all, varargin)
2
3
4
5
6
7t = max(0, t_in);
8ifloor = floor(t/dt) + 1;
9tmod = mod(t, dt)/dt;
10for i = 1:nargin-3
11 if t == t_all(end)
12 varargout{i} = varargin{i}(:, end);
13 else
14 varargout{i} = (1-tmod) * varargin{i}(:, ifloor) + tmod * varargin{i}(:, ifloor + 1);
15 end
16end
17
18end
19