The Newton-Cotes integration formulas are the most commonly used numerical integration methods in numerical analysis. The strategy involved in all these Newton cotes formulas is replacing a complicated function with an approximating function which helps us perform the integration easily. Thus, Newton cote integration formulas play an important role in solving numerical integration problems. These formulas help in deriving the simple form of given expressions in integration. Let’s learn different types of Newton cotes integration formulas here in this article.
What are Newton Cotes Integration Formulas?
In numerical analysis, the Newton–Cotes formulas, or Newton Cotes rules, are also called the Newton–Cotes quadrature rules, which comprises a group of formulas for numerical integration based on estimating the integrand at equally spaced points. Also, these formulas or rules were named after Isaac Newton and Roger Cotes. Below are a few important Newton cotes integration rules.
- Trapezoidal Rule
- Simpson’s Rules
- Simpson’s 1/3 Rule
- Simpson’s 3/8 Rule
Before learning about all these rules, let’s recall what is numerical integration.
Numerical integration:
Numerical integration is the approximate estimation of integral employing numerical methods. Sometimes, the numerical estimation of an integral is called quadrature.
Read more: |
Trapezoidal Rule
The trapezoidal rule uses the first-degree polynomial to replace the integrand, i.e., the function to be integrated.
Consider the integration of f(x) over the interval [a, b]
I = ∫ab f(x) = ∫ab f1(x)
Here, f1(x) can be written as:
f1(x) = a + [{f(b) – f(a)}/ (b – a)] (x – a)
Therefore,
I = (b – a) [f(a) + (b)]/2
The only efficient way to enhance the precision of the trapezoidal rule is to divide the integration interval [a, b] into several segments and apply the technique to each of these segments. Then, by adding the areas of all these parts, we can get the integration of the entire function. Mathematically, we can define this technique as follows.
Suppose f(x) is a continuous function on the interval [a, b]. Now divide the interval [a, b] into n equal subintervals with each of width,
Δx = h = (b – a)/n, such that a = x0 < x1< x2< x3<…..<xn = b and a = x0, b = xn.
Then the Trapezoidal Rule formula for area approximating the definite integral ∫ab f(x) is given by:
Here,
xi = a + iΔx = a + ih
If n →∞, R.H.S of the expression approaches the definite integral
This trapezoidal method of converting the interval of integration into segments is called the Multiple trapezoidal rule.
Click here to learn more about the trapezoidal rule of integration.
Simpson’s Rules
The most efficient way to get a better accurate estimation of an integral is to operate a higher-order polynomial to join the points. This way of estimating the integrals is known as Simpson’s Rule.
Simpson’s formula for n+1 equally spaced subdivision is given by;
Here,
n = even number
â–³x = h = (b – a)/n
xi = a + iâ–³x = a + ih
If the function y = f(x) is equally spaced between [a, b] and if a = x0, x1 = x0 + h, x2 = x0 + 2h …., xn = x0 + nh, where h is the difference between the terms, then the analogous values of y with each value of x will be y0 = f(x0), y1 = f(x1), y2 = f(x2),……, yn = f(xn).
Simpson’s 1/3 Rule
Simpson’s one-third rule or Simpson’s â…“ rule is generally considered as an extension of the trapezoidal rule in which a second-order polynomial approximates the integrand. Simpson’s rule can be derived from the various ways using Newton’s divided difference polynomial, Lagrange polynomial, and the coefficients method. Simpson’s 1/3 rule is defined by:
∫ab f(x) dx = h/3 [(y0 + yn) + 4(y1 + y3 + y5 + …. + yn-1) + 2(y2 + y4 + y6 + ….. + yn-2)]
This rule is known as Simpson’s One-third rule.
Simpson’s 3/8 Rule
Simpson’s 3/8 rule is based on cubic rather than quadratic interpolation. Simpson’s 3/8 rule or three-eight rule is given by the formula:
∫ab f(x) dx = 3h/8 [(y0 + yn) + 3(y1 + y2 + y4 + y5 + …. + yn-1) + 2(y3 + y6 + y9 + ….. + yn-3)]
The 3/8 rule is known as Simpson’s second rule of integration.
Get more information about Simpson’s rule here
Apart from the above formulas, we can also define some special formulas such as open and closed newton cotes.
Closed and Open Newton Cotes Formulas
Let us assume that the value of f(x) defined on the interval [a, b] is known at n + 1 equally spaced points, i.e., x0, x1, x2,…, xn such that a ≤ x0 < x1 < x2 < …. < xn ≤ b. Thus, there will be two classes of Newton cotes quadrature which are called closed and open Newton cotes. Closed Newton cotes occur when a = x0 and b = xn, i.e., they use the function values at the endpoints of the interval. Also, open cotes occur when a < x0 and xn < b, which means they do not use the values of the function at the endpoints.
Open and closed Newton cotes formulas using (n +1) points can be defined as follows.
Here,
For an open newton cotes formula, xi = a + iâ–³x = a + ih such that h = (b – a)/n.
For a closed newton cotes formula, xi = a + (i + 1)â–³x = a + (i + 1)h such that h = (b – a)/(n + 2).
Also, h represents the step size and wi represents the weights.
Open Newton Cotes Formulas
Some open Newton cotes formulas for integration are tabulated below.
n |
h = (b – a)/(n + 2) |
Common name |
Formula |
0 |
(b – a)/2 |
Midpoint rule (or) Rectangle rule |
2hf1 |
1 |
(b – a)/3 |
(3/2)h (f1 + f2) |
|
2 |
(b – a)/4 |
Milne’s rule |
(4/3)h (2f1 – f2 + 2f3) |
3 |
(b – a)/5 |
(5/24)h (11f1 + f2 + f3 + 11f4) |
Closed Newton Cotes Formulas
A few closed Newton cotes formulas for integration are given in the table below.
n |
h = (b – a)/n |
Common name |
Formula |
0 |
Undefined |
– |
– |
1 |
b – a |
Trapezoidal rule |
(h/2)(f0 + f1) |
2 |
(b – a)/2 |
Simpson’s rule |
(h/3) (f0 + 4f1 + f2) |
3 |
(b – a)/3 |
Simpson’s ⅜ rule |
(3h/8) (f0 + 3f1 + 3f2 + f3) |
4 |
(b – a)/4 |
Boole’s rule |
(2h/45) (7f0 + 32f1 + 12f2 + 32f3 + 7f4) |
Solved Examples
Example 1:
What is the Newton Cotes formula for Trapezoidal rule?
Solution:
We can write the closed Newton cotes formula for Trapezoidal rule as:
(h/2)(f0 + f1)
Here,
h = (b – a)
Example 2:
Write the Open Newton Cotes integration for n = 2.
Solution:
For n = 2, the Open Newton cotes formula for integration is given by:
(4/3)h (2f1 – f2 + 2f3)
Here,
n = 2
h = (b – a)/4
Frequently Asked Questions
What is Open Newton Cotes formula?
The Open Newton Cotes formula can be defined for a < x0 and xn < b, which means the equally spaced points do not use the values of the function f(x) at the endpoints of the interval [a, b].
What is Simpson’s ⅜ rule formula?
∫ab f(x) dx = h/3 [(y0 + yn) + 4(y1 + y3 + y5 + …. + yn-1) + 2(y2 + y4 + y6 + ….. + yn-2)]
What is the values of n in Boole’s rule?
From the closed Newton cotes formulas, we can say that the value of n is equal to 4 in Boole’s rule.
Comments