function [mu, sigma] = estimation_mu_Sigma(X)

    mu = mean(X);
    
    X_c = X - mu;
    
    sigma = X_c.' * X_c / length(X);
    
    mu = mu.';

end