% Ajuste de l�nea recta por m�nimos cuadrados 494

n = input('Escriba el n�mero de datos de la muestra: ');
sum_x = input('Escriba la sumatoria de los valores de X: ');
sum_y = input('Escriba la sumatoria de los valores de Y: ');
med_x = input('Escriba la media de los valores de X: ');
med_y = input('Escriba la media de los valores de Y: ');
sum_x_cuadrado = input('Escriba el valor de la sumator�a de los datos de X al cuadrado: ');
sum_x_y = input('Escriba el valor de la sumator�a de X y Y: ');

%Determinar a1
a1 = (n*sum_x_y - sum_x * sum_y)/(n*sum_x_cuadrado - sum_x^2);

%Determinar a0
a0 = med_y - a1*med_x;

equ = strcat(num2str(a0),' +  ', num2str(a1), 'x');

x = linspace(0,n);

xp=[1 2 3 4 5 6 7];
y=[0.5 2.5 2 4 2.4 6 5.5];
scatter(xp, y);
hold on
plot(x, a0 + a1*x);
hold off


disp(strcat('El ajuste por m�nimos cuadrados es: ', equ));