function [A,B,C] = LMSRecta(x,y)
%Ajusta por minimos cuadrados a recta Ax + By + C = 0 calcula eje inercia
% Devuelve [A,B,C]. n(A,B) unitario. Viene bien para calculo dist(P,r)
% x,y contiene las x (columnas)y las y (filas) de los pix de cont

sum_x=0;               %sumatorio x
sum_y=0;               %sumatorio j
Ixx = 0;
Iyy = 0;
Ixy = 0;
n=0;

% Media de las coord X e y (centro gravedad)
cx=mean(x);
cy=mean(y);
centro=[cx,cy];

n=length(x);
for k=1:n
    Ixx = Ixx + (x(k) - cx)^2;
    Iyy = Iyy + (y(k) - cy)^2;
    Ixy = Ixy + (x(k) - cx)*(y(k) - cy);
end

%Tita es el angulo de la recta con el eje x (horiz)
tita=0.5*atan2( 2*Ixy, Ixx-Iyy );
m=tan(tita);
b=cy-m*cx;
% mx-y+b=0. 
A=m; B=-1; C=b;

%Normalizo n(A,B)=1
n=[A,B];
normAB=norm(n);
A = A/normAB;
B = B/normAB;
C = C/normAB;

end


