function f = system1(x, alpha, betta,  gamma_ss, delta, sgma, omg, b_ss, g, n_ss , x_ss)
 
f = [x(1) - 1 + betta^sgma*(x(3)*x(11))^(sgma - 1)
    x(2) - (1 - gamma_ss*betta^sgma*x(11)^(sgma-1))/x(1)
    x(3) - omg - (1 - omg)*x(2)^(1/(1 - sgma))
    x(4) - ((1 - alpha)/(x(11) - (1 - delta)))^(1/alpha)
    x(5) - x(4)^(1 - alpha)
    x(6)*(1 - omg*(1 + x_ss)/(x(3)*x(11))) - (alpha - x(10))*x(5)
    x(7) - x(1)*((1 + (x(2) - 1)*x(8))*x(11)*x(9) + x(6))
    x(8) - (1 - omg)/(1 - omg*(1 - x(1)*x(2))*x(11)/(1 + x_ss + n_ss))
    x(9) - x(4) - b_ss*x(5)
    x(10) - g - (x(11) - (1 + x_ss + n_ss))*b_ss
    (1 - g)*x(5) - x(7) - (x_ss + n_ss + delta)*x(4)];

% x = [pi_ss epsilon_ss omega_ss k_ss y_ss h_ss c_ss lambda_ss a_ss tau_ss
% R_ss]