Put z3n+zˉ3n=fn(z). As fn(z)=fn(z), we get fn(z)∈R, for all z∈C. We should determine the set M={z∈C/fn(z)≥0}. If z∈M, then f0(z)=z+zˉ≥0, so Re(z)≥0.
For non-negative z∈R, fn(z)=2z3n≥0, for all non-negative integers n, implying [0,∞)⊂M.
If z=bi∈C, b∈R, then zˉ=−z, so fn(z)=0, for n∈N, so {bi/b∈R}⊂M. As fn(z)=fn(zˉ), we have z∈M⇔zˉ∈M.
It remains to find z∈M such that Re(z)>0 and Im(z)>0. Let z=r(cost+isint) with r>0, t∈(0,2π). Then fn(z)=2r3ncos(3nt), so z∈M⇔cos(3nt)≥0, for all n∈N.
As t∈(0,2π), there is k∈N such that 3kt≤2π<3k+1t (in fact k=⌊log32tπ⌋). If 3kt=2π, then cos(3nt)>0, for all n∈{0,1,2,…,k−1} and cos(3nt)=0, for all n≥k, that is z∈M. If 3kt<2π<3k+1t, then 2π<3k+1t<23π, implying cos(3k+1t)<0, that is z∈/M.
In conclusion
M=[0,∞)∪{r(cos2⋅3kπ±isin2⋅3kπ)/r>0,k∈N}.