x = 0:0.1:3
y1 = [];
y2 = [];

%f = 1/(a+1);

[m n] = size(x);
for k = 1:1:n
 y1 = [y1,1/(x(:,k)+1)];
 y2 = [y2,x(:,k)]; 
end;

a = [2];
b = [0];

for i = 1:1:10
a = [a,1/(a(:,i)+1)];
b = [b,1/(a(:,i)+1)]
end;


plot(x,y1,'r',x,y2,'b',a,b,'y')
%plot(x,y2,'r')