N=1000;

i=1;

adj = random_graph(N,10/N);

%plot(graph(adj));
v_eig=eig(adj);
[f,xi]=ksdensity(v_eig,'bandwidth',0.5);

plot(xi,f);
hold on;
std_f=2*(std(v_eig));

X=[-std_f:0.1:std_f];
Y=(std_f^2-X.^2).^(1/2).*2./(3.14*std_f^2);
plot(X,Y);
hold off;
% [gc,gcnode] = giant_component(adj);
% data(1,i) = length(gcnode);
% adj(gcnode,:) = 0;
% adj(:,gcnode) = 0;
% [gc,gcnode] = giant_component(adj);
% data(2,i) = length(gcnode);

