Revision:OCR Core 3 - Simpson's rule - The Student Room
The Student Room

Revision:OCR Core 3 - Simpson's rule

Contents

12. Simpson's rule

The inaccuracy of the trapezium rule

In Core 2 the trapezium rule was introduced in order to enable the candidate to give estimates of definite integrals which they were not able to integrate accurately (or at all). It is relatively easy to see the ways in which the trapezium rule is inaccurate, merely drawing a diagram can show how the line of the trapezium will leave a large area which is not accounted for, and the larger the width of the trapezium, the greater the inaccuracy.

A geometric basis for Simpson's rule

Linear approximations of curves will generally be inaccurate as the curve will change in gradient, and the line does not. The only solution is to use shorter lines, however, in the context of the trapezium rule this is inaccurate.

A parabola can be used instead of a straight line. The parabola has a quadratic equation describing it, whose three coefficients can be defined to make the parabola fit the curve in question.

Suppose there is a section of curve, from "a", to "b", whose y-coordinates are  y_{0} , and  y_{2} respectively.

Let there be a midpoint, and the approximation parabola to have an equation  px^{2} + qx + r .

From the midpoint of the region allow there to be "h" units either way (and define an origin at the mid point), such that the area under this approximation parabola is:

 \int ^{h} _{-h} \left( px^{2} + qx + r \right) \,dx .

Consider the value of this integral:

 \left[ \frac{p}{3}x^{3} + \frac{q}{2}x^{2} + rx \right] ^{h} _{-h} =  \left( \frac{p}{3}h^{3} + \frac{q}{2}h^{2} + rh \right)  - \left( \frac{p}{3} \left( -h \right) ^{3} + \frac{q}{2} \left( -h \right) ^{2} - rh \right) =  \frac{2p}{3}h^{3} + 2rh .

There is no dependence on "q" in this.

 x = -h \implies y = y_{0}

 x = 0 \implies y = y_{1}

 x = h \implies y = y_{2}

Hence:

 y_{0} = p \left( -h \right) ^{2} + q \left( -h \right) + r

 y_{1} = r

 y_{2} = ph^{2} + qh + r

Hence:

 y_{0} + y_{2} = 2ph^{2} + 2r .

 y_{0} + y_{2} - 2y_{1} = 2ph^{2}

Hence:

 \frac{2p}{3}h^{3} + 2rh = \frac{h}{3}2ph^{2} + 2rh =  \frac{h}{3} \left( y_{0} - 2y_{1} + y_{2} \right) + 2hy_{1} =  \frac{h}{3} \left( y_{0} + y_{2} - 2y_{1} + 6y_{1} \right) =  \frac{h}{3} \left( y_{0} + 4y_{1} + y_{2} \right) = A .

Hence, for a single interval Simpson's rule is defined:

 \int ^{a} _{b} f(x) \,dx \approx \frac{h}{3} \left( y_{0} + 4 y_{1} + y_{2} \right)

Where  a - b = 2h .

This implies that there must be an even number of sub-intervals, and that intervals are equal in width.

Simpson's rule: general form

As with the trapezium rule, Simpson's rule benefits from more intervals, and therefore it is common to use a number of equal intervals.

Divide the interval from "b" to "a" into "2n" equal intervals of width "h".

 x_{0} = a, \ x_{1} = a + h, ..., \ x_{2n} = x_{0} + 2nh = b

To shorten the writing (as with the trapezium rule), allow the following to be true:  y_{k} = f(x_{k}) .

Hence, using a method of simply repeating Simpson's rule:

 \int ^{a} _{b} f(x) \,dx \approx \frac{h}{3} \left( y_{0} + 4y_{1} + y_{2} \right) +  \frac{h}{3} \left( y_{2} + 4y_{3} + y_{4} \right) + ... +  \frac{h}{3} \left( y_{2n - 2} + 4y_{2n - 1} + y_{2n} \right)

