为什么不能用proc corr呢  proc iml你可以google一下

如果你只要作这一个表的
写一个iml的算法真是得不偿失
不能用proc corr的话，用Excel也能算出来
或者用matlab这样的矩阵语言也还简单

data test;
   infile cards;
   input company $ 16  x1 x2 x3;
   cards;
General Motors               126974        4224            173297
Ford                          96933        3835            160893
Exxon                         86656        3510            83219
IBM                           63438        3758             77734
General Electric              55264        3939            128344
Mobil                         50976        1809            39080
Philip Morris                 39069        2946            38528
Chrysler                      36156        359             51038
Du Pont                       35209        2480            34715
Texaco                        32416        2413             25636
;
run;
ods trace output;
PROC IML;
   use test;
   read all var _num_ into x;
   nm={x1 x2 x3};
/* Module to compute correlations */
start corr;
   n = nrow(x); /* number of observations */
   sum = x[+,] ; /* compute column sums */
   xpx = t(x)*x-t(sum)*sum/n; /* compute sscp matrix */
   s = diag(1/sqrt(vecdiag(xpx))); /* scaling matrix */
   corr = s*xpx*s; /* correlation matrix */
   print "Correlation Matrix",,corr[rowname=nm colname=nm] ;
finish corr;
/* Sample run */
run corr;
quit;

真的非常感谢大家的热心帮助。没想到这题用proc iml这么麻烦，还是proc corr方便啊，不过要求如此，咱屁民只能乖乖听话。