function simult(X,alpha)
%function simult(X,alpha)
%
% Costruisce intervalli simultanei per la media
% del campione
% 
%
% X : matrice del campione, n x p
% alpha : liv. significativita (es. alpha=0.05)
%
% Per p=2 mostra anche il grafico

[n,p]=size(X);
m=mean(X);
S=cov(X);

for k=1:p,
   a=zeros(p,1); a(k)=1;
   l = sqrt(p*(n-1)/(n-p)*finv(1-alpha,p,n-p)*(a'*S*a)/n);
   A(k)=m(k)-l; B(k)=m(k)+l;
fprintf('var. %d: intervallo simult.  %d %d\n',k,A(k),B(k));
end

if p==2,
intervs1(A(1),B(1),A(2),B(2));
end