Type a command or search the notes. Tab cycles completions. Up and down recall history when blank. Escape dismisses, then leaves the prompt.
MATH
[1251]

Reduction formulae

Recurrences that replace repeated integration by parts.

A reduction formula replaces repeated integration by parts with a single recurrence. The integral ∫xex dx\int x e^x \, dx takes one integration by parts and ∫x5ex dx\int x^5 e^x \, dx takes five of the same step; a reduction formula does the step once, in general, by writing InI_n in terms of In−1I_{n-1} or In−2I_{n-2}. Trigonometric integrands also need a Pythagorean identity. For low trigonometric powers, the direct identities in trigonometric integrals are quicker.

Reduction formulae as recurrences

[DEFINITION]

Reduction formula

A reduction formula expresses an integral InI_n in terms of In−kI_{n-k} for a smaller index, typically by integration by parts.

Here and below, power indices are non-negative integers. The step kk is 11 for xnexx^n e^x. For trigonometric powers it is 22, because sin⁡2x+cos⁡2x=1\sin^2 x+\cos^2 x=1 or sec⁡2x=1+tan⁡2x\sec^2 x=1+\tan^2 x returns a copy of InI_n that is then collected.

Reducing powers of xx times exe^x

[THEOREM]

Powers of xx times exe^x

For an integer n≥0n\ge0, choose compatible representatives InI_n of the antiderivatives of xnexx^n e^x. Then

In=xnex−nIn−1,n≥1.I_n=x^n e^x-n I_{n-1},\qquad n\ge 1.

Each application lowers the power of xx by one. Two antiderivatives can differ by a constant, so use the recurrence to construct one representative and add a single arbitrary constant when the reduction stops.

[PROOF]

Integrate by parts with u=xnu=x^n and dvdx=ex\frac{dv}{dx}=e^x, so du=nxn−1 dxdu=n x^{n-1} \, dx and v=exv=e^x. Then

In=xnex−n∫xn−1ex dx=xnex−nIn−1.I_n=x^n e^x-n\int x^{n-1} e^x \, dx=x^n e^x-n I_{n-1}.
[EXAMPLE]

∫x3ex dx\int x^3 e^x \, dx

The recurrence with n=3n=3, then n=2n=2, then n=1n=1 and I0=∫ex dxI_0=\int e^x \, dx, gives

I3=x3ex−3I2=x3ex−3(x2ex−2I1)=x3ex−3x2ex+6I1=x3ex−3x2ex+6(xex−I0)=x3ex−3x2ex+6xex−6∫ex dx=x3ex−3x2ex+6xex−6ex+C.\begin{aligned} I_3&=x^3 e^x-3 I_2 \\ &=x^3 e^x-3\bigl(x^2 e^x-2 I_1\bigr) \\ &=x^3 e^x-3x^2 e^x+6 I_1 \\ &=x^3 e^x-3x^2 e^x+6\bigl(x e^x-I_0\bigr) \\ &=x^3 e^x-3x^2 e^x+6x e^x-6\int e^x \, dx \\ &=x^3 e^x-3x^2 e^x+6x e^x-6e^x+C. \end{aligned}

Reducing sine and cosine powers

On [0,π2]\bigl[0,\frac{\pi}{2}\bigr] the boundary term from parts vanishes, so parts and a Pythagorean identity give a multiple of the integral two powers lower.

[THEOREM]

Sine and cosine powers on [0,π2]\bigl[0,\frac{\pi}{2}\bigr]

Let In=∫0π/2sin⁡nx dxI_n=\int_0^{\pi/2}\sin^n x \, dx and Jn=∫0π/2cos⁡nx dxJ_n=\int_0^{\pi/2}\cos^n x \, dx. Then

In=n−1nIn−2,Jn=n−1nJn−2,I_n=\frac{n-1}{n} I_{n-2},\qquad J_n=\frac{n-1}{n} J_{n-2},

for n≥2n\ge 2.

