diff --git a/doc/pub/week39/html/week39-reveal.html b/doc/pub/week39/html/week39-reveal.html
index a7c3d44e7..8e10ec464 100644
--- a/doc/pub/week39/html/week39-reveal.html
+++ b/doc/pub/week39/html/week39-reveal.html
@@ -206,11 +206,11 @@ which results in optimal scores (for example the smallest MSE or an \( R2=1 \)).
     <span style="color: #8B008B; font-weight: bold">return</span> np.sum((y_data-y_model)**<span style="color: #B452CD">2</span>)/n
 <span style="color: #228B22"># A seed just to ensure that the random numbers are the same for every run.</span>
 <span style="color: #228B22"># Useful for eventual debugging.</span>
-np.random.seed(<span style="color: #B452CD">315</span>)
+np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
@@ -276,7 +276,7 @@ np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
@@ -340,7 +340,7 @@ np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
diff --git a/doc/pub/week39/html/week39-solarized.html b/doc/pub/week39/html/week39-solarized.html
index a7faaa359..f8a11d320 100644
--- a/doc/pub/week39/html/week39-solarized.html
+++ b/doc/pub/week39/html/week39-solarized.html
@@ -290,11 +290,11 @@ which results in optimal scores (for example the smallest MSE or an \( R2=1 \)).
     <span style="color: #8B008B; font-weight: bold">return</span> np.sum((y_data-y_model)**<span style="color: #B452CD">2</span>)/n
 <span style="color: #228B22"># A seed just to ensure that the random numbers are the same for every run.</span>
 <span style="color: #228B22"># Useful for eventual debugging.</span>
-np.random.seed(<span style="color: #B452CD">315</span>)
+np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
@@ -360,7 +360,7 @@ np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
@@ -424,7 +424,7 @@ np.random.seed(<span style="color: #B452CD">2021</span>)
 
 n = <span style="color: #B452CD">100</span>
 x = np.random.rand(n)
-y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)
+y = np.exp(-x**<span style="color: #B452CD">2</span>) + <span style="color: #B452CD">1.5</span> * np.exp(-(x-<span style="color: #B452CD">2</span>)**<span style="color: #B452CD">2</span>)+ np.random.randn(n)
 
 Maxpolydegree = <span style="color: #B452CD">5</span>
 X = np.zeros((n,Maxpolydegree-<span style="color: #B452CD">1</span>))
diff --git a/doc/pub/week39/html/week39.html b/doc/pub/week39/html/week39.html
index 6c0ca6be3..42321f373 100644
--- a/doc/pub/week39/html/week39.html
+++ b/doc/pub/week39/html/week39.html
@@ -295,11 +295,11 @@ which results in optimal scores (for example the smallest MSE or an \( R2=1 \)).
     <span style="color: #008000; font-weight: bold">return</span> np<span style="color: #666666">.</span>sum((y_data<span style="color: #666666">-</span>y_model)<span style="color: #666666">**2</span>)<span style="color: #666666">/</span>n
 <span style="color: #408080; font-style: italic"># A seed just to ensure that the random numbers are the same for every run.</span>
 <span style="color: #408080; font-style: italic"># Useful for eventual debugging.</span>
-np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>seed(<span style="color: #666666">315</span>)
+np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>seed(<span style="color: #666666">2021</span>)
 
 n <span style="color: #666666">=</span> <span style="color: #666666">100</span>
 x <span style="color: #666666">=</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>rand(n)
-y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)
+y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)<span style="color: #666666">+</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>randn(n)
 
 Maxpolydegree <span style="color: #666666">=</span> <span style="color: #666666">5</span>
 X <span style="color: #666666">=</span> np<span style="color: #666666">.</span>zeros((n,Maxpolydegree<span style="color: #666666">-1</span>))
@@ -365,7 +365,7 @@ np<span style="color: #666666">.</span>random<span style="color: #666666">.</spa
 
 n <span style="color: #666666">=</span> <span style="color: #666666">100</span>
 x <span style="color: #666666">=</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>rand(n)
-y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)
+y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)<span style="color: #666666">+</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>randn(n)
 
 Maxpolydegree <span style="color: #666666">=</span> <span style="color: #666666">5</span>
 X <span style="color: #666666">=</span> np<span style="color: #666666">.</span>zeros((n,Maxpolydegree<span style="color: #666666">-1</span>))
