%% Input: adjacency matrix A

%D=diag(1./sum(A));

%An=sqrt(D)*A*sqrt(D);
%[v,d]=eigs(An,4,'largestreal');
%[v,d]=eigs(A,4,'largestreal');

L=lap(A);
%Ln=sqrt(D)*lap(A)*sqrt(D);
[v,d]=eigs(L,4,'smallestreal');

gplot(A,v(:,2:3));
