function [X,Y]= circ(a,b,r,rn);
z=linspace(0,pi,100);
X=r*sin(z)+a+(rand(1,100)-0.5)*rn;
Y=r*cos(z)+b+(rand(1,100)-0.5)*rn;
end;