@@ -429,7 +429,7 @@ np<span style="color: #666666">.</span>random<span style="color: #666666">.</spa
 
 n <span style="color: #666666">=</span> <span style="color: #666666">100</span>
 x <span style="color: #666666">=</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>rand(n)
-y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)
+y <span style="color: #666666">=</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>x<span style="color: #666666">**2</span>) <span style="color: #666666">+</span> <span style="color: #666666">1.5</span> <span style="color: #666666">*</span> np<span style="color: #666666">.</span>exp(<span style="color: #666666">-</span>(x<span style="color: #666666">-2</span>)<span style="color: #666666">**2</span>)<span style="color: #666666">+</span> np<span style="color: #666666">.</span>random<span style="color: #666666">.</span>randn(n)
 
 Maxpolydegree <span style="color: #666666">=</span> <span style="color: #666666">5</span>
 X <span style="color: #666666">=</span> np<span style="color: #666666">.</span>zeros((n,Maxpolydegree<span style="color: #666666">-1</span>))
diff --git a/doc/pub/week39/ipynb/ipynb-week39-src.tar.gz b/doc/pub/week39/ipynb/ipynb-week39-src.tar.gz
index 54492e2af..ac5189fac 100644
Binary files a/doc/pub/week39/ipynb/ipynb-week39-src.tar.gz and b/doc/pub/week39/ipynb/ipynb-week39-src.tar.gz differ
diff --git a/doc/pub/week39/ipynb/week39.ipynb b/doc/pub/week39/ipynb/week39.ipynb
index ffd705dee..65c5ab2bf 100644
--- a/doc/pub/week39/ipynb/week39.ipynb
+++ b/doc/pub/week39/ipynb/week39.ipynb
@@ -65,11 +65,11 @@
     "    return np.sum((y_data-y_model)**2)/n\n",
     "# A seed just to ensure that the random numbers are the same for every run.\n",
     "# Useful for eventual debugging.\n",
-    "np.random.seed(315)\n",
+    "np.random.seed(2021)\n",
     "\n",
     "n = 100\n",
     "x = np.random.rand(n)\n",
-    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)\n",
+    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)\n",
     "\n",
     "Maxpolydegree = 5\n",
     "X = np.zeros((n,Maxpolydegree-1))\n",
@@ -144,7 +144,7 @@
     "\n",
     "n = 100\n",
     "x = np.random.rand(n)\n",
-    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)\n",
+    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)\n",
     "\n",
     "Maxpolydegree = 5\n",
     "X = np.zeros((n,Maxpolydegree-1))\n",
@@ -215,7 +215,7 @@
     "\n",
     "n = 100\n",
     "x = np.random.rand(n)\n",
-    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)\n",
+    "y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)\n",
     "\n",
     "Maxpolydegree = 5\n",
     "X = np.zeros((n,Maxpolydegree-1))\n",
diff --git a/doc/src/week39/week39.do.txt b/doc/src/week39/week39.do.txt
index 816605db6..e01b937df 100644
--- a/doc/src/week39/week39.do.txt
+++ b/doc/src/week39/week39.do.txt
@@ -42,11 +42,11 @@ def MSE(y_data,y_model):
     return np.sum((y_data-y_model)**2)/n
 # A seed just to ensure that the random numbers are the same for every run.
 # Useful for eventual debugging.
-np.random.seed(315)
+np.random.seed(2021)
 
 n = 100
 x = np.random.rand(n)
-y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)
+y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)
 
 Maxpolydegree = 5
 X = np.zeros((n,Maxpolydegree-1))
@@ -110,7 +110,7 @@ np.random.seed(2021)
 
 n = 100
 x = np.random.rand(n)
-y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)
+y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)
 
 Maxpolydegree = 5
 X = np.zeros((n,Maxpolydegree-1))
@@ -172,7 +172,7 @@ np.random.seed(2021)
 
 n = 100
 x = np.random.rand(n)
-y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)
+y = np.exp(-x**2) + 1.5 * np.exp(-(x-2)**2)+ np.random.randn(n)
 
 Maxpolydegree = 5
 X = np.zeros((n,Maxpolydegree-1))