% Numero di valori
n = 50 ;

% Valori variabile indipendente
x = linspace(0, 1, n)' ;

% Valori variabile dipendente, misura
c1 = 3.0 ;
c2 = 9.0 ;
r  = 0.2 ;
y  = misura_con_rumore(x, c1, c2, r) ;

% Minimi quadrati
A = [ ones(n,1) x ] ;
c = A\y ;

% Plot
plot( x , y , 'ob' ) ; % dati
hold on ;
plot( x , c(1)+c(2)*x , '-r' ) ; % retta di regressione

% Output testuale dei coefficienti
testo = [ 'Regressione lineare di %d valori:\n' ...
          ' c      = (%.2f,%.2f)\n' ...
          ' c_regr = (%.2f,%.2f)\n' ] ;
fprintf(testo, n,         ...
               c1  , c2 , ...
               c(1), c(2) ) ;
           
% Minimi quadrati espliciti
c_lsq = (A'*A)\(A'*y) ;