function output=affine_plot(P0,r,b,n)

output=zeros(1,n+1);
output(1)=P0;
for i =1:n
    output(i+1)=r*output(i)+b;
end

