beta calculado como (X'X)^(-1) X'y:
beta_matricial (3 x 1)

      4,4307 
      2,1139 
    -0,98078 


beta calculado por 'ols' de Gretl:
beta_ols (3 x 1)

      4,4307 
      2,1139 
    -0,98078 


Diferencia (debe ser practicamente cero):
diferencia (3 x 1)

  2,6645e-15 
  4,4409e-16 
 -1,3323e-15