var f; 
varexo  a; 
parameters b; 
b=0.7;

model(block, bytecode); 
f=a+b*f(-1);
end; 

initval; 
a = 0.50000000;
f = 1.66666667;
end; 

steady(solve_algo=5);

shocks;
var a;
periods 1:20;
values 0.8;
end;

  
simul(periods=20, stack_solve_algo=5);

plot(f);