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

# Numerical integration can be described as set of algorithms for calculating the numerical value of a definite integral. Definite integrals arise in many different areas and calculus is a tool

Extracts from this document...

Introduction Aim:

Numerical integration can be described as set of algorithms for calculating the numerical value of a definite integral. Definite integrals arise in many different areas and calculus is a tool for evaluating them; with numerous applications in science and engineering as well as mathematical analysis. However, calculus cannot always be applied; there are functions which do not have antiderivaties. One such example is ;this is an important function since it will be used in this coursework.

Using my knowledge of numerical integration I shall produce an approximation to the function ,seeing as it cannot be integrated. The area will be evaluated between the values of 0 and 1 (radians). The graph below gives a visual representation of the area which I shall be calculating.

Explanation:

It is often complicated to find the analytic solution to many differential equations. However, to our benefit there are many methods for finding the approximate solutions to differential equations. These methods are referred to as polynomials: the mid-point rule, trapezium rule and Simpson’s rule. Before explaining the methods in detail, we should note that all of these methods presented do not produce exact solutions, only approximate ones.

Lissaman R. (2004), suggests that the midpoint rule uses rectangles to approximate the area underneath a curve. Below is a diagram which makes use of the mid-point rule: In figure 1.0 five rectangles, each with the same width, are used to approximate the area under the graph of a function f(x) between x = 0 and x = 1.

The widths of the rectangles are donated h. The height of the first rectangle in the example is at the mid-point of the interval 0 to 0.2; which is 0.1 (represented by dashed line in figure 1.0).

Middle   On the other hand, the trapezium rule is an underestimate and each gap is a quarter to the one on its left. Likewise, the trapezium rule is also a second order polynomial and only underestimates when the curve is convex.       Notice, that the mid-point rule is a mirror image of the trapezium rule and vice versa. This is because the mid-point overestimates whilst the trapezium rule underestimates. Therefore, by using both methods we can confine the exact value between the two different polynomials. This means by obtaining values of Tn and Mn where n is a large number, we can find an accurate approximation to the integral as we know the exact value is between that of Tn and Mn.

Technology:

Most of this coursework will be carried out using Microsoft Excel; an electronic spreadsheet program used for organizing and manipulating data. The program is capable of working accurately up to 16 decimals places; this is also the amount I have chosen to work with in my calculations. Due to its manipulative ability Excel makes the handling of data very easy and saves an enormous amount of time. This is due to the fact that you don’t need to calculate everything; once the formula is entered and two consecutive calculations are complete, the cells can be dragged down and the answers required appear.

Excel is able to do this as it follows the formula and is judicious. The reason I chosen to use a program instead of a calculator is due to its ability to use 16 decimal places whereas a calculator can only give answers accurate to 9 decimal places.

Conclusion

n+Tn)/3. The average is twice as close to Mn as it is to Tn. This can be justified by the differences in the errors associated with the two polynomials. This is explained below:

Estimate the area of I=  using the trapezium rule and mid-point rule, only using 1 strip.

Trapezium rule:         Mid-point Rule:

h= (4-0)/1= 4           h= (4-0)/1= 4

 x 0 4 f(x) 0 16 x 2 f(x) 4

T1 = ½ * 4 * (f(0)+f(4)) = 32        M1= 4* (f(2)) = 16

The area of can be worked out from integration and is exactly 21 1/3. Therefore the absolute error in T1 is (32-21 1/3) = 10 2/3, whereas that in M1 is (21 1/3- 16) = 5 1/3. From this we can see that the absolute error in T1 is twice as much as that in M1. You can see from inspecting the errors in the trapezium sum and the midpoint sum that the midpoint sum is about twice as accurate as the trapezoidal sum and opposite in sign. This explains the weighting of the formula (2Mn+Tn)/3. By using the Simpson’s rule we will have cancellations of the errors and should thus get a much more accurate approximation.

The advantage of using the Simpson’s rule is that it’s a fourth- order polynomial. As mentioned before Simpson’s fits a parabola between successive triples of pairs. The absolute error is proportional to h4 so it is able to achieve a more accurate approximation to the curve, .

Absolute error ∝ h4 Absolute error = kh4

As a result the absolute error is Sn= kh4 whereas the absolute error in S2n= k(h/2)4. This means the errors reduces by a scale factor of 1/16 between each successive Sn values. In theory, the results I have obtained for Sn should be more accurate and reach the exact value in fewer calculations. This is seen to be the case by looking at the spreadsheet attached to this coursework. 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

# Related AS and A Level Core & Pure Mathematics essays

1. ## MEI numerical Methods

and g(x) and drag it down to 9 significant figures. A spreadsheet can be made so this method can once again be used to find the root/s of the equation between 0 and ?/2. Secant method: Earlier I claimed that using the Newton raphson method can't be used because differentiation

2. ## Investigation of circumference ratio - finding the value of pi.

Segment AC and Segment BC are radius r of circle. Angle ACD dependent by the n, n decide how many equilateral triangle in circle. So we use n to divide 360 degree (the degree of central angle degree), then we will got the degree of angle ACD.

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

Microsoft Excel can also helpful and therefore a spread sheet is set up: n xn g(xn) g'(x) 1 0.5 =0.5*B2^3-B2^2+1 =(3/2)*B2^2-2*B2 2 =C2 =0.5*B3^3-B3^2+1 =(3/2)*B3^2-2*B3 3 =C3 =0.5*B4^3-B4^2+1 =(3/2)*B4^2-2*B4 -number of iteration done. - value of x - output of x when - gradient of xn n xn g(xn) g'(xn)

2. ## Functions. Mappings transform one set of numbers into another set of numbers. We could ...

Adding and Subtracting ==> To add or subtract fractions, they must have the same denominator ==> Be careful when subtracting (especially where there are double negatives) ==> When the denominators have no common factor, their product gives the new denominator (multiply the numerator accordingly)

1. ## C3 Coursework: Numerical Methods

As you can see from the table there is no change of sign. As the graph of y=x�+10x�+4.8x+0.576 is a repeated root the change of sign method, and thus Change of Sign, cannot be used to find the root. Repeated roots such as y=x�+10x�+4.8x+0.576 cause the change of sign method to fail as there will be no change of sign.

2. ## Functions Coursework - A2 Maths

Therefore taking x1=3, the following iteration is constructed: x1 3.00000000 x2 2.90000000 x3 2.89511794 x4 2.89510652 x5 2.89510652 x6 2.89510652 x7 2.89510652 x8 2.89510652 x9 2.89510652 x10 2.89510652 Each value is done to eight decimal places by Microsoft Excel. As one can see, the iteration has appeared to converge to 2.89510652 in three steps.

1. ## Solving Equations Using Numerical Methods

Here is the table of results. I can now see that the root lies between -1.526 and -1.525. As I only intend to find the root to 3 decimal places I don't need to continue the method. To be certain that the root is between -1.526 and -1.525, I will calculate f(x)

2. ## Numerical integration coursework

You then multiply this by the height, which you can get by doing the function of the mid-point of each rectangle, and finally add up the area found of them. Hence you get this formula: Mn= h (f(mâ)+f(mâ)+…+f(mn)) I used this formula and many algorithms when working on Excel to • Over 160,000 pieces
of student written work
• Annotated by
experienced teachers
• Ideas and feedback to 