PROC IML; 

S = {2 1 1, 0 4 2, 0 0 5};  S = S*S`; 
S = 10*S/8; 
IN = inv(S); 
m1 = {2,-1,1}; m2 = {-2,0,1}; m3 = {1,-1,1}; 
print S in m1 m2 m3; 
D1 = 0.5*m1`*in*m1||( - m1`*IN);
D2 = 0.5*m2`*in*m2||( - m2`*IN);
D3 = 0.5*m3`*in*m3||( - m3`*IN);
D = D1//D2//D3; 

Y = {1,2,3};
distance = D*({1}//Y); 
print D Y distance;