clear;N=1000;
x=rand(N,1)+i*rand(N,1);
length(find(abs(x)<1))/N*4
xin=x(find(abs(x)<1));
xout=x(find(abs(x)>=1));
plot(real(xin),imag(xin),'b.',real(xout),imag(xout),'r.');
hold('on');
plot(cos(deg2rad(0:90)),sin(deg2rad(0:90)),'k:');
axis('square');
axis([0,1,0,1]);

