function [skew, kurtosis] = pdf_stats(x,p)

    % Define some useful bits
    dx = x(2) - x(1);

    % The expected value (From LibreTexts 4.1)
    E = @(x) sum(x .* p) * dx;

    % First moment (mean) (From LibreTexts 4.3)
    mu = E(x);
    
    % Second moment (variance) (From LibreTexts 4.3)
    sigma = sqrt(E((x - mu).^2));

    % Third moment (skewness) (From LibreTexts 4.4)
    skew = E(((x - mu)/sigma).^3);

    % Fourth moment (excess kurtosis) (From LibreTexts 4.4)
    kurtosis = E(((x - mu)/sigma).^4) - 3;

end