更郁闷的是，我连求反函数的INV（）函数都不能用，sas同样提示说没这个函数。是不是我什么东西没打开啊？

在IML里面你不能直接用data step,如果你要把矩阵存为data set,你可以

proc IML;
I130=I(130);
create a from I130;
append from I130;
quit;
proc print data=a;run;

谢谢您的回复！我现在有遇到问题了，请您帮助解答。
我想计算leontif拟矩阵，其计算公式为L=inv(E-A);E是单位阵，就是下面的I130，A是中间消耗矩阵，其值在work.test中
下面是用的命令：

proc iml;
I130=I(130);
create a from I130;
append from I130;
quit;
proc print data=E;
run;

proc iml;
use work.test;   
read all into A;
quit;

proc iml;
b=E-A;


但是到最后一步出错了，请问是怎么回事，下面是出错的信息：
474  proc iml;
NOTE: IML Ready
475  b=E-A;
ERROR: (execution) Matrix has not been set to a value.

operation : - at line 475 column 4
operands  : E, A

E      0 row       0 col     (type ?, size 0)


A      0 row       0 col     (type ?, size 0)


statement : ASSIGN at line 475 column 1

你前面退出iml了，E和A就没了
在用iml，要重新建立E和A

不好意思，我把quit都去掉了还是不行。
proc iml;
I130=I(130);
create E from I130;
append from I130;
proc print data=E;
run;

proc iml;
use peter.A;   
read all into A;

proc iml;
b=E-A;


请问该怎么处理啊？
41   proc iml;
NOTE: IML Ready
42   I130=I(130);
43   create E from I130;
44   append from I130;
NOTE: Exiting IML.
NOTE: The data set WORK.E has 130 observations and 130 variables.
NOTE: PROCEDURE IML used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


45   proc print data=E;
46   run;

NOTE: There were 130 observations read from the data set WORK.E.
NOTE: PROCEDURE PRINT used (Total process time):
      real time           0.01 seconds
      cpu time            0.00 seconds


47   proc iml;
NOTE: IML Ready
48   use peter.A;
49   read all into A;
NOTE: Exiting IML.
NOTE: PROCEDURE IML used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds


50   proc iml;
NOTE: IML Ready
51   b=E-A;
ERROR: (execution) Matrix has not been set to a value.

operation : - at line 51 column 4
operands  : E, A

E      0 row       0 col     (type ?, size 0)


A      0 row       0 col     (type ?, size 0)


statement : ASSIGN at line 51 column 1

不理解你为什么要用那么多的proc iml
你把所有要用iml算的放在一个iml不久行了
而且你的E，是个dataset，不是矩阵
你的A，是矩阵，但是用了一个新的iml，之后内存里面就没有A了
所以在你的最后一个iml里面，是没有E和A的
你要重新定义E和A

用fcmp试试，不过不知道里面的矩阵函数是否合适你用。
data sample1;
input x1-x2;
cards;
1 2
3 4
;
run;
data sample2;
input x1-x2;
cards;
1 2
1 2
;
run;
proc fcmp;
        array mat1[1]/nosymbols;
        array mat2[1]/nosymbols;
        array result[2,2]/nosymbols;
        rc=read_array('work.sample1',mat1);
        rc=read_array('work.sample2',mat2);
        call mult(mat1,mat2,result);
        rc=write_array('work.out',result);
quit;

13# 275769263 谢谢您的回复！

正解