%
%  funcion randT2
%
%  Esta funcion genera una muestra de tamano N de una ley T^2 de
%  Hotelling con p y n grados de libertad.
%
   function t=randT2(p,n,N)
%
   n=n+1;
   for i=1:N
      X=randn(n,p);
      m=mean(X);
      S=cov(X,1);
      t(i,1)=(n-1)*m*inv(S)*m';
   end
%  numero de intervalos en el hsitograma
   if N<10^4
      k=sqrt(N);
   else
      k=1+3.22*log10(N);
   end
   int=0:max(t)/k:max(t);
   hist(t,int)
   h=findobj(gca,'Type','patch');
   set(h,'FaceColor','c','EdgeColor','b')