习题5.2.sas

来自「SAS 源码 《应用多元分析》王学民编著的SAS程序」· SAS 代码 · 共 10 行

proc iml;
x1={4,2};
x2={3,-1};
s={6.5 1.1,1.1 8.4};
x=(x1+x2)/2;
a=inv(s)*(x1-x2);
x0={2,1};
w=t(a)*(x0-x);
print a x w;