function F = us_block_steady_state(x)
    % Extract variables from input vector x
    cchi_w = x(1);
	k_w = x(2)
	s_inv_w = x(3);
	inv_w = x(4);
	lambdda_w = x(5);
	c_w = x(6);
	w_w = x(7);
	h_w = x(8);
	B_w = x(9);
	mmu_w = x(10);
	rk_w = x(11);
	Ak_w = x(12);
	r_w = x(13);
	R_w = x(14);
	pi_w = x(15);
	y_w = x(16);
	z_w = x(17);
	A_w = x(18);
	xxi_w = x(19);
	pi_target_w = x(20); 
	r_w_n = x(21);
	thetta_w = x(22);
	ggamma_A_w = x(23);
	ggamma_Ak_w = x(24);
	ggamma_B_w = x(25);
    % Parameters (example values, adjust as needed)
	v_w							= 1.5;
	deltta_w					= 0.02;
	eppsilon_w					= 6;
	betta_w						= 0.9869;
	ggammac_w					= 0.619;
	alppha_w					= 0.251;
	ommega_w					= 1.272;
	kapppa_w					= 13.130;
	ggamma_w					= 0.489;
	pphi_R_w					= 0.596;
	pphi_pi_w					= 2.696;
	pphi_y_w					= 0.251;
	rrho_a_w					= 0.271;
	rrho_ak_w					= 0.557;
	rrho_b_w					= 0.065;
	rrho_cchi_w					= 0.556;
	rrho_pi_w					= 0.934;
	rrho_i_w					= 0.651;
	rrho_z_w					= 0.748;
	rrho_thetta_w				= 0.685;
	sigmma_a_w					= 0.018;
	sigmma_ak_w					= 0.012;
	sigmma_b_w					= 0.031;
	sigmma_cchi_w				= 0.006;
	sigmma_r_w					= 0.015;
	sigmma_eppsilon_w			= 0.009;
	sigmma_pi_w					= 0.005;
	sigmma_i_w					= 0.095;
	sigmma_z_w					= 0.009;
	sigmma_thetta_w				= 0.007;
    steady_state = @(x) x; % Simplified steady-state assumption

    % Define equations as residuals
    F(1)  = cchi_w - ((1 - rrho_cchi_w) * steady_state(cchi_w) + rrho_cchi_w * cchi_w);
    F(2)  = k_w - (( (1 - deltta_w)* k_w) + s_inv_w * (1 - (ommega_w / 2) * ( (inv_w / inv_w) -1)^2)* inv_w);
    F(3)  = s_inv_w - ((1 - rrho_i_w) * steady_state(s_inv_w) + rrho_i_w * s_inv_w);
    F(4)  = lambdda_w - (cchi_w / (c_w - (ggammac_w * c_w ) ) ) ;
    F(5)  = lambdda_w * w_w - (cchi_w* ( (h_w / B_w)^(1 / v_w) ) );
    F(6)  = mmu_w - (betta_w * ( lambdda_w * rk_w + (1 - deltta_w) * mmu_w) );
    F(7)  = lambdda_w / Ak_w - ( mmu_w * s_inv_w * (1 - (ommega_w / 2) * ( ( (inv_w / inv_w ) - 1)^2) - ommega_w * ( ( inv_w / inv_w ) -1 ) * ( inv_w / inv_w ) ) + betta_w * ( mmu_w * s_inv_w * ommega_w *  ( ( inv_w / inv_w ) -1 )* ( inv_w / inv_w )^2) );
    F(8)  = lambdda_w - ( betta_w * (r_w * lambdda_w ) );
    F(9)  = r_w - ( R_w / pi_w );
    F(10) = y_w - ( z_w * (k_w ^ alppha_w) * (A_w * h_w) ^ (1 - alppha_w) );
    F(11) = z_w - ( (1 - rrho_z_w) * steady_state(z_w) + rrho_z_w * z_w );
    F(12) = lambdda_w * w_w - ( ( ( 1 - alppha_w ) * ( xxi_w * ( y_w / h_w ) ) ) );
    F(13) = lambdda_w * rk_w - ( ( alppha_w * xxi_w * ( y_w / k_w ) ) );
    F(14) = 1 + ( xxi_w / lambdda_w ) - kapppa_w * ( pi_w / ( ( pi_target_w ^ ggamma_w ) * (pi_w ^ ( 1 - ggamma_w ) ) ) - 1 ) * ( pi_w / ( ( pi_target_w ^ ggamma_w ) * ( pi_w ^ (1 - ggamma_w ) ) ) ) + betta_w * kapppa_w * ( lambdda_w / lambdda_w ) * ( pi_w / ( ( pi_target_w ^ ggamma_w ) * ( pi_w ^ (1 - ggamma_w) ) ) - 1 ) * ( pi_w / ( ( pi_target_w ^ ggamma_w ) * ( pi_w ^ ( 1 - ggamma_w ) ) ) * ( y_w / y_w ) );
    F(15) = R_w - ( ( R_w^pphi_R_w ) * ( ( steady_state(R_w) * ( ( ( pi_w / pi_target_w ) ^ ( pphi_pi_w ) ) * ( ( y_w / steady_state(y_w) ) ^ ( pphi_y_w ) ) ) ) ^ ( 1 - pphi_R_w ) ) ) ;
    F(16) = pi_target_w - ( ( 1 - rrho_pi_w ) * steady_state(pi_w) + rrho_pi_w * pi_target_w );
    F(17) = lambdda_w - ( betta_w * r_w_n * lambdda_w * thetta_w );
    F(18) = thetta_w - ( ( 1 - rrho_thetta_w ) * steady_state(thetta_w) + rrho_thetta_w * thetta_w );
    F(19) = y_w - ( c_w + inv_w + ( kapppa_w /2 ) * ( ( ( ( pi_w / ( ( pi_target_w ^ ggamma_w ) * ( pi_w ^ ( 1 - ggamma_w ) ) ) ) -1 ) ) ^ 2 ) * y_w );
    F(20) = log(A_w) - ( log(A_w) + ggamma_A_w );
    F(21) = log(Ak_w) - ( log(Ak_w) + ggamma_Ak_w );
    F(22) = log(B_w) - ( log(B_w) + ggamma_B_w );
    F(23) = ggamma_A_w - ( ( 1 - rrho_a_w ) * steady_state(ggamma_A_w) + rrho_a_w * ggamma_A_w );
    F(24) = ggamma_Ak_w - ( ( 1 - rrho_ak_w ) * steady_state(ggamma_Ak_w) + rrho_ak_w * ggamma_Ak_w );
    F(25) = ggamma_B_w - ( ( 1 - rrho_b_w ) * steady_state(ggamma_B_w) + rrho_b_w * ggamma_B_w ) ;
end
