fprintf("Ajuste de Curvas por Mínimos Quadrados\n\n");

%Dados Iniciais 
x_tabs=[0.3 2.7 4.5 5.9 7.8];
y_tab=[1.8 1.9 3.1 3.9 3.3];

n=length(x_tab);
p=2;

S_x=sum(x_tab);

S_y=sum(y_tab);

S_xy=sum(x_tab .* y_tab);

S_xx= sum(x_tab .* x_tab);

S_yy= sum(y_tab .* y_tab);

%Calculo dos coeficientes
b1= (S_xx*S_yy - n*S_xy)/((S_x)^2 - n*S_xx)

b0= (S_y-b1*S_x)/n;

%Calculo do desvil Residual
  u_tab= b0 + b1*x_tab;

  D= sum((y_tab - u_tab) .^2 );

  %Indicadores Estatísticos de Qualidade de Regressão

  r2 = 1- (D/(S_yy-(1/n)*S_yy^2));
  sigma2 = D/(n-p);

  u=@(x)*b0+b1*x;

  x_grid=min(x_tab)-0.5:0.1:max(x_tab)+0.5;

   figure(1)

   plot(x_grid, u_grid);