r = 40; n=500;
S2 = randn(3,n)';
for i = 1:n;
S2(i,:) = r*(S2(i,:)./norm(S2(i,:)));
end
X = S2;
GD=zeros(n);
for i =1:n
for j=1:n
GD(i,j)=(2*r*asin( 0.5*( norm( X(i,:)-X(j,:) ) )/r ))^2*(1+0.2*randn);