% Input
xy = zeros(3,2) ;
caratteri_coord = 'xy' ;  
for punto = 1 : 3
  for coord =  1 : 2  
    messaggio = sprintf('Inserire la coordinata %c del punto %d: ', caratteri_coord(coord), punto) ;  
    xy(punto,coord) = input( messaggio ) ;
  end
end

% Coefficienti
x = xy(:,1) ;  
e = 0 : 2 ;  
V = x .^ e ;
y = xy(:,2) ;
a = V\y ;

% Output coefficienti
fprintf( 'a%d = %.2f\n' , [ 1:3 ; a' ] ) ;
%fprintf( 'y =' ) ;
%fprintf( ' %+.2fx^%d ' , [ a' ; 0:2 ] ) ;
%fprintf( '\n' ) ;

% Plot
x_min  = min(x) ;
x_max  = max(x) ;  
x_med  = ( x_min + x_max ) / 2 ;  
delta  = x_max - x_min ;  
x_plot = linspace( x_med-delta , x_med+delta , 1000 )' ; % vettore colonna
y_plot = ( x_plot .^ e ) * a ;
plot( x_plot , y_plot , 'LineWidth' , 1 ) ;
hold on ;
plot( x , y , 'ro' ) ;