function  [wel,n,w_c_out,w_u_out,r_cc_out,r_cu_out,r_uu_out,r_uc_out,X_c_out,X_u_out,l_cc_out,l_uu_out,l_cu_out,l_uc_out] ...
        = func_eqm_update(m,w_c,w_u,r_cc,r_cu,r_uu,r_uc,X_c,X_u,l_cc,l_uu,l_cu,l_uc)

%% 不同数据中间商市场结构下的均衡求解
switch m.MarketStruct
%%%%%%%%%%%%%%%%%%%%%%%%% Case1: Monopoly Competition%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    case 'competition'              
    
    % 商品生产者的成本、定价与贸易份额
    R_c = (m.chi_cc.*r_cc.^(1-m.o)+m.chi_cu.*r_cu.^(1-m.o)).^(1./(1-m.o));
    R_u = (m.chi_uu.*r_uu.^(1-m.o)+m.chi_uc.*r_uc.^(1-m.o)).^(1./(1-m.o));
      
    A_c = m.beta_c.^(-m.beta_c).*(1-m.beta_c).^(-(1-m.beta_c));
    A_u = m.beta_u.^(-m.beta_u).*(1-m.beta_u).^(-(1-m.beta_u));
    C_c = A_c.*R_c.^m.beta_c.*w_c.^(1-m.beta_c);
    C_u = A_u.*R_u.^m.beta_u.*w_u.^(1-m.beta_u);
    
    kappa = gamma((1-m.sigma+m.theta)./m.theta);
    kappa = m.sigma./(m.sigma-1).*kappa.^(1./(1-m.sigma));
    temp_cc = m.T_c.*C_c.^(-m.theta);
    temp_uu = m.T_u.*C_u.^(-m.theta);
    temp_cu = m.T_u.*(C_u.*m.tau_cu).^(-m.theta);
    temp_uc = m.T_c.*(C_c.*m.tau_uc).^(-m.theta);
    P_c = kappa.*(temp_cc+temp_cu).^(-1./m.theta);
    P_u = kappa.*(temp_uu+temp_uc).^(-1./m.theta);
    
    pi_cc = temp_cc./(temp_cc+temp_cu);
    pi_uu = temp_uu./(temp_uu+temp_uc);
    pi_cu = temp_cu./(temp_cc+temp_cu);
    pi_uc = temp_uc./(temp_uu+temp_uc);
    
    % 跨境数据流动政策
    switch m.DataPolicy
        case 'prohibit'          % 禁止原始数据流动
            D_c = X_c/P_c; 
            D_u = X_u/P_u;
            DShare_c = X_c/P_c;
            DShare_u = X_u/P_u;
        case 'free'              % 原始数据自由流动
            D_c = X_c/P_c + m.rho_u.*X_u/P_u;
            D_u = X_u/P_u + m.rho_c.*X_c/P_c;
            DShare_c = X_c/P_c;
            DShare_u = X_u/P_u;
    end

    % 数据市场出清
    RPhi_c = m.beta_c.*(m.sigma-1)./m.sigma.*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc));
    RPhi_u = m.beta_u.*(m.sigma-1)./m.sigma.*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu));
    rphi_cc = m.chi_cc.*(r_cc./R_c).^(1-m.o).*RPhi_c;
    rphi_cu = m.chi_cu.*(r_cu./R_c).^(1-m.o).*RPhi_c;
    rphi_uu = m.chi_uu.*(r_uu./R_u).^(1-m.o).*RPhi_u;
    rphi_uc = m.chi_uc.*(r_uc./R_u).^(1-m.o).*RPhi_u;

    % 加工数据量
    phi_cc = m.B_c.*D_c.^m.gamma_c.*l_cc.^(1-m.gamma_c);
    phi_uc = m.B_c.*D_c.^m.gamma_c.*l_uc.^(1-m.gamma_c);
    phi_uu = m.B_u.*D_u.^m.gamma_u.*l_uu.^(1-m.gamma_u);
    phi_cu = m.B_u.*D_u.^m.gamma_u.*l_cu.^(1-m.gamma_u);

    % 数据要素价格
    r_cc_out = rphi_cc./phi_cc;    
    r_cu_out = rphi_cu./phi_cu;
    r_uu_out = rphi_uu./phi_uu;
    r_uc_out = rphi_uc./phi_uc;  
    
    % 数据中间商市场势力
    epsilon_cc = m.o./(m.o-1);
    epsilon_uu = m.o./(m.o-1);
    epsilon_cu = m.o./(m.o-1);
    epsilon_uc = m.o./(m.o-1);

    % 数据部门劳动力分配
    l_cc_out = 1./epsilon_cc.*(1-m.gamma_c).*rphi_cc./w_c;
    l_uc_out = 1./epsilon_uc.*(1-m.gamma_c).*rphi_uc./w_c;
    l_cu_out = 1./epsilon_cu.*(1-m.gamma_u).*rphi_cu./w_u;
    l_uu_out = 1./epsilon_uu.*(1-m.gamma_u).*rphi_uu./w_u;

    % 劳动力市场出清
    wL_c = (m.sigma-1)./m.sigma.*(1-m.beta_c).*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc)) ...
           + (1-m.gamma_c).*(rphi_cc./epsilon_cc+rphi_uc./epsilon_uc) ...
           + pi_cc.*m.f_d.*w_c + pi_uc.*m.f_x.*w_c;
    wL_u = (m.sigma-1)./m.sigma.*(1-m.beta_u).*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu)) ...
           + (1-m.gamma_u).*(rphi_uu./epsilon_uu+rphi_cu./epsilon_cu) ...
           + pi_uu.*m.f_d.*w_u + pi_cu.*m.f_x.*w_u;
    w_c_out = wL_c./m.L_c;
    w_u_out = wL_u./m.L_u;

    % 商品市场出清
    Pi_c = 1./m.sigma./(1+m.t_uc).*pi_uc.*X_u - pi_uc.*m.f_x.*w_c ...
           + 1./m.sigma.*pi_cc.*X_c - pi_cc.*m.f_d.*w_c;
    Pi_u = 1./m.sigma./(1+m.t_cu).*pi_cu.*X_c - pi_cu.*m.f_x.*w_u ...
           + 1./m.sigma.*pi_uu.*X_u - pi_uu.*m.f_d.*w_u;
    M_c = (epsilon_cc-1)./epsilon_cc.*rphi_cc + (epsilon_uc-1)./epsilon_uc.*rphi_uc ...
          + m.gamma_c.*(1./epsilon_cc.*rphi_cc+1./epsilon_uc.*rphi_uc);
    M_u = (epsilon_uu-1)./epsilon_uu.*rphi_uu + (epsilon_cu-1)./epsilon_cu.*rphi_cu ...
          + m.gamma_u.*(1./epsilon_uu.*rphi_uu+1./epsilon_cu.*rphi_cu);
    tariff_c = m.t_cu./(1+m.t_cu).*pi_cu.*X_c;
    tariff_u = m.t_uc./(1+m.t_uc).*pi_uc.*X_u;
    upsilon_u = m.nu.*X_u;

    X_c_out = wL_c + Pi_c + M_c + tariff_c - upsilon_u;
    X_u_out = wL_u + Pi_u + M_u + tariff_u + upsilon_u;
      
    
    %%%%%%%%%%%%%%%%%%%%%%%%% Case2: Oligopoly Competition %%%%%%%%%%%%%%%%%%%%%%%%%%%%
    case 'oligopoly'            
    
    % 商品生产者的成本、定价与贸易份额
    R_c = (m.chi_cc.*r_cc.^(1-m.o)+m.chi_cu.*r_cu.^(1-m.o)).^(1./(1-m.o));
    R_u = (m.chi_uu.*r_uu.^(1-m.o)+m.chi_uc.*r_uc.^(1-m.o)).^(1./(1-m.o));
      
    A_c = m.beta_c.^(-m.beta_c).*(1-m.beta_c).^(-(1-m.beta_c));
    A_u = m.beta_u.^(-m.beta_u).*(1-m.beta_u).^(-(1-m.beta_u));
    C_c = A_c.*R_c.^m.beta_c.*w_c.^(1-m.beta_c);
    C_u = A_u.*R_u.^m.beta_u.*w_u.^(1-m.beta_u);
    
    kappa = gamma((1-m.sigma+m.theta)./m.theta);
    kappa = m.sigma./(m.sigma-1).*kappa.^(1./(1-m.sigma));
    temp_cc = m.T_c.*C_c.^(-m.theta);
    temp_uu = m.T_u.*C_u.^(-m.theta);
    temp_cu = m.T_u.*(C_u.*m.tau_cu).^(-m.theta);
    temp_uc = m.T_c.*(C_c.*m.tau_uc).^(-m.theta);
    P_c = kappa.*(temp_cc+temp_cu).^(-1./m.theta);
    P_u = kappa.*(temp_uu+temp_uc).^(-1./m.theta);
    
    pi_cc = temp_cc./(temp_cc+temp_cu);
    pi_uu = temp_uu./(temp_uu+temp_uc);
    pi_cu = temp_cu./(temp_cc+temp_cu);
    pi_uc = temp_uc./(temp_uu+temp_uc);
    
    % 跨境数据流动政策
    switch m.DataPolicy
        case 'prohibit'          % 禁止原始数据流动
            D_c = X_c/P_c;
            D_u = X_u/P_u;
            DShare_c = X_c/P_c;
            DShare_u = X_u/P_u;
        case 'free'              % 原始数据自由流动
            D_c = X_c/P_c + m.rho_u.*X_u/P_u;
            D_u = X_u/P_u + m.rho_c.*X_c/P_c;
            DShare_c = X_c/P_c;
            DShare_u = X_u/P_u;
    end

    % 数据市场出清
    RPhi_c = m.beta_c.*(m.sigma-1)./m.sigma.*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc));
    RPhi_u = m.beta_u.*(m.sigma-1)./m.sigma.*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu));
    rphi_cc = m.chi_cc.*(r_cc./R_c).^(1-m.o).*RPhi_c;
    rphi_cu = m.chi_cu.*(r_cu./R_c).^(1-m.o).*RPhi_c;
    rphi_uu = m.chi_uu.*(r_uu./R_u).^(1-m.o).*RPhi_u;
    rphi_uc = m.chi_uc.*(r_uc./R_u).^(1-m.o).*RPhi_u;

    % 加工数据量
    phi_cc = m.B_c.*D_c.^m.gamma_c.*l_cc.^(1-m.gamma_c);
    phi_uc = m.B_c.*D_c.^m.gamma_c.*l_uc.^(1-m.gamma_c);
    phi_uu = m.B_u.*D_u.^m.gamma_u.*l_uu.^(1-m.gamma_u);
    phi_cu = m.B_u.*D_u.^m.gamma_u.*l_cu.^(1-m.gamma_u);

    % 数据要素价格
    r_cc_out = rphi_cc./phi_cc;    
    r_cu_out = rphi_cu./phi_cu;
    r_uu_out = rphi_uu./phi_uu;
    r_uc_out = rphi_uc./phi_uc;  

    % 数据中间商市场份额和市场势力
    s_cc = rphi_cc./(rphi_cc+rphi_cu);
    s_cu = rphi_cu./(rphi_cc+rphi_cu);
    s_uu = rphi_uu./(rphi_uu+rphi_uc);
    s_uc = rphi_uc./(rphi_uu+rphi_uc);
    
    epsilon_cc = (m.o-m.o.*s_cc)./(m.o-m.o.*s_cc-1);
    epsilon_uu = (m.o-m.o.*s_uu)./(m.o-m.o.*s_uu-1);
    epsilon_cu = (m.o-m.o.*s_cu)./(m.o-m.o.*s_cu-1);
    epsilon_uc = (m.o-m.o.*s_uc)./(m.o-m.o.*s_uc-1);

    % 数据部门劳动力分配
    l_cc_out = 1./epsilon_cc.*(1-m.gamma_c).*rphi_cc./w_c;
    l_uc_out = 1./epsilon_uc.*(1-m.gamma_c).*rphi_uc./w_c;
    l_cu_out = 1./epsilon_cu.*(1-m.gamma_u).*rphi_cu./w_u;
    l_uu_out = 1./epsilon_uu.*(1-m.gamma_u).*rphi_uu./w_u;

    % 劳动力市场出清
    wL_c = (m.sigma-1)./m.sigma.*(1-m.beta_c).*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc)) ...
           + (1-m.gamma_c).*(rphi_cc./epsilon_cc+rphi_uc./epsilon_uc) ...
           + pi_cc.*m.f_d.*w_c + pi_uc.*m.f_x.*w_c;
    wL_u = (m.sigma-1)./m.sigma.*(1-m.beta_u).*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu)) ...
           + (1-m.gamma_u).*(rphi_uu./epsilon_uu+rphi_cu./epsilon_cu) ...
           + pi_uu.*m.f_d.*w_u + pi_cu.*m.f_x.*w_u;
    w_c_out = wL_c./m.L_c;
    w_u_out = wL_u./m.L_u;

    % 商品市场出清
    Pi_c = 1./m.sigma./(1+m.t_uc).*pi_uc.*X_u - pi_uc.*m.f_x.*w_c ...
           + 1./m.sigma.*pi_cc.*X_c - pi_cc.*m.f_d.*w_c;
    Pi_u = 1./m.sigma./(1+m.t_cu).*pi_cu.*X_c - pi_cu.*m.f_x.*w_u ...
           + 1./m.sigma.*pi_uu.*X_u - pi_uu.*m.f_d.*w_u;
    M_c = (epsilon_cc-1)./epsilon_cc.*rphi_cc + (epsilon_uc-1)./epsilon_uc.*rphi_uc ...
          + m.gamma_c.*(1./epsilon_cc.*rphi_cc+1./epsilon_uc.*rphi_uc);
    M_u = (epsilon_uu-1)./epsilon_uu.*rphi_uu + (epsilon_cu-1)./epsilon_cu.*rphi_cu ...
          + m.gamma_u.*(1./epsilon_uu.*rphi_uu+1./epsilon_cu.*rphi_cu);
    tariff_c = m.t_cu./(1+m.t_cu).*pi_cu.*X_c;
    tariff_u = m.t_uc./(1+m.t_uc).*pi_uc.*X_u;
    upsilon_u = m.nu.*X_u;

    X_c_out = wL_c + Pi_c + M_c + tariff_c - upsilon_u;
    X_u_out = wL_u + Pi_u + M_u + tariff_u + upsilon_u;

    
    %%%%%%%%%%%%%%%%%%%%%%%%% Case3: Monopoly %%%%%%%%%%%%%%%%%%%%%%%%%%%%
    case 'monopoly'            
    
    % 商品生产者的成本、定价与贸易份额
    R_c = (m.chi_cu.*r_cu.^(1-m.o)).^(1./(1-m.o));
    R_u = (m.chi_uu.*r_uu.^(1-m.o)).^(1./(1-m.o));
      
    A_c = m.beta_c.^(-m.beta_c).*(1-m.beta_c).^(-(1-m.beta_c));
    A_u = m.beta_u.^(-m.beta_u).*(1-m.beta_u).^(-(1-m.beta_u));
    C_c = A_c.*R_c.^m.beta_c.*w_c.^(1-m.beta_c);
    C_u = A_u.*R_u.^m.beta_u.*w_u.^(1-m.beta_u);
    
    kappa = gamma((1-m.sigma+m.theta)./m.theta);
    kappa = m.sigma./(m.sigma-1).*kappa.^(1./(1-m.sigma));
    temp_cc = m.T_c.*C_c.^(-m.theta);
    temp_uu = m.T_u.*C_u.^(-m.theta);
    temp_cu = m.T_u.*(C_u.*m.tau_cu).^(-m.theta);
    temp_uc = m.T_c.*(C_c.*m.tau_uc).^(-m.theta);
    P_c = kappa.*(temp_cc+temp_cu).^(-1./m.theta);
    P_u = kappa.*(temp_uu+temp_uc).^(-1./m.theta);
    
    pi_cc = temp_cc./(temp_cc+temp_cu);
    pi_uu = temp_uu./(temp_uu+temp_uc);
    pi_cu = temp_cu./(temp_cc+temp_cu);
    pi_uc = temp_uc./(temp_uu+temp_uc);
    
    % 跨境数据流动政策
    switch m.DataPolicy
        case 'prohibit'          % 禁止原始数据流动
            D_c = 0;
            D_u = X_u/P_u;
            DShare_c = 0;
            DShare_u = X_u/P_u;
        case 'free'              % 原始数据自由流动
            D_c = 0;
            D_u = X_u/P_u + m.rho_c.*X_c/P_c;
            DShare_c = X_c/P_c;
            DShare_u = X_u/P_u;
    end

    % 数据市场出清
    RPhi_c = m.beta_c.*(m.sigma-1)./m.sigma.*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc));
    RPhi_u = m.beta_u.*(m.sigma-1)./m.sigma.*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu));
    rphi_cu = RPhi_c;
    rphi_uu = RPhi_u;

    % 加工数据量
    phi_uu = m.B_u.*D_u.^m.gamma_u.*l_uu.^(1-m.gamma_u);
    phi_cu = m.B_u.*D_u.^m.gamma_u.*l_cu.^(1-m.gamma_u);

    % 数据要素价格
    r_cc_out = 0;            % useless
    r_uc_out = 0;
    r_uu_out = rphi_uu./phi_uu;
    r_cu_out = rphi_cu./phi_cu;
    
    % 数据中间商市场势力
    epsilon_uu = m.mu_u;
    epsilon_cu = m.mu_c; 

    % 数据部门劳动力分配
    l_cc_out = 0;            % useless
    l_uc_out = 0;
    l_cu_out = 1./epsilon_cu.*(1-m.gamma_u).*rphi_cu./w_u;
    l_uu_out = 1./epsilon_cu.*(1-m.gamma_u).*rphi_uu./w_u;
    
    % 劳动力市场出清
    wL_c = (m.sigma-1)./m.sigma.*(1-m.beta_c).*(pi_cc.*X_c+pi_uc.*X_u./(1+m.t_uc)) ...
           + pi_cc.*m.f_d.*w_c + pi_uc.*m.f_x.*w_c;
    wL_u = (m.sigma-1)./m.sigma.*(1-m.beta_u).*(pi_uu.*X_u+pi_cu.*X_c./(1+m.t_cu)) ...
           + (1-m.gamma_u).*(rphi_uu./epsilon_uu+rphi_cu./epsilon_cu) ...
           + pi_uu.*m.f_d.*w_u + pi_cu.*m.f_x.*w_u;
    w_c_out = wL_c./m.L_c;
    w_u_out = wL_u./m.L_u;

    % 商品市场出清
    Pi_c = 1./m.sigma./(1+m.t_uc).*pi_uc.*X_u - pi_uc.*m.f_x.*w_c ...
           + 1./m.sigma.*pi_cc.*X_c - pi_cc.*m.f_d.*w_c;
    Pi_u = 1./m.sigma./(1+m.t_cu).*pi_cu.*X_c - pi_cu.*m.f_x.*w_u ...
           + 1./m.sigma.*pi_uu.*X_u - pi_uu.*m.f_d.*w_u;
    M_c = 0;
    M_u = (epsilon_uu-1)./epsilon_uu.*rphi_uu + (epsilon_cu-1)./epsilon_cu.*rphi_cu ...
          + m.gamma_u.*(1./epsilon_uu.*rphi_uu+1./epsilon_cu.*rphi_cu);
    tariff_c = m.t_cu./(1+m.t_cu).*pi_cu.*X_c;
    tariff_u = m.t_uc./(1+m.t_uc).*pi_uc.*X_u;
    upsilon_u = m.nu.*X_u;

    X_c_out = wL_c + Pi_c + M_c + tariff_c - upsilon_u;
    X_u_out = wL_u + Pi_u + M_u + tariff_u + upsilon_u;
    
