One issue seems to be
for example in
log(A/steady_state(A))=rho_A*log(A(-1)/steady_state(A))+sigma_A*e_A;