function [x,y] = randomdiskpoint(n)

diskradius=10;

rtheta=rand(n,2);
rtheta(:,1)=rtheta(:,1).*diskradius;
rtheta(:,2)=rtheta(:,2).*(2*pi);
[xx,yy]=pol2cart(rtheta(:,2),rtheta(:,1));
xx = xx+10;
yy = yy+10;

diskradius1=10;

rtheta1=rand(n,2);
rtheta1(:,1)=rtheta1(:,1).*diskradius1;
rtheta1(:,2)=rtheta1(:,2).*(2*pi);
[xx1,yy1]=pol2cart(rtheta1(:,2),rtheta1(:,1));
xx1 = xx1-10;
yy1 = yy1-10;

x = [xx;xx1];
y = [yy;yy1];
%scatter(xx,yy);