end

% 标准化：国家U
w_num = w_u_out;
w_c_out = w_c_out./w_num;
w_u_out = w_u_out./w_num;
r_cc_out = r_cc_out./w_num;
r_cu_out = r_cu_out./w_num;
r_uu_out = r_uu_out./w_num;
r_uc_out = r_uc_out./w_num;
X_c_out = X_c_out./w_num;
X_u_out = X_u_out./w_num;

%% 结果导出
% 福利：是否考虑隐私负外部性
switch m.Privacy
    case 'nonePrivacy' 
    wel_c = X_c./P_c;
    wel_u = X_u./P_u;
    case 'privacyExternality' 
    wel_c = X_c./P_c./(1+(1./m.zeta.*DShare_c).^2);
    wel_u = X_u./P_u./(1+(1./m.zeta.*DShare_u).^2);
end
wel = [wel_c wel_u]';

% 结果导出
n.wel_c = wel_c;
n.wel_u = wel_u;
n.Pi_c = Pi_c;
n.Pi_u = Pi_u;
n.M_c = M_c;
n.M_u = M_u;
n.w_c = w_c;
n.w_u = w_u;
n.P_c = P_c;
n.P_u = P_u;
n.R_c = R_c;
n.R_u = R_u;
n.X_c = X_c;
n.X_u = X_u;
n.r_cc = r_cc;
n.r_uu = r_uu;
n.r_uc = r_uc;
n.r_cu = r_cu;

