%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= %f', s_x);
s_y = sum(y_tab);
fprintf('\n s_y = %f', s_y);
s_xy = sum(x_tab .* y_tab);
fprintf('\n s_produtos = %f', s_xy);
s_x2 = sum(x_tab.^2);
fprintf('\n s_square_x = %f', s_x2);
s_y2 = sum(y_tab.^2);
fprintf('\n s_square_y = %f', s_y2 );

%Coeficientes da reta ajustada
b_1 = (s_x * s_y - n * s_xy)/(s_x^2 - n *s_x2);
b_0 = (s_y - b_1 * s_x)/n;
fprintf('\n b_0 = %f', b_0);
fprintf('\n b_1 %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