diff --git a/A6-simscape-nass/nass_2_hac.m b/A6-simscape-nass/nass_2_hac.m
index 0481a94..81082bd 100644
--- a/A6-simscape-nass/nass_2_hac.m
+++ b/A6-simscape-nass/nass_2_hac.m
@@ -507,11 +507,11 @@ tiledlayout(3, 1, 'TileSpacing', 'Compact', 'Padding', 'None');
 ax1 = nexttile([2,1]);
 hold on;
 plot(f, abs(squeeze(freqresp(G_m1_vc(1,1), f, 'Hz'))), 'color', colors(1,:), ...
-     'DisplayName', '$f_{ni}/f_i$ - $\Omega_z = 0$')
+     'DisplayName', '$\epsilon_{\mathcal{L}i}/f_i$ - $\Omega_z = 0$')
 plot(f, abs(squeeze(freqresp(G_m1_vc_Rz_slow(1,1), f, 'Hz'))), 'color', colors(2,:), ...
-     'DisplayName', '$f_{ni}/f_i$ - $\Omega_z = 36$ deg/s')
+     'DisplayName', '$\epsilon_{\mathcal{L}i}/f_i$ - $\Omega_z = 36$ deg/s')
 plot(f, abs(squeeze(freqresp(G_m1_vc_Rz_fast(1,1), f, 'Hz'))), 'color', colors(3,:), ...
-     'DisplayName', '$f_{ni}/f_i$ - $\Omega_z = 360$ deg/s')
+     'DisplayName', '$\epsilon_{\mathcal{L}i}/f_i$ - $\Omega_z = 360$ deg/s')
 for i = 1:5
     for j = i+1:6
         plot(f, abs(squeeze(freqresp(G_m1_vc(i,j), f, 'Hz'))), 'color', [colors(1,:), 0.2], ...