Simpson's Rule Calculator
Approximate definite integrals using Simpson's 1/3 rule and 3/8 rule for f(x) = ax² + bx + c.
Compares to the exact value and shows error reduction with more subdivisions.
Simpson’s rule approximates a definite integral by fitting parabolas through each pair of subintervals rather than rectangles. It is far more accurate than simple Riemann sums for smooth functions.
Simpson’s 1/3 rule (n must be even): integral ≈ (dx/3) * [f(x0) + 4f(x1) + 2f(x2) + 4f(x3) + … + 4f(x_{n-1}) + f(x_n)]
The alternating 4-2-4-2 pattern of coefficients is the signature of this rule. Each group of three points (two subintervals) fits one parabola. Requiring n to be even ensures the groups pair up correctly.
Simpson’s 3/8 rule (n must be divisible by 3): integral ≈ (3dx/8) * [f(x0) + 3f(x1) + 3f(x2) + 2f(x3) + 3*f(x4) + … + f(x_n)]
For a polynomial of degree 3 or less, both Simpson rules give the exact answer regardless of n. For a quadratic f(x) = ax^2 + bx + c, even n=2 (three points) is exact.
The error in Simpson’s 1/3 rule is proportional to dx^4 (proportional to 1/n^4). That is dramatically better than the trapezoidal rule (1/n^2) – halving the step size reduces the error by a factor of 16 rather than 4.
Gauss-Legendre quadrature and adaptive quadrature can do even better for smooth functions, but Simpson’s rule is what most engineering and physics courses use as the go-to numerical integration method.