%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
u_tab = b_0 + B_1 * x_tab;
fprintf('\n u_tab %f', u_tab);
D = sum((y_tab - u_tab).^2);
fprintf('\n D %f', D);
r2 = 1 - (D/(s_y2 - (1/n)* S_y2));
fprintf('\n r2  %f', r2);