习题5.2.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;