
function [pm_pn_ss] = SS_LS1_aux(X1, X2, alp, psi)
% Define the anonymous function for fsolve
fun = @(pm_pn_ss) ps2pn_function(pm_pn_ss, X1, X2, alp, psi);

% Initial guess
x0 = 1;

% Solve the system
options = optimoptions('fsolve','Display','none'); % Use 'iter' for iteration display
[solution, ~, ~, ~] = fsolve(fun, x0, options);
pm_pn_ss = solution;
% Display result
% disp('pm_pn_ss:');
% disp(pm_pn_ss);

