close all
hold off
axis off
y=.5^(1/2);
p=input('input occupation probability')
hold off
for j=0:10
for k=0:7
	hold on
plot(j+.5,2*k*y,'o')
end
end
for j=0:10
for k=1:7
plot(j,(2*k-1)*y,'o')
end
end
pause(2)
for j=0:10
for k=0:7
  x=rand;
if x<p
	hold on
plot(j+.5,2*k*y,'*')
end
end
end
for j=0:10
for k=1:7
	x=rand;
	if x<p
plot(j,(2*k-1)*y,'*')
end
end
end
