前面是data步，
  proc iml;
   n=20; p=3;
   m0={4 50 10};
   use d321;
   xa={x1 x2 x3};
   read all  var xa into x; print x;
   ln={[20] 1} ; print ln;
   x0=(ln*x)/n; print x0;
   xm=x0-m0; print xm;
   mm=i(20)-j(20,20,1)/n;print mm;
   a=x`*mm*x; print a;
a应该是样本协方差阵，那倒数第二步的i（20）-j（20,20,1）是什么意思啊？