The recurrence stops at I0=π2I_0=\frac{\pi}{2} when nn is even and at I1=1I_1=1 when nn is odd and likewise for JnJ_n. In the proof, collect the reintroduced copy of InI_n before dividing by nn.

[PROOF]

For InI_n, integrate by parts with u=sin⁡n−1xu=\sin^{n-1}x and dvdx=sin⁡x\frac{dv}{dx}=\sin x, so du=(n−1)sin⁡n−2xcos⁡x dxdu=(n-1)\sin^{n-2}x\cos x \, dx and v=−cos⁡xv=-\cos x. Then

In=[−sin⁡n−1xcos⁡x]0π/2+(n−1)∫0π/2sin⁡n−2xcos⁡2x dx.I_n=\bigl[-\sin^{n-1}x\cos x\bigr]_0^{\pi/2}+(n-1)\int_0^{\pi/2}\sin^{n-2}x\cos^2 x \, dx.

The boundary term is 00, since cos⁡π2=0\cos\frac{\pi}{2}=0 and sin⁡0=0\sin 0=0. Replace cos⁡2x=1−sin⁡2x\cos^2 x=1-\sin^2 x:

In=(n−1)∫0π/2sin⁡n−2x(1−sin⁡2x) dx=(n−1)(In−2−In).\begin{aligned} I_n&=(n-1)\int_0^{\pi/2}\sin^{n-2}x\bigl(1-\sin^2 x\bigr) \, dx \\ &=(n-1)\bigl(I_{n-2}-I_n\bigr). \end{aligned}

Hence In=(n−1)In−2−(n−1)InI_n=(n-1)I_{n-2}-(n-1)I_n, so nIn=(n−1)In−2n I_n=(n-1)I_{n-2} and therefore

In=n−1nIn−2,n≥2.I_n=\frac{n-1}{n} I_{n-2},\qquad n\ge 2.

For JnJ_n, take u=cos⁡n−1xu=\cos^{n-1}x and dvdx=cos⁡x\frac{dv}{dx}=\cos x, so du=−(n−1)cos⁡n−2xsin⁡x dxdu=-(n-1)\cos^{n-2}x\sin x \, dx and v=sin⁡xv=\sin x. Then

Jn=[cos⁡n−1xsin⁡x]0π/2+(n−1)∫0π/2cos⁡n−2xsin⁡2x dx.J_n=\bigl[\cos^{n-1}x\sin x\bigr]_0^{\pi/2}+(n-1)\int_0^{\pi/2}\cos^{n-2}x\sin^2 x \, dx.

The boundary term is 00. Replace sin⁡2x=1−cos⁡2x\sin^2 x=1-\cos^2 x:

Jn=(n−1)(Jn−2−Jn).J_n=(n-1)\bigl(J_{n-2}-J_n\bigr).

The same rearrangement yields Jn=n−1nJn−2J_n=\frac{n-1}{n} J_{n-2} for n≥2n\ge 2.

[EXAMPLE]

∫0π/2sin⁡5x dx\int_0^{\pi/2}\sin^5 x \, dx

I5=45I3=45⋅23I1=815∫0π/2sin⁡x dx=815[−cos⁡x]0π/2=815.I_5=\frac{4}{5} I_3=\frac{4}{5}\cdot\frac{2}{3} I_1=\frac{8}{15}\int_0^{\pi/2}\sin x \, dx=\frac{8}{15}\bigl[-\cos x\bigr]_0^{\pi/2}=\frac{8}{15}.
[EXAMPLE]

∫0π/2cos⁡7x dx\int_0^{\pi/2}\cos^7 x \, dx

J7=67J5=67⋅45J3=2435⋅23J1=1635∫0π/2cos⁡x dx=1635[sin⁡x]0π/2=1635.J_7=\frac{6}{7} J_5=\frac{6}{7}\cdot\frac{4}{5} J_3=\frac{24}{35}\cdot\frac{2}{3} J_1=\frac{16}{35}\int_0^{\pi/2}\cos x \, dx=\frac{16}{35}\bigl[\sin x\bigr]_0^{\pi/2}=\frac{16}{35}.

