% Oc5_1b_RegressaoLinearAjusteParabola
clc;
% Regressão linear com reta como modelo
x=[0.25; 0.30; 0.36; 0.40; 0.50; 0.60; 0.70];
y=[96; 60; 35; 24; 10; 5; 3];
plot(x, y, 'or');
hold on;
F=[ones(length(x),1)   x   1./x];
p=inv(F'*F)*(F'*y);
residual=norm(F*p-y);
xn=[0.25:0.01:0.70];
yn=p(1)+ p(2)*xn + p(3)./xn;
plot(xn, yn);
axis([0.2 1.0 0 100]);