axes[1].plot(t, n[:, i_trajectory], color=(gray, gray, gray))