run = [68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81];
gap = [9.0, 9.2, 9.4, 9.6, 9.8, 10.0, 7.72, 12.12, 13.0, 15.0, 20.0, 25.0, 8.5, 8.0];
k = [3.053077, 2.978438, 2.906377, 2.836604, 2.769246, 2.703887, 2.795831, 2.119360];
phase_match_err_enter = [-36.534376, 8.816163, 54.951577, 102.197338, 150.453645, 199.524725, 130.835626, 773.194763];

%Plot the results
figure;
hold on;
plot(k, phase_match_err_enter, '*b');
title('Phase Match Err vs K');
xlabel('K');
ylabel('Phase Match Err (deg)');
grid on;
zoom on;
hold off;

%Plot the results
figure;
hold on;
plot(k, phase_match_err_enter .* (1 + 0.5 * k.^2), '*b');
title('Phase Match Err * (1 + K**2 / 2) vs K');
xlabel('K');
ylabel('Phase Match Err  * (1 + K**2 / 2) (deg)');
grid on;
zoom on;
hold off;

