% Métodos Numéricos - Interpolação Polinomial
% Método da base monómica
% Exercicio 3 da ficha nº 8
% Carlos Balsa 2007/2008

% Resolução da alínea a)
clear all
t=[1 2 3 4]';
y=[11 29 65 125]';
c=ones(4,1);
A=[c t t.^2 t.^3];
x=A\y
tt=[1:0.1:4]';
p4=x(1)+x(2)*tt+x(3)*tt.^2+x(4)*tt.^3;

plot(t,y,'*r'); hold('on')
plot(tt,p4)

% Alínea d)

t=3.5;
p4=x(1)+x(2)*t+x(3)*t^2+x(4)*t^3