In the current simulation functions (initial_response, step_response, forced_response) the time values must be equally spaced. For discrete time systems, SciPy’s dlsim function allows arbitrary time instants to be passed as inputs and interpolates between them. Perhaps we should do the same (with equally spaced time points as a special case that uses a more efficient computation?