Input: XY(n,2) is the array of coordinates of n points x(i)=XY(i,1), y(i)=XY(i,2) Output: Par = [a b R] is the fitting circle: center (a,b) and radius R Note: this fit does not use built-in matrix functions (except "mean"), so it can be