Hence:

 \int ^{a} _{b} f(x) \,dx \approx \frac{h}{3} \left( y_{0} + 4y_{1} + 2y_{2} + 4y_{3} + 2y_{4} + 4y_{5} + ... + y_{2n} \right) .

Due to the nature of Simpson's rule, the method of looking at the graph to see whether an underestimation or an overestimation has been made (as done with the trapezium rule) cannot be replicated.

As there is a lot of calculation involved it is important to be organised; using a table is helpful:

 \begin{array}{r|l|c} x_{i} & y_{i} & ky_{i} \\ \hline x_{0} & y_{0} & y_{0} \\ x_{1} & y_{1} & 4y_{1} \\ x_{2} & y_{2} & 2y_{2} \\ x_{3} & y_{3} & 4y_{3} \\ ... & ... & ... \\ x_{2n} & y_{2n} & y_{2n} \end{array} .

Make a note of "h" also.

Discussions Toggle
I have become nocturnal
started by: Jim_Reid
forum: Health
replies: 19
last post: 2 Minutes Ago
Vegetarianism/veganism and ' animals eating animals...'
started by: grey_chair
forum: Philosophy
replies: 109
last post: 3 Minutes Ago
Football Manager Society II
started by: Dalimyr
forum: Gaming
replies: 5166
last post: 3 Minutes Ago
Most annoying thing people do at University?
started by: The Phelps
forum: Student Life
replies: 246
last post: 5 Minutes Ago
Who are your idols from history?
started by: Gurmeet.Kapoor
forum: History
replies: 58
last post: 6 Minutes Ago
How to ask someone to live with you
started by: laparisienne
forum: Student Accommodation
replies: 12
last post: 10 Minutes Ago
Living in lots of different countries to learn languages
started by: language_man
forum: Travel
replies: 50
last post: 10 Minutes Ago
Edinburgh Accommodation 2012/2013
started by: gpostleth
forum: University of Edinburgh
replies: 5
last post: 10 Minutes Ago
The "I have an offer from Edinburgh" thread
started by: oxymoronic
forum: University of Edinburgh
replies: 83
last post: 12 Minutes Ago
Top 5 regrets of the dying
started by: Raving_Hippy
forum: Advice on Everyday Issues
replies: 0
last post: 15 Minutes Ago
Housemate is denying he told me he cheated
started by: Anonymous
forum: Friends, Family and Work
replies: 37
last post: 15 Minutes Ago
Are atheists good people?
started by: Lizzeraptor
forum: Religion
replies: 36
last post: 17 Minutes Ago
Increasing Tension Between the UK and Argentina over the Falkland Islands
started by: Craig_D
forum: News and Current Affairs
replies: 686
last post: 25 Minutes Ago
Vampire Diaries - Back tonight!
started by: Phelps
forum: TV
replies: 1111
last post: 26 Minutes Ago
Is this a good way to send a message to your kids?
started by: Agenda Suicide
forum: Society
replies: 38
last post: 29 Minutes Ago
Too nervous to try counselling again because I can't open up properly
started by: Anonymous
forum: Mental Health
replies: 8
last post: 30 Minutes Ago
*Updated Mk1* Project Opel Speedster, a fast car to beat an expensive junk.
started by: Herr
forum: Motoring
replies: 41
last post: 30 Minutes Ago
Good experiences with Citalopram?
started by: AJ2890
forum: Mental Health
replies: 15
last post: 31 Minutes Ago
High on coffee/caffeine all the time, could it affect my health?
started by: Raving_Hippy
forum: Food and Drink
replies: 7
last post: 31 Minutes Ago
Malia 2012!!
started by: CharlotteM
forum: Travel
replies: 11
last post: 32 Minutes Ago
Article Updates Toggle
Contact Us | Site Rules | Staying Safe on TSR | Advertising | Staff Blog | Essays & Coursework | Terms & Conditions | Top
Customise your TSR | Life Advice | Hobbies and Interests | Debate and Current Affairs | Study Help | University and University courses
Universities and HE Colleges | Careers, Employment and Gap Years | General Discussion

Customise your TSR