close all;
clear; clc;

N = [1750, 500];
L = [128, 2048, 4096];

for n = N
    for l = L
        % Read the data
        data = load(["N",num2str(n),"_L",num2str(l),".txt"]);
        [skew, kurt] = pdf_stats(data(:,1),data(:,2));
        fprintf("N=%i, L=%i: skew=%f, kurt=%f\n",n,l,skew,kurt)
    end
end