Something strange about the return shape of the frequency response. For a square 2x2 system:

sys(np.array(0.0), squeeze=False).shape
(2, 2, 1)

sys(0.0, squeeze=False).shape
(2, 2)

Originally posted with the wrong result mapping by @forgi86 in #593 (comment)