diff --git a/code/code_cw.py b/code/code_cw.py
new file mode 100644
index 0000000..23f2af2
--- /dev/null
+++ b/code/code_cw.py
@@ -0,0 +1,33 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+freq  = 800
+freq2 = 820
+dt = 0.00001
+x = np.arange(0.0, 1.0, dt)
+eod = np.sin(x * 2 * np.pi * freq) + np.sin(x * 2 * np.pi * freq * 2) * 0.5
+eod2 = np.sin(x * 2 * np.pi * freq2) + np.sin(x * 2 * np.pi * freq2 * 2) * 0.5
+  
+fig = plt.figure(figsize=(5., 7.5))
+ax= fig.add_subplot(311)
+ax.plot(x, eod, color="darkgreen", linewidth = 1.0)
+ax.set_xlim(0.0, 0.1)
+ax.set_ylabel("voltage [mV]")
+
+
+ax= fig.add_subplot(312)
+ax.plot(x, eod2, color="crimson", linewidth = 1.0)
+ax.set_xlim(0.0, 0.1)
+ax.set_ylabel("voltage [mV]")
+
+ax= fig.add_subplot(313)
+ax.plot(x, eod + eod2 * 0.05, color="lightblue", linewidth = 1.0)
+ax.set_xlim(0.0, 0.1)
+ax.set_xlabel("time [s]")
+ax.set_ylabel("voltage [mV]")
+
+plt.tight_layout()
+plt.savefig("eods.pdf")
+plt.show()
+
+
diff --git a/code/eod_cl.py b/code/eod_cl.py
new file mode 100644
index 0000000..e69de29
diff --git a/code/eod_sim_rs.py b/code/eod_sim_rs.py
new file mode 100644
index 0000000..e69de29