function [] = tcl(n, p, nc)
  xbasc();
  X=rand(n,p);
  Z=sqrt(12/p)*(sum(X,'c') - p/2); // somme des colonnes de X, centrage et
                                 // renormalisation 
  histplot(nc,Z)
  C=[-5:1/1000:+5]; 
  plot2d(C,exp(-C.^2/2)/sqrt(2*%pi),3) // densité de la loi N(0,1) 
endfunction