Reducing secant and tangent powers

On [0,π4]\bigl[0,\frac{\pi}{4}\bigr] the boundary term from parts does not vanish: tan⁡π4=1\tan\frac{\pi}{4}=1 and sec⁡π4=2\sec\frac{\pi}{4}=\sqrt{2} leave a power of 22 in the recurrence for sec⁡nx\sec^n x. Tangent powers need no parts at all.

[THEOREM]

Secant powers on [0,π4]\bigl[0,\frac{\pi}{4}\bigr]

Let In=∫0π/4sec⁡nx dxI_n=\int_0^{\pi/4}\sec^n x \, dx. Then

In=1n−1(2(n−2)/2+(n−2)In−2)I_n=\frac{1}{n-1}\bigl(2^{(n-2)/2}+(n-2)I_{n-2}\bigr)

for n≥2n\ge 2.

The base values are I0=∫0π/41 dx=π4I_0=\int_0^{\pi/4}1\,dx=\frac{\pi}{4} and I1=∫0π/4sec⁡x dx=ln⁡(1+2)I_1=\int_0^{\pi/4}\sec x\,dx=\ln(1+\sqrt{2}).

Without the limits, the same splitting sec⁡nx=sec⁡n−2xsec⁡2x\sec^n x=\sec^{n-2}x\sec^2 x produces the antiderivative of sec⁡3x\sec^3 x.

[PROOF]

Integrate by parts with u=sec⁡n−2xu=\sec^{n-2}x and dvdx=sec⁡2x\frac{dv}{dx}=\sec^2 x, so

dudx=(n−2)sec⁡n−2xtan⁡x,v=tan⁡x.\frac{du}{dx}=(n-2)\sec^{n-2}x\tan x,\qquad v=\tan x.

Then

In=[sec⁡n−2xtan⁡x]0π/4−(n−2)∫0π/4sec⁡n−2xtan⁡2x dx.I_n=\bigl[\sec^{n-2}x\tan x\bigr]_0^{\pi/4}-(n-2)\int_0^{\pi/4}\sec^{n-2}x\tan^2 x \, dx.

The boundary term is 2(n−2)/22^{(n-2)/2}, since sec⁡π4=2\sec\frac{\pi}{4}=\sqrt{2} and tan⁡π4=1\tan\frac{\pi}{4}=1, while the lower limit is 00. Replace tan⁡2x=sec⁡2x−1\tan^2 x=\sec^2 x-1:

In=2(n−2)/2−(n−2)∫0π/4sec⁡n−2x(sec⁡2x−1) dx=2(n−2)/2−(n−2)(In−In−2)=2(n−2)/2−(n−2)In+(n−2)In−2.\begin{aligned} I_n&=2^{(n-2)/2}-(n-2)\int_0^{\pi/4}\sec^{n-2}x\bigl(\sec^2 x-1\bigr) \, dx \\ &=2^{(n-2)/2}-(n-2)\bigl(I_n-I_{n-2}\bigr) \\ &=2^{(n-2)/2}-(n-2)I_n+(n-2)I_{n-2}. \end{aligned}

Hence (n−1)In=2(n−2)/2+(n−2)In−2(n-1)I_n=2^{(n-2)/2}+(n-2)I_{n-2} and therefore

In=1n−1(2(n−2)/2+(n−2)In−2),n≥2.I_n=\frac{1}{n-1}\bigl(2^{(n-2)/2}+(n-2)I_{n-2}\bigr),\qquad n\ge 2.
[THEOREM]

Tangent powers on [0,π4]\bigl[0,\frac{\pi}{4}\bigr]

Let Jn=∫0π/4tan⁡nx dxJ_n=\int_0^{\pi/4}\tan^n x \, dx. Then

