function SteadyState
 
b=0.991; % Parameter
v=0.357;
ceta=0.3;
delta=0.0196;

x0=[1 1 1 1 1 0]; % Startwerte

R = fsolve(@eqn,x0)
 
    function F = eqn(x)
        % V_ss=x(1)
        % c_ss=x(2)
        % l_ss=x(3)
        % k_ss=x(4)
        % R_ss=x(5)
        % i_ss=x(6)
        F=[ ((x(2)^(v))*((1 - x(3))^(1-v)))- x(1);
            ((1/b)) - (ceta * (x(4)^(ceta - 1)) * (x(3)^(1-ceta)) + 1 - delta);
            ((1-ceta) * (x(4)^(ceta)) * (x(3)^(-ceta))) - (((1-v)/v) * (x(2) / (1 - x(3))));
            ((1/b)) - (x(5));
            ((x(4)^(ceta)) * (x(3)^(1-ceta))) - (x(2)+ x(6));
            (delta * x(4)) - (x(6));          
            ];
    end
 
end