for i=1:100
   X=rand(50,4); %Randomly choose 50 points
   A=edm(X,X);
   V=eig(A);
   x(i)=min(min(A+eye(50)));  %This is epsilon.
   y(i)=min(abs(V));          %Smallest eigenvalue.
end
plot(x,y,'.',x,x/2);