x0 = np.deg2rad([5.0, -3.0, 0.0, 0.0])

model.plot_simulation(times, x0, v=6.0)