diff --git a/roi_pac.m b/roi_pac.m
index 7a333c5..9f65e02 100644
--- a/roi_pac.m
+++ b/roi_pac.m
@@ -1,4 +1,4 @@
-% roi_pac() - Compute phase-amplitude coupling (PAC) between ROIs  (cf. Zandvoort and Nolte, 2021). 
+% roi_pac() - Computes phase-amplitude coupling (PAC) between ROIs  (cf. Zandvoort and Nolte, 2021). 
 % Wrapper function for pac_bispec().
 %
 % Usage:
diff --git a/simulations/sim_tde_bispec.m b/simulations/sim_tde_bispec.m
index feece70..06e4a43 100644
--- a/simulations/sim_tde_bispec.m
+++ b/simulations/sim_tde_bispec.m
@@ -68,9 +68,10 @@
 [peak_val, peak_idx] = max(aT); % extract estimated delay/peak
 est_delay = delay_scale(peak_idx);
 
-figure; plot(delay_scale, aT)
+figure; plot(delay_scale, aT, 'black')
 xline(est_delay, '--r')
 xlabel('Time (s)')
+ylabel('a.u.')
 title(sprintf('TDE | Method %d', method))
 subtitle("\tau = " + num2str(est_delay) + " (s)")
 grid on