diff --git a/code/eod_chirp_beat.py b/code/eod_chirp_beat.py
index 7f825df..072cca7 100644
--- a/code/eod_chirp_beat.py
+++ b/code/eod_chirp_beat.py
@@ -47,7 +47,7 @@ def plot_chirp(eodf, eodf1, phase, axis):
 
     y = chirp_eod * 0.4 + eod
     p, t = pd.detect_peaks(y, 0.1)
-    axis.plot(time*1000, y)
+    axis.plot(time*1000, y, color = 'royalblue')
     axis.plot(time[p]*1000, (y)[p], lw=2, color='k')
     axis.plot(time[t]*1000, (y)[t], lw=2, color='k')
     axis.spines["top"].set_visible(False)
@@ -63,18 +63,18 @@ ax2 = fig.add_subplot(222)
 ax3 = fig.add_subplot(223)
 ax4 = fig.add_subplot(224)
 
-plot_chirp(600, 650, 0, ax1)
-plot_chirp(600, 650, np.pi, ax3)
+plot_chirp(600, 650, 0, ax2)
+plot_chirp(600, 650, np.pi, ax4)
 
-plot_chirp(600, 620, 0, ax2)
-plot_chirp(600, 620, np.pi, ax4)
+plot_chirp(600, 620, 0, ax1)
+plot_chirp(600, 620, np.pi, ax3)
 
 ax1.set_ylabel('EOD [mV]', fontsize=22)
-ax1.set_title('$\Delta$f = 50 Hz', fontsize = 18)
+ax1.set_title('$\Delta$f = 20 Hz', fontsize = 18)
 ax1.yaxis.set_tick_params(labelsize=18)
 ax1.set_xticklabels([])
 
-ax2.set_title('$\Delta$f = 20 Hz', fontsize = 18)
+ax2.set_title('$\Delta$f = 50 Hz', fontsize = 18)
 ax2.set_xticklabels([])
 ax2.set_yticklabels([])
 ax3.set_ylabel('EOD [mV]', fontsize=22)