Jn=1n−1−Jn−2,n≥2.J_n=\frac{1}{n-1}-J_{n-2},\qquad n\ge 2.

This recurrence needs no parts: the identity tan⁡2x=sec⁡2x−1\tan^2 x=\sec^2 x-1 turns the leftover into an elementary integral in u=tan⁡xu=\tan x.

Its base values are J0=π4J_0=\frac{\pi}{4} and J1=∫0π/4tan⁡x dx=ln⁡2J_1=\int_0^{\pi/4}\tan x\,dx=\ln\sqrt{2}.

[PROOF]

Write tan⁡nx=tan⁡n−2xtan⁡2x\tan^n x=\tan^{n-2}x\tan^2 x and replace tan⁡2x=sec⁡2x−1\tan^2 x=\sec^2 x-1:

Jn=∫0π/4tan⁡n−2x(sec⁡2x−1) dx=∫0π/4tan⁡n−2xsec⁡2x dx−Jn−2.J_n=\int_0^{\pi/4}\tan^{n-2}x\bigl(\sec^2 x-1\bigr) \, dx=\int_0^{\pi/4}\tan^{n-2}x\sec^2 x \, dx-J_{n-2}.

Let u=tan⁡xu=\tan x, so du=sec⁡2x dxdu=\sec^2 x \, dx. The limits change from x=0x=0 to x=π4x=\frac{\pi}{4} into u=0u=0 to u=1u=1. Then

Jn=∫01un−2 du−Jn−2=[un−1n−1]01−Jn−2=1n−1−Jn−2J_n=\int_0^1 u^{n-2} \, du-J_{n-2}=\biggl[\frac{u^{n-1}}{n-1}\biggr]_0^1-J_{n-2}=\frac{1}{n-1}-J_{n-2}

for n≥2n\ge 2.

The mixed integrand tan⁡nθsec⁡θ\tan^n\theta\sec\theta on the same interval also drops by two, now with a factor 2\sqrt{2} from the upper limit. Its integrals are written KnK_n to keep them apart from the secant powers InI_n above.

[THEOREM]

Tangent powers times secant

Let Kn=∫0π/4tan⁡nθsec⁡θ dθK_n=\int_0^{\pi/4}\tan^n\theta\sec\theta \, d\theta. Then

Kn=1n(2−(n−1)Kn−2)K_n=\frac{1}{n}\bigl(\sqrt{2}-(n-1)K_{n-2}\bigr)

for n≥2n\ge 2.

Its base values are K0=∫0π/4sec⁡θ dθ=ln⁡(1+2)K_0=\int_0^{\pi/4}\sec\theta\,d\theta=\ln(1+\sqrt{2}) and K1=∫0π/4tan⁡θsec⁡θ dθ=2−1K_1=\int_0^{\pi/4}\tan\theta\sec\theta\,d\theta=\sqrt{2}-1.

[PROOF]

Integrate by parts with u=tan⁡n−1θu=\tan^{n-1}\theta and dvdθ=tan⁡θsec⁡θ\frac{dv}{d\theta}=\tan\theta\sec\theta, so du=(n−1)tan⁡n−2θsec⁡2θ dθdu=(n-1)\tan^{n-2}\theta\sec^2\theta \, d\theta and v=sec⁡θv=\sec\theta. Then

Kn=[tan⁡n−1θsec⁡θ]0π/4−(n−1)∫0π/4tan⁡n−2θsec⁡3θ dθ.K_n=\bigl[\tan^{n-1}\theta\sec\theta\bigr]_0^{\pi/4}-(n-1)\int_0^{\pi/4}\tan^{n-2}\theta\sec^3\theta \, d\theta.

The boundary term is 2\sqrt{2}. Replace sec⁡2θ=1+tan⁡2θ\sec^2\theta=1+\tan^2\theta in the remaining integral:

