function plot_exp_lu()
    
    vareps = [1e-14, 1e-12, 1e-10, 1e-08, 1e-06, 1e-04, 1e-02];
    
    n_eps = size(vareps, 2);
    lineWidth = 2;
    markerSize = 15;
    
    %% matrix 1
    load("results/lu1_1.mat");
    load("results/lu1_2.mat");
    load("results/lu1_3.mat");
    
    load("results/lu1_1_bound.mat");
    load("results/lu1_2_bound.mat");
    load("results/lu1_3_bound.mat");
    
    fontSize = 19;
    %% plot 1
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 1),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 1),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 1),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 1),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=2');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu1_d2.pdf')
    hold off
    
    %% plot 2
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 2),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 2),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 2),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 2),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=5');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);

    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu1_d5.pdf')
    hold off
    
    %% plot 3
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 3),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 3),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 3),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 3),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=8');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);

    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel'); 
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu1_d8.pdf')
    hold off
    
    
    
    % ---------------------------------------------------------------------------------
    %% matrix 2
    load("results/lu2_depth2.mat");
    load("results/lu2_depth5.mat");
    load("results/lu2_depth8.mat");
    
    fontSize = 19;
    %% plot 1
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 1),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 1),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 1),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 1),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=2');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu2_d2.pdf')
    hold off
    
    %% plot 2
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 2),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 2),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 2),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 2),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=5');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu2_d5.pdf')
    hold off
    
    %% plot 3
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 3),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 3),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 3),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 3),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=8');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu2_d8.pdf')
    hold off
    
    
    % -----------------------------------------------------
    %% matrix 3
    load("results/lu3_depth2.mat");
    load("results/lu3_depth5.mat");
    load("results/lu3_depth8.mat");
    
    
    fontSize = 19;
    %% plot 1
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 1),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 1),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 1),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 1),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=2');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu3_d2.pdf')
    hold off
    
    %% plot 2
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 2),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 2),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 2),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 2),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=5');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu3_d5.pdf')
    hold off
    
    %% plot 3
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 3),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 3),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 3),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 3),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=8');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-16, 1e4]);
    yticks([ 1e-16, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu3_d8.pdf')
    hold off
    
    
    % ---------------------------------------------------------------------------------
    %% matrix 4
    load("results/lu4_depth2.mat");
    load("results/lu4_depth5.mat");
    load("results/lu4_depth8.mat");
    
    
    fontSize = 19;
    %% plot 1
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 1),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 1),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 1),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 1),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=2');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu4_d2.pdf')
    hold off
    
    %% plot 2
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 2),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 2),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 2),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 2),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=5');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu4_d5.pdf')
    hold off
    
    %% plot 3
    figure()
    set(gcf, 'Position',  [10 10 800 650]);
    
    h = gobjects(6, 1);
    
    h(1) = semilogy(1:n_eps, err_back_list1(:, 3),':b', 'Marker', 'o', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp64');
    hold on 
    h(2) = semilogy(1:n_eps, err_back_list2(:, 3),'-.r', 'Marker', '*', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'fp32');
    hold on 
    h(3) = semilogy(1:n_eps, err_back_list3(:, 3),'-.g', 'Marker', 's', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'bf16');
    hold on 
    
    h(4) = semilogy(1:n_eps, err_bound_list1(:, 3),':k', 'Marker', '.', ...
        'MarkerSize',markerSize, 'Linewidth',lineWidth, 'DisplayName', 'error bound');
    
    
    ld = legend(h([1,2,3,4]), ...
         'NumColumns',4, 'Location', 'Best', 'FontSize', fontSize, BackgroundAlpha=.3);
    legend boxoff
    rect = [0.5, 0, 0, 0.05];
    set(ld, 'Position', rect);
    
    [l, ~] = title('$\ell$=8');
    set(l,'interpreter','latex');
    
    xticklabels(split(num2str(vareps,'%.e ')));
    ylim([1e-15, 1e4]);
    yticks([ 1e-15, 1e-12, 1e-8, 1e-4, 1e-0, 1e4]);
    %extraInputs = {'interpreter','latex','fontsize',fontSize};
    %lx = xlabel('$\varepsilon$', extraInputs{:});
    
    a = get(gca,'XTickLabel');  
    set(gca,'XTickLabel',a,'fontsize',fontSize); % ,'FontWeight','bold'
    l.FontSize = fontSize+12;
    grid on;
    exportgraphics(gca, 'figures/lu4_d8.pdf')
    hold off
    
end