a=0.5;
b=2;

N=10;

X=b;

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

figure(1)
plot(0:N,X)

figure(2)
plot(xosa,pom)