Kn=2−(n−1)∫0π/4tan⁡n−2θ(1+tan⁡2θ)sec⁡θ dθ=2−(n−1)(Kn−2+Kn).\begin{aligned} K_n&=\sqrt{2}-(n-1)\int_0^{\pi/4}\tan^{n-2}\theta\bigl(1+\tan^2\theta\bigr)\sec\theta \, d\theta \\ &=\sqrt{2}-(n-1)\bigl(K_{n-2}+K_n\bigr). \end{aligned}

Collecting the KnK_n terms gives nKn=2−(n−1)Kn−2nK_n=\sqrt{2}-(n-1)K_{n-2}, so

Kn=1n(2−(n−1)Kn−2),n≥2.K_n=\frac{1}{n}\bigl(\sqrt{2}-(n-1)K_{n-2}\bigr),\qquad n\ge 2.
[EXAMPLE]

Secant and tangent recurrences on [0,π4]\bigl[0,\frac{\pi}{4}\bigr]

The base value I2=∫0π/4sec⁡2x dx=[tan⁡x]0π/4=1I_2=\int_0^{\pi/4}\sec^2 x \, dx=\bigl[\tan x\bigr]_0^{\pi/4}=1 gives

∫0π/4sec⁡4x dx=I4=13(2(4−2)/2+2I2)=13(2+2)=43.\int_0^{\pi/4}\sec^4 x \, dx=I_4=\frac{1}{3}\bigl(2^{(4-2)/2}+2I_2\bigr)=\frac{1}{3}(2+2)=\frac{4}{3}.

The base value J2=∫0π/4tan⁡2x dx=∫0π/4(sec⁡2x−1) dx=[tan⁡x−x]0π/4=1−π4J_2=\int_0^{\pi/4}\tan^2 x \, dx=\int_0^{\pi/4}(\sec^2 x-1)\, dx=\bigl[\tan x-x\bigr]_0^{\pi/4}=1-\frac{\pi}{4} gives

∫0π/4tan⁡4x dx=J4=13−J2=13−(1−π4)=π4−23.\int_0^{\pi/4}\tan^4 x \, dx=J_4=\frac{1}{3}-J_2=\frac{1}{3}-\bigl(1-\frac{\pi}{4}\bigr)=\frac{\pi}{4}-\frac{2}{3}.

For mixed powers, K1=∫0π/4tan⁡θsec⁡θ dθ=[sec⁡θ]0π/4=2−1K_1=\int_0^{\pi/4}\tan\theta\sec\theta \, d\theta=\bigl[\sec\theta\bigr]_0^{\pi/4}=\sqrt{2}-1, so

∫0π/4tan⁡3θsec⁡θ dθ=K3=13(2−2K1)=13(2−2(2−1))=2−23.\int_0^{\pi/4}\tan^3\theta\sec\theta \, d\theta=K_3=\frac{1}{3}\bigl(\sqrt{2}-2K_1\bigr)=\frac{1}{3}\bigl(\sqrt{2}-2(\sqrt{2}-1)\bigr)=\frac{2-\sqrt{2}}{3}.

Reducing two powers together

When both exponents vary, each reduction drops one index by two, as in the one-parameter trigonometric recurrences. Reduce whichever index leads to a known base integral.

[DEFINITION]

Reduction formula

A reduction formula expresses an integral InI_n in terms of In−kI_{n-k} for a smaller index, typically by integration by parts.

[PROPOSITION]

A two-parameter recurrence

Let Im,n=∫0π/2cos⁡mxsin⁡nx dxI_{m,n}=\int_0^{\pi/2}\cos^m x\sin^n x \, dx. Then

Im,n=m−1m+nIm−2,nI_{m,n}=\frac{m-1}{m+n} I_{m-2,n}

if m≥2m\ge 2 and

Im,n=n−1m+nIm,n−2I_{m,n}=\frac{n-1}{m+n} I_{m,n-2}

if n≥2n\ge 2.

