x=[0:0.01:10];
n=size(x,2);
y=zeros(n,1);
alpha=zeros(10,1);
alpha(1)=2*f(0);
for k=1:9
  alpha(k+1)=2*f(k)-alpha(k);
end;
for i=1:n
  for j=1:10
    y(i)=y(i)+alpha(j)*f1(x(i)+2-j);
  end;
end;
plot(x,y,'black');
axis([-1 11 -2 3]);
axis on;
