Common Integrals

A short reference table of antiderivatives you will use constantly, plus the rules that let you combine them. These standard forms show up throughout biology: ektdt\int e^{-kt}\,dt gives a drug’s total exposure or the cumulative decay of a labeled tracer, the reciprocal integral 1xdx\int \tfrac1x\,dx produces the ln\ln behind log-scaled assays, and the Gaussian integral at the end normalizes the normal density of biological measurements.

The table

Each row gives an indefinite integral (add a constant CC to any of them).

xndx=xn+1n+1+C(n1)1xdx=lnx+Cexdx=ex+Caxdx=axlna+C(a>0, a1)sinxdx=cosx+Ccosxdx=sinx+C\begin{aligned} \int x^n\,dx &= \frac{x^{n+1}}{n+1} + C \quad (n \neq -1) \\[4pt] \int \frac{1}{x}\,dx &= \ln|x| + C \\[4pt] \int e^x\,dx &= e^x + C \\[4pt] \int a^x\,dx &= \frac{a^x}{\ln a} + C \quad (a > 0,\ a \neq 1) \\[4pt] \int \sin x\,dx &= -\cos x + C \\[4pt] \int \cos x\,dx &= \sin x + C \end{aligned}

The 1/x1/x case is the special exception to the power rule: it is the antiderivative that the formula xn+1/(n+1)x^{n+1}/(n+1) cannot produce (division by zero when n=1n=-1).

Linearity rules

Integration is linear, which lets you break big integrals into small ones:

The Gaussian integral

The area under the “bell curve” kernel ex2e^{-x^2} over the whole real line is

ex2dx=π.\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi} .

This has no elementary antiderivative, yet the total is finite and exact. It is why the normal density 12πσexp ⁣((xμ)22σ2)\frac{1}{\sqrt{2\pi}\,\sigma}\exp\!\big(-\frac{(x-\mu)^2}{2\sigma^2}\big) integrates to 11 — the 2π\sqrt{2\pi} is precisely the normalizing constant.

Worked example

Integrate the polynomial 02(3x2+2x+1)dx\int_0^2 (3x^2 + 2x + 1)\,dx. Using the sum and constant-multiple rules with the power rule:

(3x2+2x+1)dx=x3+x2+x+C.\int (3x^2 + 2x + 1)\,dx = x^3 + x^2 + x + C .

Evaluating from 00 to 22:

[x3+x2+x]02=(8+4+2)0=14.\big[x^3 + x^2 + x\big]_0^2 = (8 + 4 + 2) - 0 = 14 .

Computing it

R

# Numeric check of the polynomial example
integrate(function(x) 3*x^2 + 2*x + 1, 0, 2)$value   # 14
$
# Gaussian integral
integrate(function(x) exp(-x^2), -Inf, Inf)$value    # 1.772454 = sqrt(pi)
$sqrt(pi)                                              # 1.772454

Python

import sympy as sp
from scipy.integrate import quad
import numpy as np

x = sp.symbols("x")
print(sp.integrate(3*x**2 + 2*x + 1, (x, 0, 2)))     # 14
print(sp.integrate(sp.exp(-x**2), (x, -sp.oo, sp.oo)))  # sqrt(pi)

val, _ = quad(lambda t: np.exp(-t**2), -np.inf, np.inf)
print(val, np.sqrt(np.pi))                           # 1.7724538509... 1.7724538509...
14
sqrt(pi)
1.7724538509055159 1.7724538509055159

Julia

using Symbolics, QuadGK

@variables x
D = Differential(x)  # (not needed here, shown for context)

# Numeric checks
println(quadgk(t -> 3t^2 + 2t + 1, 0, 2)[1])         # 14.0
println(quadgk(t -> exp(-t^2), -Inf, Inf)[1])        # 1.7724538509055159
println(sqrt(pi))                                    # 1.7724538509055159

Why it matters for statistics

Densities are built from these functions: the exponential density uses eλxdx\int e^{-\lambda x}dx, the normal uses the Gaussian integral, and moments (means, variances) reduce to power-rule integrals against a density. Knowing the table by sight makes normalizing constants and expected values fast to derive.