function out=radius_error(guess_centre,data)

x=data(:,1)-guess_centre(1);
y=data(:,2)-guess_centre(2);

r=sqrt(x.^2+y.^2);
%theta=acos(x./r);
theta=atan2(y,x);

out=sum((r-mean(r)).^2);