The current method for computing the default simulation time for discrete-time simulations uses a call to scipy.signal._default_response_times(), but this uses an eigenvalue calculation that assumes continuous time. For discrete time systems, need to figure out what the right default values are for the length of a simulation for impulse, step, and initial condition responses. See also discussion in issue #239.