Páginas

mìnimos cuadrados

clc
clear all
disp('datos ingresados')
x=[0 2 3 6 7]
y=[0.12 0.153 0.170 0.225 0.260]
n=length(x);
a=[n sum(x)
sum(x) sum(x.^2)];
b=[sum(y) sum(x.*y)]';
c=inv(a)*b;
disp('la ecuacion es:')
fprintf(' %5.4f *x + %5.4f\n ',c(2),c(1))
%comprobacion
polyfit(x,y,1);