function    value=L(D)
%  Operator L of D (distance matrix) - not used in program
n=size(D,1);
J=eye(n)-ones(n)/n;
value=-(1/2)*J*D*J;
