%% Retta per 2 punti nella forma implicita ax+by+c=0

function [a , b , c] = retta(x1,y1,x2,y2)

% a = y1-y2;
% b = x2-x1;
% c = y1*(x1-x2)+x1*(y2-y1);

m = length(x1);% numero di centri della telecamera
% calcolo il numero solo su x poichè anche il numero di y sarà uguale
a = zeros(1,m);
b = zeros(1,m);
c = zeros(1,m);

for i = 1 : m
    a(i) = y1(i) - y2;
    b(i) = x2 - x1(i);
    c(i) = y1(i)*(x1(i) - x2) + x1(i)*(y2 - y1(i));
end