May I ask how I can find 95% confidence interval of skew_FP() from these codes?

rng(1000)

nn=143;

for n=1:1000;

r=rand(nn,1);

ss(n)=std(r);

m(n)=mean(r);

med(n)=median(r);

nu(n)=numel(r)

skew_FP(n) = sum((r-m(n)).^3)/(nn-1)/ss(n).^3

end

hist(skew_FP())

xlabel('Skew_FP()')

ylabel('Frequency')

Answer 1

[Jeff Miller]

LBUB = prctile(skew_FP(),[2.5,97.5]);