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?


1 Answer 1


You can use a zero equation for the implicit definition and, if possible, solve it, but only if you have a finite summatino number

n = 5;
G[n] = g*(F/F^2 + e^2);
eq = F - Sum[G[i], {i, 1, n}] == 0;
Fs = Solve[eq, F]

enter image description here

  • $\begingroup$ Thanks, but still I get imaginary numbers in some regions! and I am not sure if it's Mathematica's error or my function just turns out to be like this! $\endgroup$
    – ShS
    Dec 20, 2016 at 18:48
  • $\begingroup$ But there is no method for solve command that can solve it like point by point, or Monte Carlo? $\endgroup$
    – ShS
    Dec 20, 2016 at 18:51
  • 1
    $\begingroup$ 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]. $\endgroup$ Dec 21, 2016 at 10:09

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.