
t = [time variables]';
p = [values        ]';

n = max(size(t));

for j=1:n
A(:,j) = t.^(n-j);
end

s = A\p;

plot(t',p','o');
hold on
x = cumsum(ones(1,timerange));
plot(x,polyval(s,x),'g');

