4

Those integral techniques I thought I'd never need to remember are coming back to bite me.

Is there a way to compute integrals of the form $$ \int_0^\infty \frac1{\sqrt{1 + a x^2 + b x^4}} dx $$

I am feeling hopeless about this, but any reference/hints will be appreciated.

2 Answers2

6

Perhaps you still don't need to remember them, because they don't help much with this integral. OK, substitution helps a bit. But, the integral cannot be expressed in terms of functions that you are familiar with. It can be only expressed in terms of the complete elliptic integral of first kind, denoted by $K$ below.

If you remember the quadratic formula, you will not be surprised to learn that the sign of $a^2-4b$ matters. (By the way, I assume that $a,b$ in your post are positive numbers.)

Case 1: $a^2\ge 4b$. The substitution $u=b^{1/4}x$ changes the integral to $$b^{-1/4}\int_0^\infty \frac1{\sqrt{1+cu^2+u^4}}\,du \label{1}\tag{1}$$ where $c=a/\sqrt{b}\ge 2$. We can write $c=k+k^{-1}$ for some $0<k\le 1$ and get $$\int_0^\infty \frac1{\sqrt{1+cu^2+u^4}}\,du = \sqrt{k}\,K\left(\sqrt{1-k^2}\right) \label{2}\tag{2}$$ courtesy of Maple. (I don't see right away which substitution changes the left side of $(\ref{2})$ to the definition of $K$.)

Case 2: $a^2\le 4b$. Substitute $u=a^{1/2}x$ to change the integral to $$a^{-1/2}\int_0^\infty \frac1{\sqrt{1+u^2+cu^4}}\,du \label{3}\tag{3}$$ where $c=b/a^2\ge 1/4$. Here, $$\int_0^\infty \frac1{\sqrt{1+u^2+cu^4}}\,du = c^{-1/4}K\left(\frac12 \sqrt{2-1/\sqrt{c}}\right) \label{4}\tag{4}$$ Again, I fail to justify $(\ref{4})$ with an explicit transformation.

Sanity check: when $a^2=4b$, the first formula yields $b^{-1/4}K(0)$, while the second gives $a^{-1/2}(b/a^2)^{-1/4}K(0)$. Same thing.

  • Thank you for your very thorough answer. It turns out, that for the case I care about, the coefficients were special case such that the K function evaluates to something simple.

    As a side note, I tried doing this in mathematica but it timed out. Is maple better for these sorts of things?

    – user63235 Feb 21 '13 at 05:57
  • @user63235 Maple would not do anything either, if given the integral as stated. It succeeded only after I helped by rewriting the integral and splitting into two cases as above. For example, in the second case the command was int(...) assuming c>1/4. Mathematica would probably perform similarly, but I don't have it. –  Feb 21 '13 at 11:39
1

"I don't see right away which substitution changes the left side of (2) to the definition of $K$."

Consider first the indefinite integral

$$\int\frac{\mathrm du}{\sqrt{au^4+bu^2+1}}$$

Letting

$$h=\frac{b\pm\sqrt{b^2-4 a}}{2}$$

we have the factorization

$$\frac1{\sqrt{a}}\int\frac{\mathrm du}{\sqrt{\left(u^2+\frac{h}{a}\right)\left(u^2+\frac1{h}\right)}}$$

We now proceed to make a trigonometric substitution, $u=\frac{\tan\,v}{\sqrt{h}}$, to yield

$$\frac1{\sqrt{h}}\int\frac{\mathrm dv}{\sqrt{1-\left(1-\frac{a}{h^2}\right)\sin^2 v}}$$

At this point, we now pay attention to the limits; since $\arctan(0)=0$ and $\arctan(\infty)=\pi/2$, we then have

$$\frac1{\sqrt{h}}\int_0^{\pi/2}\frac{\mathrm dv}{\sqrt{1-\left(1-\frac{a}{h^2}\right)\sin^2 v}}=\frac1{\sqrt{h}}K\left(1-\frac{a}{h^2}\right)$$

In general, $h$ will be complex, necessitating the evaluation of the complete elliptic integral of the first kind for complex arguments. Thankfully, if you're using the AGM to evaluate the complete elliptic integral, it will still work; all that's needed is complex arithmetic.