Simpson’s Rule: Definition and Formula

One of the numerical approaches for evaluating the definite integral is Simpson's rule. To get the definite integral, we usually employ the fundamental theorem of calculus, which requires us to use anti-derivative integration techniques. However, in other cases, such as in Scientific Experiments, where the function must be calculated from observed data, finding the anti-derivative of an integral is difficult. In such situations, numerical approaches are utilized to approximate the integral. Trapezoidal rule, midpoint rule, and left or right approximation using Riemann sums are some of the other numerical methods used. We will go through Simpson's rule formula, the 1/3 rule, the 3/8 rule, and some examples in this section. 

Key takeaways: Simpson’s Rule, 1/3 rule, 3/8 rule, Simpson’s Rule Error, definite integrals, trapezoidal rule, integral, calculus, midpoint rule

Also read: Isosceles Triangle Theorems

---

Simpson’s Rule Formula

[Click Here for Sample Questions]

Simpson's rule can be used to approximate integrals, according to multiple sources. Quadratic polynomials are used to do this. In this case, parabolic arcs have taken the place of the trapezoidal rule's straight line segments. When it comes to finding approximate polynomials, Simpson's one-third rule can provide certain results. This is possible up to a cubic degree. It's crucial to note that the trapezoidal formula can help you find the area of objects that are shaped beneath a curve. Simpson's formula, on the other hand, is the way to go if one wishes to make those approximations better and more precise. It's also worth noting that Simpson’s rule parabolas are used to discover curve components. This indicates that the approximate area under the curve can be determined using the following formula, according to the Simpson definition:

Also Read:

---

Simpson’s Rule Derivation

[Click Here for Sample Questions]

Let us derive Simpson's rule, which involves splitting the area under the curve f(x) into parabolas to approximate the value of the definite integral ba f(x) d_x. Divide the interval [a, b] into n subintervals [x₀, x₁], [x₁, x₂], [x₂, x₃],..., [x _n−2, x _n−1], [x _n−1, x _n] each of width 'h', where xâ‚€ = a and xâ‚^TM= b.

Let us now approximate the area under the curve by assuming that every three consecutive points lie on a parabola. By sketching a parabola between the points x₀, x₁, and x₂, we may approximate the area under the curve between x₀ and x₂. Naturally, all three may not appear on the same parabola. However, let us attempt to draw a parabola through these three locations.

Let's make the y-axis of this parabola symmetric. After that, it looks like this:

Let's pretend that the parabola's equation is y = ax²+ bx + c. The definite integral is then used to approximate the area between x₀ and x₂:

Area between xâ‚€ and xâ‚‚ ≈ â‚‹â‚•∫Ê° (ax² + bx + c) dx

= (ax³/3 + bx²/2 + cx) â‚‹â‚•|Ê°

= (2ah³/3 + 0 + 2ch)

= h/3 (2ah² + 6c) ... (1)

Let's look at another point from the diagram above.

• f(xâ‚€) = a(-h)² + b(-h) + c = ah² - bh + c
• f(x₁) = a(0)² + b(0) + c = c
• f(xâ‚‚) = a(h)² + b(h) + c = ah² + bh + c

Now, f(xâ‚€) + 4f(x₁) + f(xâ‚‚) = (ah² - bh + c) + 4c + (ah² + bh + c) = 2ah² + 6c.

Substitute this in (1):

Area between xâ‚€ and xâ‚‚ ≈ h/3 (f(xâ‚€) + 4f(x₁) + f(xâ‚‚))

Similarly, we can see that:

Area between xâ‚‚ and xâ‚„ ≈ h/3 (f(xâ‚‚) + 4f(x₃) + f(xâ‚„))

Calculating the other areas in a similar way, we get

^b∫ₐ f(x) d^x

= h/3 (f(x₀) + 4f(x₁) + f(x₂))

+ h/3 (f(x₂) + 4f(x₃) + f(x₄))

+ ...

+ h/3 (f(xn−2n−2) + 4f(xn−1n−1) + f(xn))

≈h₃[f(x₀)+4f(x₁)+2f(x₂)+⋯+2f(x_n−2)+4f(x_n−1)+f(x_n)]h₃[f(x₀)+4f(x₁)+2f(x₂)+⋯+2f(x_n−2)+4f(x_n−1)+f(x_n)]

The like terms are combined here.

Hence we have derived Simpson's rule formula.

---

Simpson’s 1/3 Formula

[Click Here for Sample Questions]

The trapezoidal rule is extended by Simpson's 1/3rd rule, in which the integrand is approximated by a second-order polynomial. The Simpson rule can be determined in a variety of ways, including utilizing Newton's divided difference polynomial, Lagrange polynomial, and the coefficients technique. The Simpson's 1/3 rule is defined as follows:

∫_a^b f(x) d^x = h/3 [(y₀ + y_n) + 4(y₁ + y₃ + y₅ + …. + y_n-1) + 2(y₂ + y₄ + y₆ + ….. + y_n-2)]

Simpson's One-third Rule is the name of this rule.

Simpson’s 1/3 Rule for Integration

When we divide a tiny interval [a, b] into two halves, we can get a rapid approximation for definite integrals. As a result of splitting the interval, we get:

x₀= a, x₁= a + b, x₂ = b

Hence, the approximation can be written as;

∫_a^b f(x) d^x ≈ S₂ = h/3[f(x₀) + 4f(x₁) + f(x₂)]

S₂ = h/3 [f(a) + 4 f((a+b)/2) + f(b)]

Where h = (b – a)/2

This is the Simpson’s â…“ rule for integration.

---

Simpson’s 3/8 Rule

[Click Here for Sample Questions]

∫_a^b f(x) d^x = 3h/8 [(y₀ + y_n) + 3(y₁ + y₂ + y₄ + y₅ + …. + y_n-1) + 2(y₃ + y₆ + y₉ + ….. + y_n-3)]

Because it employs one extra functional value, this rule is more accurate than the normal technique. The composite Simpson's 3/8 rule, which is comparable to the generalized form, also exists for the 3/8 rule. Simpson's second rule of integration is known as the 3/8 rule.

---

Simpson’s Rule Error

[Click Here for Sample Questions]

It's important to note that, while employing Simpson's approach for definite integral calculation yields a more precise approximation, errors still occur. When this occurs, this is specified by n = 2; -(1/ 90) (b - 1 / 2) 5f (4) (ξ)

Here, ξ is some number that exists between a and b.

---

Things to Remember

• Simpson's rules are numerous approximations for definite integrals in numerical integration, named after Thomas Simpson (1710–1761). 
• It is called after Johannes Kepler in German and various other languages, who derived it in 1615 after seeing it used for wine barrels (barrel rule, Keplersche Fassregel). If f is a polynomial of up to a third-degree, the rule's approximation equality becomes precise.
• The composite Simpson's rule is obtained by applying the 1/3 rule to n equal subdivisions of the integration range [a, b]. Weights of 4/3 and 2/3 are alternated for points within the integration range. Simpson's 3/8 rule, sometimes known as Simpson's second rule, necessitates one extra function evaluation within the integration range and results in reduced error limits, but does not improve error order.
• Simpson's rules are substantially less efficient than trapezoidal rules when it comes to estimating the complete area of narrow peak-like functions. 
• To put it another way, the composite Simpson's 1/3 rule requires 1.8 times as many points as the trapezoidal rule to obtain the same level of accuracy. Simpson's 3/8 rule for composites is much less precise.

Also Read: