% Métodos Numéricos - Engs. Civil e Electrotécnica 
% exercico 1 da ficha 4
% Carlos Balsa - 2009/2010

t=[1 2 3 4 5 6 7]';
y=[0.5 2.5 2 4 3.5 6 5.5]';

% alínea a)
m=length(t);

A=[t.^0 t];

H=A'*T
b=A'*y
x=H\b
z=x(1)+x(2)*t; % eq. da recta de regressão

% alínea b)
t1=2.23;
z1=x(1)+x(2)*t1
t2=4.5;
z2=x(1)+x(2)*t2

% alínea c)
ybar=mean(y);
St=norm(y-ybar,2)^2;
Sr=norm(y-A*x,2)^2;
r=sqrt((St-Sr)/St)

% alínea d)
clf
plot(t,y,'o'); hold on;
plot(t,z,'-r');
