• Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month

Numerical Methods coursework

Extracts from this document...

Introduction

Numerical Methods Coursework

Numerical Integration

The Problem

Integration means finding the area underneath a particular region of a function. At my current knowledge of maths I am not able to integrate various functions. Therefore I am going to use knowledge of numerical methods to produce an approximation to an area which does not have an analytic solution.

In my coursework I will integrate the function:

image29.png

The graph below shows the graph of the function.

image30.jpg

The red arrow determines the region between 0.25 and 1.25, which then leads to the integral:

image37.png

I can not solve this problem using the knowledge of C1 and C2, because I am not able to integrate cos(x) yet. Due to this I suggest that this problem will be appropriate for numerical solution.

The Approximation Rules

To solve this problem, I am going to use knowledge of numerical integration studied in the “Numerical Methods” textbook. The approximate methods of definite integrals may be determined by numerical integration using:

1. The Trapezium Rule:

The Trapezium Rule divides the area underneath the curve into trapeziums. We can then use the formula image45.png (where a, b are the bases and h is the height of the trapezium) to estimate the area.

image61.jpgimage54.pngimage00.png

...read more.

Middle

image27.pngimage28.pngimage26.png

image03.pngimage02.png

image04.png

image39.jpgimage07.pngimage08.pngimage05.pngimage06.pngimage10.pngimage09.png

image13.pngimage13.pngimage11.png

image15.pngimage14.pngimage16.png

Error Analysis

Because all three rules are just estimates, there is always an error involved. In this section I will give an estimate of how much error there is in my solution of the problem.

Firstly I want to refer to my graph of my function on page 1. The section of the curve is concave that means that the Trapezium Rule gives an overestimate and the Midpoint Rule gives an underestimate. For that reason the estimates of the Trapezium Rule are decreasing and they will reach the solution from above. On the other side the Midpoint Rule estimates are increasing and they will reach the solution from below.

Error in the Trapezium Rule

In the Trapezium Rule the error is proportional to h2. We add a constant k so that:

image40.png

In general, if image41.png is the absolute error in image33.png then there is a constant (k) so that:

image42.png

Where h is the strip width corresponding to n strips.

So the Trapezium Rule with 2n strips has a strip width of image43.png, such that:

image44.png

That shows that halving h or doubling n will reduce the error by a factor of 0.25.

Therefore the “error multiplier” is 0.25.

...read more.

Conclusion

This can be proved by looking on the diagram below.

image19.pngimage18.pngimage17.pngimage19.png

image20.pngimage20.pngimage20.pngimage20.png

image22.png

The distance between image50.png and image51.png is 4 times smaller than the distance between image52.png and image53.png, where X represents the actual solution. To get the next distance between image51.png and image55.png we multiply the “multiplier” 0.25 with the distance between image50.png and image51.png.

This development leads to improved solution by extrapolation, as stated in the “Formula Application” section.

Error in my solution

The error gets smaller the more strips are used as shown above. Therefore:

Rule

Error

T2, M2

image56.png

T4, M4

image57.png

T8, M8

image58.png

T4096, M4096

image59.png

---

---

S2

image57.png

S4

image60.png

S256

image62.png

Interpretation

 My solution to the integral image63.png is: 0.885909(6 decimal places)

The solution refers to:        T4096=        0.885909120684554

M4096=        0.885909049759261

S256=        0.885909075135303

The fact that the Trapezium Rule and the Midpoint Rule give the same solution till the sixth decimal place guaranties that my solution is valid. Simpson’s Rule, the weighted average, proves the solution as well.

Therefore: My solution is proved by 3 different rules.

I can say with guarantee that there is an error involved, because all three rules are just approximation rules. But I can say that my solution is not far away from the actual solution as the error for the rules with 4096 stripes is very small (actual value stated in the “Error Analysis” section).

...read more.

This student written piece of work is one of many that can be found in our AS and A Level Core & Pure Mathematics section.

Found what you're looking for?

  • Start learning 29% faster today
  • 150,000+ documents available
  • Just £6.99 a month

Not the one? Search for your essay title...
  • Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month

See related essaysSee related essays

Related AS and A Level Core & Pure Mathematics essays

  1. Marked by a teacher

    The Gradient Function

    5 star(s)

    � - 3x� = 3(x� + 2hx + h�) - 3x� = 3x� + 6hx + 3h� - 3x�= 6hx + 3h� x +h-x h h h = 6x + 3h = 6x axn value gradient function x2 2x 2x2 4x 3x2 6x Here, I have noticed a trend in these data.

  2. MEI numerical Methods

    The difference between the two methods is that the secant method uses two approximations which are the most recent approximations of the root. The false position method also uses two approximations; it uses the most recent approximation of the root and the most recent root which has the opposite sign of the most recent approximation of the root.

  1. Solving Equations. Three numerical methods are discussed in this investigation. There are advantages and ...

    0.03125 6 0.6875 0.002319336 0.71875 -0.049697876 0.703125 -0.023702621 0.015625 7 0.6875 0.002319336 0.703125 -0.023702621 0.6953125 -0.010694265 0.0078125 8 0.6875 0.002319336 0.6953125 -0.010694265 0.69140625 -0.004188031 0.00390625 9 0.6875 0.002319336 0.69140625 -0.004188031 0.689453125 -0.000934478 0.001953125 10 0.6875 0.002319336 0.689453125 -0.000934478 0.688476563 0.000692398 0.000976563 11 0.688476563 0.000692398 0.689453125 -0.000934478 0.688964844 -0.000121048 0.000488281 12

  2. Numerical integration can be described as set of algorithms for calculating the numerical value ...

    The trapezium rule makes use of trapeziums instead of rectangles. Like the mid-point rule it is easier to describe by using a diagram: Figure 1.1 shows 4 trapezia, each one of length one (h), approximating the area under the curve, between 1 and 5.

  1. C3 Coursework: Numerical Methods

    This is illustrated by the graph below. This graph shows the equation of the line y=x�+10x�+4.8x+0.576 The roots appear to be between -1 and 0. By taking increments of 0.1 between -1 and 0 it will be possible to use decimal search to attempt to look for a change in sign.

  2. Functions Coursework - A2 Maths

    x f(x) 1.87930 -0.0006474767 1.87931 -0.0005715231 1.87932 -0.0004955684 1.87933 -0.0004196125 1.87934 -0.0003436555 1.87935 -0.0002676974 1.87936 -0.0001917381 1.87937 -0.0001157777 1.87938 -0.0000398162 1.87939 0.0000361464 1.87940 0.0001121102 The root therefore lies in the interval [1.87938,1.87939]. The next table show values of f(x) at values of x from 1.87938 to 1.87939, with intervals of 0.000001.

  1. Numerical integration coursework

    make an approximation from Mâ up to Mââ, these algorithms made it much easier and quicker to come to a final answer thanks to the ability to just simply pull down the box with the algorithm in to transfer it to other boxes.

  2. C3 COURSEWORK - comparing methods of solving functions

    22 0.876 -0.02565 23 0.877 -0.01809 The root lies between 0.879 and 0.880 24 0.878 -0.01051 25 0.879 -0.00293 26 0.8791 -0.00217 27 0.8792 -0.00141 28 0.8793 -0.00065 29 0.87931 -0.00057 30 0.87932 -0.0005 31 0.87933 -0.00042 The root lies between 0.8793 and 0.8794 32 0.87934 -0.00034 33 0.87935 -0.00027

  • Over 160,000 pieces
    of student written work
  • Annotated by
    experienced teachers
  • Ideas and feedback to
    improve your own work