0
$\begingroup$

I have some problems in solving this equation numerically in Mathematica. Mathematica gives me the error message "NSolve::nsmet: This system cannot be solved with the methods available to NSolve". This is the code:

f[x_, d_] := NIntegrate[z/(Exp[z] - 1), {z, x, x Sqrt[1 + d^2]}]

NSolve[f[x, d] == e x^2 d, d]

where e and x are two parameters. Can anyone help me?

$\endgroup$
2
  • $\begingroup$ You can't use NIntegrate when there are unknowns in the integrand or in the limits as in your case. It is a numerical integration after all. You can use Integrate. But the resulting equation can't be solved for d, too complicated. If you have some numerical values for x and e, there could be a better chance. $\endgroup$
    – Nasser
    Oct 17, 2021 at 8:46
  • $\begingroup$ I simplified the code in this: f[d_] := Integrate[z/(Exp[z] - 1), {z, 1, Sqrt[1 + d^2]}] NSolve[f[r] == r, r] But give me the same error $\endgroup$
    – Tony Stack
    Oct 17, 2021 at 8:57

2 Answers 2

0
$\begingroup$

Lets solve the integral with Integrate and restrict to x>0,d>0. A plot shows at which d the one sided equation f[x, d] - e x^2 d is zero in dependence of x (and e, later set to 3).

f[x_, d_] = 
Integrate[z/(Exp[z] - 1), {z, x, x Sqrt[1 + d^2]}, 
   Assumptions -> {x > 0, d > 0}]

(*   -(1/2) x (-2 I (-1 + Sqrt[1 + d^2]) \[Pi] + d^2 x + 2   Log[-1 + E^x] - 
2 Sqrt[1 + d^2] Log[-1 + E^(Sqrt[1 + d^2] x)]) - PolyLog[2, E^x] +
PolyLog[2, E^(Sqrt[1 + d^2] x)]   *)

Manipulate[
Plot[f[x, d] - e x^2 d // Chop, {d, 0, 30}, 
PlotRange -> 1], {{x, .15}, 0, 1, Appearance -> "Labeled"}, {{e, 3},
0, 5, Appearance -> "Labeled"}]

In order to get the two d values where equation is zero (here called dzero1 and dzero2), regard d as a function of x, d[x], differentiate for x ( you get also d'[x]) and solve with NDSolve with two initial conditions for d (her for d[.15]) and plot result (attention: LogPlot)

dd = D[0 == f[x, d] - e x^2 d /. {e -> 3} /. d -> d[x], x] // 
     ExpandAll // Simplify

d01 = d /. 
FindRoot[f[x, d] - e x^2 d /. {e -> 3} /. x -> .15, {d, 3}] // Chop

d02 = d /. 
FindRoot[f[x, d] - e x^2 d /. {e -> 3} /. x -> .15, {d, 15}] // Chop

dzero1 = d /. First@NDSolve[{dd, d[15/100] == d01}, d, {x, 10^-6, 1}]

dzero2 = d /. First@NDSolve[{dd, d[15/100] == d02}, d, {x, 10^-6, 1}]

LogPlot[Evaluate[{dzero1[x], dzero2[x]}], {x, 10^-6, .21}, 
  PlotRange -> 100]

enter image description here

$\endgroup$
1
$\begingroup$

With your updated/simplified problem, Reduce and Solve can do it

Clear["Global`*"]
f[d_] := Integrate[z/(Exp[z] - 1), {z, 1, Sqrt[1 + d^2]}, GenerateConditions -> False];
eq = f[r] == r

Mathematica graphics

Reduce[eq, r, Reals]

Mathematica graphics

Verify solution

 Simplify[eq /. r -> 0]

Mathematica graphics

And Solve

Solve[eq, r, Reals]

Mathematica graphics

And NSolve

NSolve[eq, r, Reals]

Mathematica graphics

But need to use the Reals options in all of them.

Update

Fixing only "e" and making "x" a free variable, can't be solved?

A quick experiment shows that solution is always d=0?

enter image description here

It could not solve it without having specific values for x and e. At least I could not find a way. Tried Solve and Reduce. May be someone else can find a better method to tackle this.

Clear["Global`*"];
f = Integrate[z/(Exp[z] - 1), {z, x, x Sqrt[1 + d^2]}, 
   GenerateConditions -> False];
Manipulate[
 Module[{eq},
  eq = f == e *x^2* d;
  Print[eq];
  eq = eq /. {x -> x0, e -> e0};
  Print[eq];
  Grid[{{Row[{"Equation =", eq}]}, {Row[{"solution = ", 
       Reduce[eq, d, Reals]}]}}]
  ],
 {{e0, 1, "e"}, 0, 10, 1, Appearance -> "Labeled"},
 {{x0, 1, "x"}, 0, 10, 1, Appearance -> "Labeled"},
 TrackedSymbols :> {e0, x0},
 SynchronousUpdating -> False 
 ]
$\endgroup$
1
  • $\begingroup$ Thanks! Fixing only "e" and making "x" a free variable, can't be solved? $\endgroup$
    – Tony Stack
    Oct 17, 2021 at 9:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

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