xr=[0.9 0.7]';


S=[ 0.05^2 0 ; 0 0.1^2];

x1=[0.85 0.72]';
x2=[0.8  0.7]';

%S=eye(2);
dM1=(x1 -xr)'*inv(S)*(x1-xr)

dM2=(x2 -xr)'*inv(S)*(x2-xr)
