irflength           = 5000;
irfburninlength     = 0;

% oo_.stochastic_steady_state = sss(oo_.dr, irfburninlength, options_.order);

for ivariable = 1:M_.endo_nbr
   dum_name = [char(M_.endo_names(ivariable,:)), '_ss'];
    assignin('base',dum_name,oo_.steady_state(ivariable,:)');                     
end


for ishock = 1:M_.exo_nbr    
    variables_irfsss  =  irfss(oo_.dr,sqrt(M_.Sigma_e(:,ishock)),irflength,irfburninlength,options_.order, M_, options_)';
    for ivariable = 1:M_.endo_nbr
        dum_name2 = [char(M_.endo_names(ivariable,:)), '_', char(M_.exo_names(ishock,:)), '_', text];
        assignin('base',dum_name2,variables_irfsss(:,ivariable));                         
    end
end