r=1;
r1=7;
theta=2*pi*rand();
z1=rand();
    	R1=r+(r1-r)*(rand())^.5;  %creates random uniform distribution
x1=R1*cos(theta);
y1=R1*sin(theta);
plot (x1,y1,'x')