• 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.  5 star(s)

12.6 2 12 2.01 12.1203 12.06 2 12 2.001 12.012003 12.006 1 3 1.1 3.63 6.6 1 3 1.01 3.0603 6.06 1 3 1.001 3.006003 6.006 x x2 3x2 gradient 1 1 3 6 2 4 12 12 3 9 27 18 4 16 48 24 After looking through each value, I have observed a pattern.

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 ...

-1.203364 6 -1.203364 0 3x2+6x+0.5 -2.37593 0 -1.203364 After 5 converge we can see that the 2nd target root is to , with an error bound and solution bound. Target root 3 n xn f(xn) f'(x) f'(xn) f(xn)/f'(xn) xn+1 1 -3 -3.5 3x2+6x+0.5 9.5 -0.368421053 -2.631578 2 -2.631578 -0.764397142 3x2+6x+0.5

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

x1 -1.00000000 x2 -0.50000000 x3 -0.60000000 x4 -0.60270270 x5 -0.60270493 x6 -0.60270493 x7 -0.60270493 x8 -0.60270493 x9 -0.60270493 x10 -0.60270493 x1 -3.00000000 x2 -2.50000000 x3 -2.31914894 x4 -2.29294435 x5 -2.29240182 x6 -2.29240159 x7 -2.29240159 x8 -2.29240159 x9 -2.29240159 x10 -2.29240159 Therefore the roots to the equation x3-7x-4=0 are: x= 2.8951065, -0.6027049, -2.2924016 (all to seven decimal places)

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

To estimate the area you first need to find the area of each rectangle. You can do this by thinking logically; to get the width of each rectangle you need to do the range between the limits of the integral divided by the number of rectangles. • Over 160,000 pieces
of student written work
• Annotated by
experienced teachers
• Ideas and feedback to 