import matplotlib.pyplot as plt
from dtk.bicycle import (benchmark_matrices,
                         benchmark_state_space_vs_speed)
from dtk.control import eig_of_series, sort_modes, plot_phasor

M, C1, K0, K2 = benchmark_matrices()
v, A, B = benchmark_state_space_vs_speed(M, C1, K0, K2)
evals, evecs = sort_modes(*eig_of_series(A))
plot_phasor(evals[25], evecs[25])