function   f=xi_wss(xx,omega_r,xi_r,sigmaH,betaH,R,g_n,kai)

xi_w=xx(1);

Omega= 1 - omega_r + omega_r*(xi_r/xi_w)^((sigmaH - 1)/ sigmaH);

% 求解的方程
f1=xi_w-1+(betaH*R*Omega)^(1/(1 - sigmaH))/(R*Omega)*g_n^(kai*sigmaH/(1 - sigmaH));

f=[f1];