%normalnosc
x = normrnd(0,1,100,1);
figure(1)
subplot(221)
normplot(x);
subplot(222)
histfit(x,15)
[H,P,LSTAT,CV]=lillietest(x);
title(sprintf('h: %d,p: %.2f,lstat: %.2f, cv: %.2f', H,P,LSTAT,CV))
y=x(1:50)
subplot(223)
normplot(y);
subplot(224)
histfit(y,15)
[H,P,LSTAT,CV]=lillietest(y);
title(sprintf('h: %d,p: %.2f,lstat: %.2f, cv: %.2f', H,P,LSTAT,CV))
xx=x.^3;
figure(2)
subplot(221)
normplot(xx);
subplot(222)
histfit(xx,15)
[H,P,LSTAT,CV]=lillietest(xx);
title(sprintf('h: %d,p: %.2f,lstat: %.2f, cv: %.2f', H,P,LSTAT,CV))
yy=xx(1:50);
subplot(223)
normplot(yy);
subplot(224)
histfit(yy,15)
[H,P,LSTAT,CV]=lillietest(yy);
title(sprintf('h: %d,p: %.2f,lstat: %.2f, cv: %.2f', H,P,LSTAT,CV))
