function [ys,check] = paper_0_1_steadystate(ys,exe)
    global M_


    % Import variables
    NumberOfParameters = M_.param_nbr;                            % Number of deep parameters.
    for j = 1:NumberOfParameters                                  % Loop...
      paramname = deblank(M_.param_names(j,:));                   %    Get the name of parameter i. 
      eval([ paramname ' = M_.params(' int2str(j) ');']);         %    Get the value of parameter i.
    end                                                           % End of the loop.  
    check = 0;


    % Compute Steady State Values using K/L system
    x0 = [1,0.1,1,0.1,1];
    optiops2 = optimset('Display','off');
    [ss,fval,exitflag,output] = fsolve(@ss_solver1,x0,optiops2);


    % ss = fsolve(@tester,x0,optiops);
    optiops1 = optimset('Display', 'off');
    ss2 = fsolve(@ss_solver2,ss,optiops1);
    
    % ss_kl = fsolve(@(x) main5b_ss_kl(x,M_), init, tester);

    % Back out other SS variables from k*, l*
    k = ss2(1);
    l = ss2(2);
    i = ss2(3);
    c = ss2(4);
    w = ss2(5);

    % var_z = var0z/(1-rvarz);
    % z = z0/(1-rz);
    var_bh = var_bh_0/(1-r_varbh);
    bh = bh0/(1-rbh);
    r = (1/bh)-1+d;
    
    y = z * k^a * l^(1-a);
    pr = z * k^a * l^(1-a)-w*l-r*k-g1*i^g2*k^(1-g2);
    % cost = g1*(i^g2)*k^(1-g2);
    


    % Load SS back into dynare    
    NumberOfEndogenousVariables = M_.endo_nbr;                    % Number of endogenous variables.
    ys = zeros(NumberOfEndogenousVariables,1);                    % Initialization of ys (steady state).
    for j = 1:NumberOfEndogenousVariables                         % Loop...
      varname = deblank(M_.endo_names(j,:));                      %    Get the name of endogenous variable i.                     
      eval(['ys(' int2str(j) ') = ' varname ';']);                %    Get the steady state value of this variable.
    end   
    
end

% old code
    % Compute Steady State Values using K/L system
    %init = [2; 2];
    % ss_kl = fsolve(@(x) main5b_ss_kl(x,M_), init, tester);
    %ss_kl = fsolve(@(x) paper_0_1_ss_kl(x,M_), init, tester);
    % Back out other SS variables from k*, l*