bode_plot phase wrapping incorrect for multiple systems.
Minimum example:
# two RHP's G = zpk([],[1,1], gain=1) bode_plot(G) bode_plot((G,G))
The second call shows a figure with different phases for each (identical) system.