clear all

b=2;
a=0.5;
N=10;

X=0;
%X=b/(1-a)
xosa=0;
pomX=X;
for i=1:N
    xosa=[xosa, i];
    pomX=[pomX, a*X(i)];
    X(i+1)=a*X(i)+b;
    xosa=[xosa, i];
    pomX=[pomX, X(i+1)];
end

figure(1)
plot(0:N,X)
figure(2)
plot(xosa,pomX)