clear z y
subplot(111);
for n=1:101
	j=0;
	z(n)=(n-51)/50;
	k=z(n)*pi;
	for i=1:1000,
 	if (er(i)<k)
 	 j = j+1;
	 end
	end
	y(n) = j/1000;
end
plot(z,y)
grid
title('distribution function, zeta = .707, Int&Ld')
xlabel('x in pi radians')
ylabel('probability of phase < x')

