disp('ESCUELA SUPERIOR POLITECNICA DE CHIMBORAZO')
disp('Regresion Polinomial')
disp('Daniela Zuñiga')
clear
clc
x=input('Ingrese lo datos de x: ');
y=input('Ingrese lo datos de y: ');
dim=length(y)
sumx=0;st=0;sumy=0;sr=0;
sumx2=0;sumx3=0;sumx4=0;sumxy=0;sumx2y=0;
for i=1:dim
sumx=sumx+x(i);
sumy=sumy+y(i);
sumx2=sumx2+x(i)*x(i);
sumx3=sumx3+x(i)*x(i)*x(i);
sumx4=sumx4+x(i)*x(i)*x(i)*x(i);
sumxy=sumxy+x(i)*y(i);
sumx2y=sumx2y+x(i)*x(i)*y(i);
end
xm=sumx/dim;
ym=sumy/dim;
a=[dim sumx sumx2;sumx sumx2 sumx3;sumx2 sumx3 sumx4]
b=[sumy;sumxy;sumx2y]
[m,n]=size(a);
i=1;
l(i,i)=sqrt(a(i,i));
for j=1:m
l(i,j)=a(i,j)/l(i,i);
end
i=i+1;
