function [] = prettify_axes(ax, a, b)
%% Set style for a single plot axis
% no return 
if nargin > 1
        % Shade region between a and b
        yl = ax.YLim;  % Get current y-axis limits
        
        patch([a a b b], [yl(1) yl(2) yl(2) yl(1)], 'blue', ...
            'FaceAlpha', 0.2, 'EdgeColor', 'none');
        
        % Send shading to back so it doesn't cover the plot
        %uistack(gca().Children(1), 'bottom');
    end 
    set(ax, 'FontSize', 16, 'LineWidth',1.2, ...
        'GridAlpha', 0.4, 'GridLineStyle', '-');
    set(findall(ax, 'Type', 'Line'), 'LineWidth', 2.5, 'MarkerSize', 8);
    grid(ax, "on")
    axis(ax,"square")
    
    %title(ax, [], 'FontSize', 20);
    ax.Title.FontSize = 10;
    set(ax.Title, "FontSize", 20)

 
    