>> X = rand(2,10)

X =

    0.8147    0.1270    0.6324    0.2785    0.9575    0.1576    0.9572    0.8003    0.4218    0.7922
    0.9058    0.9134    0.0975    0.5469    0.9649    0.9706    0.4854    0.1419    0.9157    0.9595

>> i=1;
>> d = sum((X - repmat(X(:,1),1,10)).^2,1)

d =

         0    0.4730    0.6865    0.4164    0.0239    0.4360    0.1970    0.5838    0.1545    0.0034

>> f = rand(1,10)

f =

    0.6557    0.0357    0.8491    0.9340    0.6787    0.7577    0.7431    0.3922    0.6555    0.1712

>> mean(d)

ans =

    0.2975

>> sh = 1 - d./mean(d)

sh =

    1.0000   -0.5903   -1.3080   -0.3997    0.9197   -0.4658    0.3376   -0.9625    0.4805    0.9886

>> f(1) = f(1)/(1+0.92+0.34+0.48+0.99)

f =

    0.1758    0.0357    0.8491    0.9340    0.6787    0.7577    0.7431    0.3922    0.6555    0.1712