diff --git a/lectures/python_essentials.md b/lectures/python_essentials.md
index 8589b06e..ab8a4a84 100644
--- a/lectures/python_essentials.md
+++ b/lectures/python_essentials.md
@@ -1068,10 +1068,11 @@ code.
 ```{code-cell} python3
 import numpy as np
 
+rng = np.random.default_rng()
 n = 100
 ϵ_values = []
 for i in range(n):
-    e = np.random.randn()
+    e = rng.standard_normal()
     ϵ_values.append(e)
 ```
 
@@ -1085,8 +1086,9 @@ for i in range(n):
 Here's one solution.
 
 ```{code-cell} python3
+rng = np.random.default_rng()
 n = 100
-ϵ_values = [np.random.randn() for i in range(n)]
+ϵ_values = [rng.standard_normal() for i in range(n)]
 ```
 
 ```{solution-end}