I couldn’t attach pic of my self consistent equations in the comment section that I am trying to solve in an earlier message threading

(Elimination of the variables)

I am trying to solve for z for given x = 0.057, y = 0.01 **self consistently**

Code from earlier threading is

`y = 1/100; x = 57/1000; gE[(n_Integer)?NonNegative] := 1/g[n] - Sqrt[1 + z^2/(Pi*(2*n + 1)*y + (44*x*Pi*g[n])/10)^2]; sce[0] = (2*Pi*1*((Pi*1*(g[0] - 1))/(100*(((Pi*1)/100 + (44*Pi*x)/10)*((Pi*1)/100 + (44*Pi*x*g[0])/10)))))/100 + 6431/10000; sce[(n_Integer)?Positive] := (2*Pi*1*(((2*n + 1)*Pi*1*(g[n] - 1))/(100*((((2*n + 1)*Pi*1)/100 + (44*Pi*x)/10)*(((2*n + 1)*Pi*1)/100 + (44*Pi*x*g[n])/10)))))/100; nmax = 3; nValues = Range[0, nmax]; ee = Flatten[{Total[sce /@ nValues], gE /@ nValues}]; xed0 = GroebnerBasis[ee, {z}, g /@ nValues, MonomialOrder -> EliminationOrder]; NSolve[xed0 == 0, z] `

For nmax = 3 it doesn’t converge and keeps on running and I am aware of @ Daniel Lichtblau comment that here complexity arises as nmax increases.

Question: Is there any other way of solving above self consistent equations?