X=[216,203,197,208,206,209,206,208,202,203,206,213,218,207,208,...
    202,194,203,213,211,193,213,208,208,204,206,204,206,208,209,...
    213,203,206,207,196,201,208,207,213,208,210,208,211,211,214,...
    220,211,203,216,224,211,209,218,214,219,211,208,221,211,218,...
    218,190,219,211,208,199,214,207,207,214,206,217,214,201,212,...
    213,211,212,216,206,210,216,204,221,208,209,214,214,199,204,...
    211,201,216,211,209,208,209,202,211,207,202,205,206,216,206,...
    213,206,207,200,198,200,202,203,208,216,206,222,213,209,219];
[H,p]=jbtest(X,0.05)

[mu1,sig1,mu_ci,sig_ci]=normfit(X,0.05); mu=[mu1,mu_ci']

sig=[sig1, sig_ci']