diff --git a/work/week12/starting_point/test_computeFrequency.txt b/work/week12/starting_point/test_computeFrequency.txt
deleted file mode 100644
index 9c94750..0000000
--- a/work/week12/starting_point/test_computeFrequency.txt
+++ /dev/null
@@ -1,10 +0,0 @@
-0.00
-0.10
-0.20
-0.30
-0.40
--0.50
--0.40
--0.30
--0.20
--0.10
diff --git a/work/week12/starting_point/test_computeFrequency_1D.txt b/work/week12/starting_point/test_computeFrequency_1D.txt
new file mode 100644
index 0000000..bf583f5
--- /dev/null
+++ b/work/week12/starting_point/test_computeFrequency_1D.txt
@@ -0,0 +1,10 @@
+0.00
+0.63
+1.26
+1.88
+2.51
+-3.14
+-2.51
+-1.88
+-1.26
+-0.63
diff --git a/work/week12/starting_point/test_computeFrequency_2D.txt b/work/week12/starting_point/test_computeFrequency_2D.txt
new file mode 100644
index 0000000..78e6d5d
--- /dev/null
+++ b/work/week12/starting_point/test_computeFrequency_2D.txt
@@ -0,0 +1,10 @@
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
+0.00 0.63 1.26 1.88 2.51 -3.14 -2.51 -1.88 -1.26 -0.63
diff --git a/work/week12/starting_point/test_fftfreq.py b/work/week12/starting_point/test_fftfreq.py
index 3017797..e9c2cf7 100644
--- a/work/week12/starting_point/test_fftfreq.py
+++ b/work/week12/starting_point/test_fftfreq.py
@@ -1,13 +1,18 @@
 import numpy as np
 
 n = 10
-timestep = 1
+timestep = 10
 
-freq = np.fft.fftfreq(n, d=timestep)
-print(freq)
+#  Generate 1D frequency coordinates
+freq_1D = np.fft.fftfreq(n, d=timestep) * 2 * np.pi * n
+print(freq_1D)
+#  Generate 2D frequency coordinates
+freq_2D = np.einsum('i,j->ij', np.ones(n), freq_1D)
+print(freq_2D)
 
-file_data = np.hstack((freq))
-np.savetxt("test_computeFrequency.txt", file_data, delimiter=" ", fmt='%.2f')
-
-# freqArray2D = np.meshgrid(freq, freq)
-# print(freqArray2D)
+#  Store 1D frequency coordinates
+np.savetxt("test_computeFrequency_1D.txt", freq_1D, delimiter=" ",
+           fmt='%.2f')
+#  Store 2D frequency coordinates
+np.savetxt("test_computeFrequency_2D.txt", freq_2D, delimiter=" ",
+           fmt='%.2f')