diff --git a/docs/tutorials/Lindblad_dynamics_simulation.rst b/docs/tutorials/Lindblad_dynamics_simulation.rst
index 9602d2299..898cd6a7a 100644
--- a/docs/tutorials/Lindblad_dynamics_simulation.rst
+++ b/docs/tutorials/Lindblad_dynamics_simulation.rst
@@ -244,8 +244,8 @@ parameters can be systematically analyzed - we encourage you to try it!
     ax.set_xlabel('$t$', fontsize = fontsize)
     ax.set_title('Mean Bloch vector vs. $t$', fontsize = fontsize)
 
-    display(plot_bloch_vector([x_mean[-1], y_mean[-1], z_mean[-1]],
-                      f'Mean Bloch vector at $t = {t_eval[-1]}$'))
+    plot_bloch_vector([x_mean[-1], y_mean[-1], z_mean[-1]],
+                      f'Mean Bloch vector at $t = {t_eval[-1]}$')
 
     if N > 1 and ((abs(x_mean[-1]) > 1e-5 and abs(x_data[0, -1] / x_mean[-1] - 1) > 1e-5 or
                   (abs(z_mean[-1]) > 1e-5 and abs(z_data[1, -1] / z_mean[-1] - 1) > 1e-5))):
diff --git a/docs/tutorials/Rabi_oscillations.rst b/docs/tutorials/Rabi_oscillations.rst
index e5a84610c..f498a56ad 100644
--- a/docs/tutorials/Rabi_oscillations.rst
+++ b/docs/tutorials/Rabi_oscillations.rst
@@ -130,8 +130,8 @@ gate.
         ax.set_title('Bloch vector vs. $t$', fontsize = fontsize)
         plt.show()
 
-        display(plot_bloch_vector([x_data[-1], y_data[-1], z_data[-1]],
-                                  f'Bloch vector at $t = {t_eval[-1]}$'))
+        plot_bloch_vector([x_data[-1], y_data[-1], z_data[-1]],
+                                  f'Bloch vector at $t = {t_eval[-1]}$')
 
     plot_qubit_dynamics(sol, t_eval, X, Y, Z)