function [Y, U, N] = my_function2(TAU, A, Pi, R, I, realinterest, Pi_star, S, nu, siggma, varphi, phi_pi, phi_y, theta, rho_nu, rho_a, betta, alppha, epsilon, rho_tau, psi_1, eta_1, eta_2, psi_2)
    % fsolve 옵션 설정
    options = optimoptions('fsolve', 'MaxFunEvals', 10000, 'MaxIter', 10000, 'TolFun', 1e-8, 'TolX', 1e-8);

    % 초기값 설정
    X0 = [1, 0.3, 0.3];

    % fsolve로 해 구하기
    [XX, ~, exitflag] = fsolve(@(XX) equations_to_solve(XX, TAU, A, S, epsilon, psi_1, eta_1, eta_2, varphi), X0, options);

    % fsolve가 성공적으로 해를 구했는지 확인
    if exitflag <= 0
        warning('fsolve did not converge.');
    end

    % 해 결과를 출력
    Y = XX(1);
    U = XX(2);
    N = XX(3);

    % 내부 방정식 함수 정의
    function F = equations_to_solve(ZZ, TAU, A, S, epsilon, psi_1, eta_1, eta_2, varphi)
        Y = ZZ(1);
        U = ZZ(2);
        N = ZZ(3);

        % 연립방정식 정의
        F = [
            Y - A * N / S;
            U - (psi_1 * TAU * Y^(-psi_2) / (eta_1 * eta_2))^(1 / (eta_2 - 1));
            N - ((epsilon - 1) / epsilon - psi_1 * TAU * (1 - U) * (1 - psi_2) * Y^(-psi_2) - eta_1 * U^eta_2)^(1 / (varphi + 1)) * (S / (1 - eta_1 * U^eta_2 * S))^(1 / (varphi + 1))
        ];
    end
end