The case m=0m=0 recovers the sine-power recurrence; the case n=0n=0 recovers the cosine-power recurrence. The base values needed to stop the reduction are I0,0=π2I_{0,0}=\frac{\pi}{2}, I0,1=1I_{0,1}=1, I1,0=1I_{1,0}=1 and I1,1=12I_{1,1}=\frac{1}{2}.

[PROOF]

Assume m≥2m\ge 2. Integrate by parts with u=cos⁡m−1xu=\cos^{m-1}x and dv=cos⁡xsin⁡nx dxdv=\cos x\sin^n x \, dx, so du=−(m−1)cos⁡m−2xsin⁡x dxdu=-(m-1)\cos^{m-2}x\sin x \, dx and v=sin⁡n+1xn+1v=\frac{\sin^{n+1}x}{n+1}. Then

Im,n=[cos⁡m−1xsin⁡n+1xn+1]0π/2+m−1n+1∫0π/2cos⁡m−2xsin⁡n+2x dx.I_{m,n}=\biggl[\frac{\cos^{m-1}x\sin^{n+1}x}{n+1}\biggr]_0^{\pi/2}+\frac{m-1}{n+1}\int_0^{\pi/2}\cos^{m-2}x\sin^{n+2}x \, dx.

The boundary term is 00. Replace sin⁡2x=1−cos⁡2x\sin^2 x=1-\cos^2 x:

Im,n=m−1n+1∫0π/2cos⁡m−2xsin⁡nx(1−cos⁡2x) dx=m−1n+1(Im−2,n−Im,n).\begin{aligned} I_{m,n}&=\frac{m-1}{n+1}\int_0^{\pi/2}\cos^{m-2}x\sin^n x\bigl(1-\cos^2 x\bigr) \, dx \\ &=\frac{m-1}{n+1}\bigl(I_{m-2,n}-I_{m,n}\bigr). \end{aligned}

Hence m+nn+1Im,n=m−1n+1Im−2,n\frac{m+n}{n+1} I_{m,n}=\frac{m-1}{n+1} I_{m-2,n} and therefore

Im,n=m−1m+nIm−2,n.I_{m,n}=\frac{m-1}{m+n} I_{m-2,n}.

For n≥2n\ge 2 reverse the roles. Take u=sin⁡n−1xu=\sin^{n-1}x and dv=sin⁡xcos⁡mx dxdv=\sin x\cos^m x\,dx, so v=−cos⁡m+1xm+1v=-\frac{\cos^{m+1}x}{m+1}. The boundary term vanishes and cos⁡2x=1−sin⁡2x\cos^2x=1-\sin^2x gives

Im,n=n−1m+1(Im,n−2−Im,n).I_{m,n}=\frac{n-1}{m+1}\bigl(I_{m,n-2}-I_{m,n}\bigr).

Rearranging this equation yields Im,n=n−1m+nIm,n−2I_{m,n}=\frac{n-1}{m+n}I_{m,n-2}.

[EXAMPLE]

∫0π/2cos⁡5xsin⁡4x dx\int_0^{\pi/2}\cos^5 x\sin^4 x \, dx

This is the definite version of the odd-power evaluation of ∫cos⁡5xsin⁡4x dx\int\cos^5 x\sin^4 x \, dx. Reduce the cosine power twice, then the sine power twice and use I1,0=1I_{1,0}=1:

I5,4=49I3,4=49⋅27I1,4=863I1,4=863⋅35I1,2=8105I1,2=8105⋅13I1,0=8315.\begin{aligned} I_{5,4}&=\frac{4}{9} I_{3,4}=\frac{4}{9}\cdot\frac{2}{7} I_{1,4}=\frac{8}{63} I_{1,4} \\ &=\frac{8}{63}\cdot\frac{3}{5} I_{1,2}=\frac{8}{105} I_{1,2} \\ &=\frac{8}{105}\cdot\frac{1}{3} I_{1,0}=\frac{8}{315}. \end{aligned}

On this page