Die Integrale ∫∞0dx(a+bxc)d(e+fxg)h∫0∞dx(a+bxc)d(e+fxg)h\int_0^\infty \frac{\mathrm{d}x}{( a+bx^c)^d(e+fx^g)^h} und ∫∞0xe(a+bxc)ddx∫0∞xe(a+bxc)ddx\int_0^\infty \frac{x^e} {(a+bx^c)^d}\mathrm{d}x

Dank des Kommentars von Zacky habe ich die Antwort auf das zweite Integral gefunden, indem ich die Beta-Funktion verwendet habe, insbesondere die Integraldefinition

$$B(x,y)=\int_0^\infty \frac{t^{x-1}}{(1+t)^{x+y}}\mathrm{d}t$$ Lösung geben $$\large{\int_0^\infty \frac{x^e}{(a+bx^c)^d}\mathrm{d}x=\frac{a^{\frac{e+1}{c}}\Gamma(\frac{e+1}{c})\Gamma(d-\frac{e+1}{c})}{a^db^{\frac{e+1}{c}}c\,\Gamma(d)}}$$

Konzentrieren wir uns auf das erste Integral. Wir nehmen an, dass$a$,$b$,$e$Und$f$sind positive reelle Zahlen und das$c$,$g$,$d$Und$h$sind alle positiven ganzen Zahlen größer gleich eins. Ausgangspunkt ist die Partialbruchzerlegung:

$$\begin{eqnarray} \frac{1}{a+b x^c} = \frac{1}{b c} \cdot \left( \frac{b}{a}\right)^{\frac{c-1}{c}} e^{-\imath \frac{\pi}{c} (c-1)} \sum\limits_{p=0}^{c-1} \frac{(-1)^{\frac{2 p}{c}}}{x - (\frac{a}{b})^{1/c} e^{\imath \pi \frac{2 p+1}{c}}} \end{eqnarray}$$ Schreiben Sie nun eine weitere Kopie des obigen mit$(a,b,c)$ersetzt werden durch$(e,f,g)$jeweils und dann Multiplizieren des Originals mit der Kopie ergibt: $$\begin{eqnarray} &&\frac{1}{a+b x^c}\frac{1}{e+f x^g} = \frac{e^{-i \pi \left(-\frac{1}{c}-\frac{1}{g}+2\right)} \left(\frac{b}{a}\right)^{1-\frac{1}{c}} \left(\frac{f}{e}\right)^{1-\frac{1}{g}}}{b c f g} \cdot \\ && \sum\limits_{p=0}^{c-1} \sum\limits_{q=0}^{g-1} \left( \frac{1}{x-e^{\frac{i \pi (2 p+1)}{c}} \left(\frac{a}{b}\right)^{\frac{1}{c}}}-\frac{1}{x-e^{\frac{i \pi (2 q+1)}{g}} \left(\frac{e}{f}\right)^{\frac{1}{g}}} \right) \frac{(-1)^{\frac{2 p}{c}+\frac{2 q}{g}}}{e^{\frac{i \pi (2 p+1)}{c}} \left(\frac{a}{b}\right)^{\frac{1}{c}}-e^{\frac{i \pi (2 q+1)}{g}} \left(\frac{e}{f}\right)^{\frac{1}{g}}} \end{eqnarray}$$ Der obige Ausdruck eignet sich perfekt, um die Stammfunktion in Bezug auf zu finden$x$. Nachdem wir dies getan haben, werten wir das Integral bezüglich aus$x$von null bis unendlich und dann berechnen wir die$a$-Derivat$(d-1)$mal und dann die$e$-Derivat$(h-1)$mal. Nachdem wir dies getan haben, multiplizieren wir das Ganze mit$(-1)^{d+h}/((d-1)!(h-1)!)$. Das gibt: $$\begin{eqnarray} &&\int\limits_0^\infty \frac{1}{(a+b x^c)^d} \frac{1}{(e+f x^g)^h} dx=\frac{(-1)^{d+h}}{b c f g (d-1)! (h-1)!} e^{-i \pi \left(-\frac{1}{c}-\frac{1}{g}+2\right)}\cdot \\ && \frac{\partial^{d-1} }{\partial \xi^{d-1}} \frac{\partial^{h-1} }{\partial \eta^{h-1}} \left.\left[ \left(\frac{b}{\xi }\right)^{1-\frac{1}{c}} \left(\frac{f}{\eta }\right)^{1-\frac{1}{g}} \sum\limits_{p=0}^{c-1} \sum\limits_{q=0}^{g-1} \left( -\frac{\log \left(\frac{\xi }{b}\right)}{c}+\frac{i \pi (2 c q+c-2 g p-g)}{c g}+\frac{\log \left(\frac{\eta }{f}\right)}{g} \right) \frac{(-1)^{\frac{2 p}{c}+\frac{2 q}{g}}}{e^{\frac{i \pi (2 p+1)}{c}} \left(\frac{\xi }{b}\right)^{\frac{1}{c}}-e^{\frac{i \pi (2 q+1)}{g}} \left(\frac{\eta }{f}\right)^{\frac{1}{g}}} \right]\right|_{(\xi,\eta)=(a,e)} \end{eqnarray}$$ Es ist klar, dass das Ergebnis noch vereinfacht werden kann, indem man die Kettenregel der Differentiation anwendet.

For[count = 1, count <= 50, count++,
  {a, b, e, f} = RandomReal[{0, 10}, 4, WorkingPrecision -> 50];
  {c, g, d, h} = RandomInteger[{1, 10}, 4];
  I1 = NIntegrate[1/(a + b x^c)^d 1/(e + f x^g)^h, {x, 0, Infinity}, 
    WorkingPrecision -> 30];
  I2 = 1/(b c ) 1/(f g ) (-1)^(
    d + h)/((d - 1)! (h - 1)!) (D[
       Exp[-I Pi (2 - 1/c - 1/g)] (b/xi)^(1 - 1/c) (f/eta)^(1 - 1/g)
         Sum[( -1/c Log[ xi/b ] + 1/g Log[eta/f ] + (
            I \[Pi] (c - g - 2 g p + 2 c q))/(c g)) (-1)^(
          2 p/c + 2 q/g)/((xi/b)^(1/c)
             Exp[(I Pi (2 p + 1))/c] - (eta/f)^(1/g)
             Exp[(I Pi (2 q + 1))/g]), {p, 0, c - 1}, {q, 0, 
          g - 1}], {xi, d - 1}, {eta, h - 1}] /. {xi :> a, eta :> e});
  If[Abs[I2/I1 - 1] > 10^(-3), 
   Print["mismatch..", {a, b, e, f, c, g, d, h}]; Break[]];
  PrintTemporary[count];
  ];