% 检查国际收支平衡是否有问题
switch m.MarketStruct
    case 'competition'          % 垄断竞争市场结构
    % 检查国际收支平衡
    n.TB_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc) ...
                 + rphi_cu - rphi_uc + upsilon_u;
    n.TB_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu) ...
                 + rphi_uc - rphi_cu - upsilon_u;
    % Good trade deficit
    n.Upsilon_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc);
    n.Upsilon_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu);
    % Data trade deficit
    n.Data_c = rphi_cu - rphi_uc;
    n.Data_u = rphi_uc - rphi_cu;

    % 检查数据要素
    n.check1 = RPhi_c - rphi_cc - rphi_cu;
    n.check2 = RPhi_u - rphi_uu - rphi_uc;
    % 检查劳动力是否出清
    n.l_cc = l_cc;
    n.l_uu = l_uu;
    n.l_cu = l_cu;
    n.l_uc = l_uc;    
    n.l_cc1 = (rphi_cc./r_cc./m.B_c./(D_c.^m.gamma_c)).^(1./(1-m.gamma_c));
    n.l_uu1 = (rphi_uu./r_uu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
    n.l_cu1 = (rphi_cu./r_cu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
    n.l_uc1 = (rphi_uc./r_uc./m.B_c./(D_c.^m.gamma_c)).^(1./(1-m.gamma_c));   
   
    case 'oligopoly'           % 寡头竞争市场结构
    % 检查国际收支平衡
    n.TB_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc) ...
                 + rphi_cu - rphi_uc + upsilon_u;
    n.TB_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu) ...
                 + rphi_uc - rphi_cu - upsilon_u;
    % Good trade deficit
    n.Upsilon_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc);
    n.Upsilon_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu);
    % Data trade deficit
    n.Data_c = rphi_cu - rphi_uc;
    n.Data_u = rphi_uc - rphi_cu;

    % 检查数据要素
    n.check1 = RPhi_c - rphi_cc - rphi_cu;
    n.check2 = RPhi_u - rphi_uu - rphi_uc;
    % 检查劳动力是否出清    
    n.l_cc = l_cc;
    n.l_uu = l_uu;
    n.l_cu = l_cu;
    n.l_uc = l_uc;    
    n.l_cc1 = (rphi_cc./r_cc./m.B_c./(D_c.^m.gamma_c)).^(1./(1-m.gamma_c));
    n.l_uu1 = (rphi_uu./r_uu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
    n.l_cu1 = (rphi_cu./r_cu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
    n.l_uc1 = (rphi_uc./r_uc./m.B_c./(D_c.^m.gamma_c)).^(1./(1-m.gamma_c));   

    case 'monopoly'            % 垄断市场结构
    % 检查国际收支平衡
    n.TB_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc) ...
                 + rphi_cu + upsilon_u;
    n.TB_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu) ...
                 - rphi_cu - upsilon_u;
    % Good trade deficit
    n.Upsilon_c = pi_cu.*X_c./(1+m.t_cu) - pi_uc.*X_u./(1+m.t_uc);
    n.Upsilon_u = pi_uc.*X_u./(1+m.t_uc) - pi_cu.*X_c./(1+m.t_cu);
    % Data trade deficit
    n.Data_c = rphi_cu;
    n.Data_u = -rphi_cu;

    % 检查数据要素
    n.check1 = RPhi_c - rphi_cu;
    n.check2 = RPhi_u - rphi_uu;
    % 检查劳动力是否出清    
    n.l_uu = l_uu;
    n.l_cu = l_cu;
    n.l_uu1 = (rphi_uu./r_uu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
    n.l_cu1 = (rphi_cu./r_cu./m.B_u./(D_u.^m.gamma_u)).^(1./(1-m.gamma_u));
end


end