read in the data
center the x-values x - xbar
find the best linear fit b0 + b1 x
	X = [ 1 x ] y
find the best quadratic fit c0 + c1 x + c2 x^2  
	X = [ 1 x x^2] y
output 
[xbar;
b0;
b1;
c0;
c1;
c2]
-->proj0('sample1')                                               
 ans  =
    3.   
    15.  
    7.   
    13.  
    7.   
    1.   
-->proj0('sample2')
 ans  =
    50.5       
    0.6302007  
    3.416E-19  
    0.9794065  
    3.416E-19  
  - 0.0004191