S = 4;
n = 30;
mu=4.4311
sig2=5.365
sig = sqrt(sig2)
beta=sig2/mu;alpha = mu/beta;
x = 0:0.1:10;
plot(x,gampdf(x,alpha,beta))
X = gamrnd(alpha,beta,n,S);
format SHORT
tab = [X;mean(X);median(X);std(X)];
%printmat(tab,'gam','Obs1 Obs2 Obs3 Obs4 Obs5 Obs6 Obs7 Obs8 Obs9 Obs10 Mean Median Std a a','S1 S2 S3 S4')
printmat(tab)