範 例 二 ： 逆矩陣與方程式組解
 
 
程式檔名稱 : linEAR2.htm
 
程式檔指令

PROC IML;
premix_p =   1.00;                       預混料所佔的百分比
total    = 500.00;                       全部預定調配成 CP18%檢定飼料總量

feed     = TOTAL - TOTAL*PREMIX_P/100;   大豆(x)與玉米(y)的總量(x+y) 
                                         大豆CP 44% 玉米CP 8%，方程式組為
                   
A = { 1    1,                                x +     y  =      feed;
      0.44 0.08};                        0.44 x + 0.08y = 0.18*total;
                                                                                                                    
b={0, 
   0};                                   預開方程式組右邊的行向量

b[1,1]=1.00*feed;                        大豆與玉米的總量 
b[2,1]=0.18*total;                       全部總飼料中含有 18% CP

z=inv(A)*b;                              以逆矩陣解出大豆與玉米的量
print A, b, z; 
run;