function randomg

n = 20;     % number of nodes
p = .05;     % connection proba

x = cos(2*pi*(0:n-1)/n);
y = sin(2*pi*(0:n-1)/n);       % coordinates for (primitive) drawing

plot(x,y,'ro','MarkerSize',10)
axis([-1.2 1.2 -1.2 1.2])
axis('square')
hold on

for j=1:n
    for k=1:n
        q = rand(1);
        if (q <= p & k ~= j)
            plot([x(j) x(k)],[y(j) y(k)],'b-','LineWidth',3)
            drawnow
        end
    end
end

hold off