%Luiz Henrique Meira Andrade Leite
%Pratica Minimos Quadrados
clc;

fprintf('Ajuste de Curva por Mínimos Quadrados \n \n');

%Entradas
x_tab = [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;

%Somatorios
s_x = sum(x_tab);
fprintf('\n s_x= %2.f', s_x);
s_y = sum(y_tab);
fprintf('\n s_y = %2.f', s_y);
s_produtos = sum(x_tab .* y_tab);
fprintf('\n s_produtos = %2.f', s_produtos);
s_square_x = (sum((s_x)).^2);
fprintf('\n s_square_x = %2.f', s_square_x);
s_square_y = (sum((s_y)).^2);
fprintf('\n s_square_y = %2.f', s_square_y );

%Coeficientes da reta ajustada
b_1 = (s_x .* s_y - n .* s_y)/(s_x.^2) - n .* (s_x .^2);
b_0 = s_y - (b_1 .* s_x)/n;
fprintf('\n b_0 = %2.f', b_0);
fprintf('\n b_1 %2.f', b_1);
%Desvio Residual
%for (i = 0;i = len(x_tab);i++) 
 %   u_i = b_0 + b_1 * x[i];
  %  fprintf(u_i)
    
   %$ end