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)