The time response routines in matlab.py return values in a different order from their counterparts in timeresp.py. For instance, one does

t, y = timeresp.step_response(sys)

but

y, t = matlab.step(sys)

Worse, for matlab.impulse, matlab.initial, and matlab.lsim, the actual behavior (y,t) is inconsistent with what is in the docstring (t,y). This is confusing, and I think the order should be consistent. I guess it would make the most sense to make all of these consistent with MATLAB, which is (y,t).

This is a nasty change, though, that is likely to break existing code people may have. Any thoughts?