x=[0:0.01:10];
n=size(x,2);
k=floor(n/100);
y=zeros(n,1);
for i=0:k-1
  c1=f(x(100*i+1));
  c2=f(x(100*(i+1)+1));
  for j=1:101
    y(100*i+j)=c1+((j-1)/100)*(c2-c1);
  end;
end;
plot(x,y,'black');
axis([-1 11 -2 3]);
axis on;
