# Solving a coupled equation (non-linear)

I have a function $F[x]= \sum_i G_i[x]$ (in term of some other known functions of $G_i[x]$), except $G_a[x]=g[x](\frac{F[x]}{F[x]^2+\epsilon^2})$ (g[x] is known and well-behave)

Now, of course I cannot define my $F$[$x$_]$:=\sum_i G_i[x]$.

So, I decided to use Solve command,$Func[x\_]:= F[x]/. Solve[F[x]==\sum_i G_i [x],F[x],WorkingPrecision -> MachinePrecision][[1]]$

but it really doesn't look like a right command or at least it needs a method to avoid inconsistent results. Any idea please?

n = 5;

• You just might have bad luck with your function, I mean, imaginary solutions are just solutions. You might restrict the solutions with Solve[eq,F,Reals]. – Mauricio Fernández Dec 21 '16 at 10:09