Romberg integration in SciPy

Introduction

Romberg integration helps us find the area under a curve more accurately by improving simple methods step by step.

When you want a precise estimate of an integral without complicated formulas.

When you have a function that is smooth and you want to calculate its total value over an interval.

When you want to improve the accuracy of basic numerical integration methods like trapezoidal rule.

When you need a quick and reliable way to integrate functions in scientific or engineering problems.

Syntax

SciPy

from scipy.integrate import romberg

result = romberg(function, a, b, tol=1.48e-08, divmax=10, show=False)

function is the function you want to integrate.

a and b are the start and end points of the interval.

Examples

Integrate x² from 0 to 1 using Romberg integration.

SciPy

from scipy.integrate import romberg

def f(x):
    return x**2

result = romberg(f, 0, 1)
print(result)

Calculate the integral of sin(x) from 0 to π, which should be close to 2.

SciPy

from scipy.integrate import romberg

import math

result = romberg(math.sin, 0, math.pi)
print(result)

Sample Program

This program calculates the integral of the function e^(-x²) between 0 and 1 using Romberg integration. This function is common in probability and statistics.

SciPy

from scipy.integrate import romberg
import math

def f(x):
    return math.exp(-x**2)

# Integrate e^(-x^2) from 0 to 1
result = romberg(f, 0, 1)
print(f"Integral of e^(-x^2) from 0 to 1 is approximately {result:.6f}")

Important Notes

Romberg integration uses repeated trapezoidal approximations and improves accuracy by combining them.

It works best for smooth functions without sharp corners or discontinuities.

You can control accuracy with the tol parameter and maximum steps with divmax.

Summary

Romberg integration is a smart way to get accurate integral values by improving simple methods.

It is easy to use with scipy.integrate.romberg by giving a function and interval.

Great for smooth functions and when you want better accuracy than basic methods.