clear('x')
r=1.01;
x(1)=0.1;
for n=1:29
    x(n+1)=r*x(n)
end
plot([1:1:30],x,'-o')
