CC 9
CC 9
Simpson’s method is a numerical technique used to approximate the definite integral of a function,
particularly when the function is complex or when an analytical solution is difficult or impossible to obtain.
Named after Thomas Simpson, this method provides a more accurate result compared to basic numerical
techniques like the Trapezoidal Rule. It is especially effective when the integrand is a smooth and continuous
function.
Concept and Derivation:
Simpson’s method works by approximating the integrand 𝑓(𝑥) with a second-degree polynomial (i.e., a
parabola) on each subinterval. It uses quadratic interpolation to approximate the area under the curve. The
idea is to divide the interval [𝑎, 𝑏] into an even number 𝑛 of subintervals, each of width ℎ = 𝑏 −𝑛 𝑎 , and then apply
the formula:
∫ 𝑓(𝑥) 𝑑𝑥 ≈ ℎ [𝑓(𝑥0 ) + 4
𝑏 𝑛−1 𝑛−2
𝑎
∑ 𝑓 ( 𝑥𝑖 ) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )]
3 𝑖 = 1, 3, 5, … 𝑖 = 2, 4, 6, …
Here, 𝑥0 = 𝑎, 𝑥𝑛 = 𝑏, and the 𝑥𝑖 's are the intermediate points.
Types:
Simpson’s 1/3 Rule: As shown above, it uses a quadratic polynomial over each pair of intervals.
Simpson’s 3/8 Rule: A less common version, which fits a cubic polynomial over three intervals.
Error Analysis:
Simpson’s rule is known for its accuracy. The error term for the 1/3 rule is proportional to the fourth
derivative of the function:
5
( 𝑏 − 𝑎)
𝐸= − 𝑓(4) (𝜉), for some 𝜉 ∈ ( 𝑎, 𝑏 )
180𝑛
4
This means that the error decreases rapidly as 𝑛 increases, especially if the function is smooth.
Applications:
Simpson’s method is widely used in physics, engineering, and computational sciences for:
The Runge-Kutta (RK) methods are a family of iterative techniques used to approximate solutions of ordinary
differential equations (ODEs), particularly initial value problems of the form:
𝑑𝑦
= 𝑓(𝑥, 𝑦), 𝑦(𝑥0 ) = 𝑦
𝑑𝑥 0
Among these, the most commonly used is the fourth-order Runge-Kutta method (RK4) due to its balance
between accuracy and computational efficiency.
RK4 Formula:
Given the current point (𝑥𝑛 , 𝑦𝑛 ) and a step size ℎ, the next value 𝑦𝑛 + 1 is calculated as:
𝑘1 = ℎ𝑓(𝑥𝑛 , 𝑦𝑛 )
ℎ, 𝑦 𝑘1
𝑘2 = ℎ𝑓(𝑥𝑛 + + )
2 𝑛 2
ℎ, 𝑦 𝑘2
𝑘3 = ℎ𝑓(𝑥𝑛 + + )
2 𝑛 2
𝑘4 = ℎ𝑓(𝑥𝑛 + ℎ, 𝑦𝑛 + 𝑘3 )
1
𝑦𝑛 + 1 = 𝑦 𝑛 + (𝑘 1 + 2 𝑘 2 + 2 𝑘 3 + 𝑘 4 )
6
Each 𝑘𝑖 represents an estimate of the slope at different points within the interval.
Advantages:
High accuracy even with a relatively large step size.
No need for higher derivatives as in Taylor series methods.
Self-starting and does not require previous values as in multi-step methods.
Error:
RK4 has a local truncation error of order 𝑂(ℎ ), and a global error of 𝑂(ℎ ), making it much more accurate
5 4
Limitations:
Computationally intensive compared to simpler methods.
Not suitable for stiff equations without modification.
Monte Carlo Integration is a probabilistic method used to estimate definite integrals, especially in high-
dimensional spaces where traditional numerical methods become inefficient or infeasible. It relies on the law
of large numbers, using random sampling to approximate the value of an integral.
Basic Idea:
To evaluate an integral:
𝑏
𝐼= ∫ 𝑓(𝑥)𝑑𝑥
𝑎
Monte Carlo method approximates this as:
𝐼 ≈ 𝑏 −𝑛 𝑎 ∑ 𝑓(𝑥𝑖 )
𝑛
𝑖=1
Where 𝑥𝑖 are randomly selected points in the interval [𝑎, 𝑏].
In multiple dimensions:
𝑉(Ω) 𝑛
𝐼= ∫ 𝑓(𝑥)𝑑𝑥 ≈ 𝑛 ∑ 𝑓(𝑥𝑖 )
Ω 𝑖=1
Here, Ω is the domain of integration and 𝑉(Ω) is its volume.
Key Features:
Particularly useful for multidimensional integrals.
Works well for irregular or discontinuous functions.
Easily parallelizable, which suits high-performance computing environments.
Advantages:
Dimension-independent convergence.
Flexibility in handling complex domains.
Simple to implement.
Limitations:
Slow convergence rate.
High variance in results unless variance-reduction techniques (e.g., importance sampling, stratified
sampling) are used.
May require a large number of samples for accurate results.
Let me know if you want these formatted for a report or LaTeX document.