• Join over 1.2 million students every month
• Accelerate your learning by 29%
• Unlimited access from just £6.99 per month
Page
1. 1
1
2. 2
2
3. 3
3
4. 4
4
5. 5
5
6. 6
6
7. 7
7
8. 8
8
9. 9
9
10. 10
10
11. 11
11
12. 12
12
13. 13
13
14. 14
14
15. 15
15
16. 16
16
17. 17
17
18. 18
18
19. 19
19
20. 20
20
21. 21
21
22. 22
22
23. 23
23
24. 24
24
25. 25
25
26. 26
26
27. 27
27
28. 28
28
29. 29
29

# Pure Mathematics 2: Solution of equation by Numerical Methods

Extracts from this document...

Introduction

Pure Mathematics 2: Solution of equation by Numerical Methods Introduction: In this coursework, I am going to solve equations by using the Numerical Methods. Numerical methods are used to solve equation that cannot be solved algebraically e.g. quadratic equations ax�+bx+c=0 can be solved using this formula: x= -b� V b� - 4ac 2a Therefore numerical methods would not be used for quadratic equations. I will be working with cubic equation because there is no formula to solve it. There are three methods, which I will be using: * Change of sign method * Newton-Raphson method * Rearranging f(x) = 0 in the form x = g(x) Change of sign method: This method is concerned with when a function crosses the x-axis, and by definition changes sign (+ and -). If we are looking the root of equation f(x) = 0. The point at which the curve crosses x-axis is the root. Once an interval where f(x) changes sign then the root must be in the interval. f(a) > 0 f(b) < 0 Therefore root must be between [a,b] f(a) < 0 f(b) > 0 Root is between the interval [a,b] To find the interval of each root for the equation, I'll be doing a decimal search first. Lets take the equation y = x� - 12x + 5 x -4 -3 -2 -1 0 1 2 3 4 f(x) -11 14 21 16 5 -6 -11 -4 21 There are 3 roots in this equation and they are in these intervals: [-4, -3] [ 0 , 1 ] [ 3 , 4 ] Now, I am going to use Interval Bisection to find out one of the root. This method is similar to decimal search but instead of dividing each interval into 10 parts, using interval bisection, only need to divide into 2, which is more simpler. f(x) = x� - 12x + 5 I am going to use the root in the interval [0,1] f(0) ...read more.

Middle

f(x(n)) f'(x(n)) x(n+1) 1 2 3 9 1.666667 2 1.666667 0.62963 5.333333 1.548611 3 1.548611 0.06804 4.194589 1.53239 4 1.53239 0.001218 4.044659 1.532089 5 1.532089 4.17E-07 4.04189 1.532089 6 1.532089 4.88E-14 4.041889 1.532089 7 1.532089 0 4.041889 1.532089 This shows the root is close to 1.532089 (6 d.p.) Failure of Newton-Raphson method: To test the failure of the Newton-Raphson method, I am going to choose another equation; y = x� -4x + 2 f(x) = x� - 4x + 2 xn+1 = xn - xn� - 4xn + 2 3xn� - 4 Let x1 = 1.1 This is close to the root in interval [0,1] and [1, 2], therefore using the Newton-Raphson method, it should converge to either one of these roots. x1 = 1.1 x2 = 1.1 - (1.1� - 4(1.1) + 2) 3(1.1�) - 4 = 1.1 - (1.331 - 4.4 + 2) 3.63 - 4 = 1.1 - (-1.069) -0.37 = 1.1 - 2.889189189 = -1.789189 x3 = -1.789189 - (-1.789189� - 4(-1.789189) + 2) 3(-1.789189�) - 4 = -1.789189 - (-5.7275456956 - (-7.156756) + 2) 9.603591833 - 4 = -1.789189 - 3.429210305 5.603591833 = -1.789189 - 0.611966468 = -2.401155 x4 = -2.401155 - (-2.401155� - 4(-2.401155) + 2) 3(-2.401155�) - 4 = -2.401155 - (-13.84396801 - (-9.60462) + 2) 17.296636 - 4 = -2.401155 - (-2.23934801) 13.296636 = -2.401155 - (-0.168414628) = -2.232740 n x(n) f(x(n)) f'(x(and)) x(n+1) 1 1.1 -1.069 -0.37 -1.789189 2 -1.789189 3.429208 5.603594 -2.401155 3 -2.401155 -2.239348 13.29664 -2.23274 4 -2.23274 -0.199539 10.95539 -2.214527 5 -2.214527 -0.002216 10.71238 -2.21432 6 -2.21432 -2.84E-07 10.70964 -2.21432 7 -2.21432 -5.33E-15 10.70964 -2.21432 8 -2.21432 0 10.70964 -2.21432 As you can see, the root is diverging towards the interval [-2, -3] instead of converging towards [0, 1]. This diagram shows clearly that when x1 is 1.1, the tangent diverts to -2, the reason that this happens is because the value that I chosen to draw the tangent is too close to the turning point of the curve, therefore when the tangent is drawn, it diverts to -2. ...read more.

Conclusion

However, the root that is found using this method is just as accurate as the other two methods, but using this method, the number of steps to work the root out is 20, and is the most steps taken in all three methods. * Newton-Raphson The Newton-Raphson method gave the fastest speed of convergence, although when calculating the root using the Newton-Raphson iteration formula takes a lot of time because xn, follows one after another. When the iterative formula is used, it can fail to find the root if the calculations show divergence (when x1 is too close to the turning point). Using the iterative formula, only 5 steps are needed to work out the root. * Rearranging method Rearranging the function of f(x) = 0 to g(x), using the iteration formula, and the staircase diagram, this converges to the root in [0, 1] but have to rearrange the formula again in order to work out the other two. Therefore, using this method may take longer to work out all three roots in the equation. Only 6 steps are needed to work out the root. All three methods are at 6 d.p. Uses of hardware and software: In this investigation, I had used quite a range of software packaging, which enables me to do complicated equations or graph instantly. A graph-drawing package (Autograph) benefited all the methods. I was able to create the graphs for the functions and equations. This package also enables me to zoom into a root to draw tangents. A spreadsheet software package (Excel) allowed me to customize the tables and implement formulae that helped to calculate the approximations at each value of x. On the graphics application, I could easily revise and correct the formulae to make changes where applicable. I was then able to import the formulae into the word processor and present them alongside the graphs and explanations. For hardware, I used a calculator to help me with complicated calculations to check for answers for Newton-Raphson and Rearranging formular. . ?? ?? ?? ?? Celene Leong 13.7 1 ...read more.

The above preview is unformatted text

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

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. ## MEI numerical Methods

False position method: Evaluating this formula, lets assume f(a) = -2 and f (b) = 10, lets assume C = 3, however the f(c) = -5, this means for the new equation would involve b and a would be replaced by c.

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

0.025129 0.007813 9 1.523438 0.025129 1.53125 -0.00827 1.527344 0.008436 0.003906 10 1.527344 0.008436 1.53125 -0.00827 1.529297 8.43E-05 0.001953 11 1.529297 8.43E-05 1.53125 -0.00827 1.530273 -0.00409 0.000977 12 1.529297 8.43E-05 1.530273 -0.00409 1.529785 -0.002 0.000488 13 1.529297 8.43E-05 1.529785 -0.002 1.529541 -0.00096 0.000244 14 1.529297 8.43E-05 1.529541 -0.00096 1.529419 -0.00044 0.000122

2. ## Numerical solution of equations, Interval bisection---change of sign methods, Fixed point iteration ---the Newton-Raphson ...

=(B8+C8)/2 =(D8-1)*(D8+2)*(D8-3)-1 =ABS(C8-B8)/2 =A8+1 =IF(E8>0,D8,B8) =IF(E8<0,D8,C8) =(B9+C9)/2 =(D9-1)*(D9+2)*(D9-3)-1 =ABS(C9-B9)/2 =A9+1 =IF(E9>0,D9,B9) =IF(E9<0,D9,C9) =(B10+C10)/2 =(D10-1)*(D10+2)*(D10-3)-1 =ABS(C10-B10)/2 =A10+1 =IF(E10>0,D10,B10) =IF(E10<0,D10,C10) =(B11+C11)/2 =(D11-1)*(D11+2)*(D11-3)-1 =ABS(C11-B11)/2 =A11+1 =IF(E11>0,D11,B11) =IF(E11<0,D11,C11) =(B12+C12)/2 =(D12-1)*(D12+2)*(D12-3)-1 =ABS(C12-B12)/2 =A12+1 =IF(E12>0,D12,B12) =IF(E12<0,D12,C12) =(B13+C13)/2 =(D13-1)*(D13+2)*(D13-3)-1 =ABS(C13-B13)/2 =A13+1 =IF(E13>0,D13,B13) =IF(E13<0,D13,C13) =(B14+C14)/2 =(D14-1)*(D14+2)*(D14-3)-1 =ABS(C14-B14)/2 =A14+1 =IF(E14>0,D14,B14)

1. ## I am going to solve equations by using three different numerical methods in this ...

0.000655 -1.80072 -1.800689697 -2.47371E-05 0.000655 -1.800704956 0.000315 -1.80072 -1.800704956 -2.47371E-05 0.000315 -1.800712585 0.000145 -1.80072 -1.800712585 -2.47371E-05 0.000145 -1.8007164 6.02E-05 -1.80072 -1.8007164 -2.47371E-05 6.02E-05 -1.800718307 1.77E-05 -1.80072 -1.800718307 -2.47371E-05 1.77E-05 -1.800719261 -3.5E-06 -1.800719 -1.800718307 -3.5084E-06 1.77E-05 -1.800718784 7.11E-06 In this spreadsheet, a and b are the two values of intervals.

2. ## Different methods of solving equations compared. From the Excel tables of each method, we ...

So, although we start with a value that is quite close to it, we fail to find the particular root we want. That is why the failure occurs. Method 3: Rearrangement This method aims to rearrange f(x)=0 in the form x=g(x)

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

+ f2n-1) + 2(f2 + f4 + f6 + .... + f2n-2)] However, seeing as I am carrying out the mid-point rule then the trapezium rule, I will be using a different version of Simpson's rule (2Mn + Tn)/3. This is simply due to the fact that it involves less arithmetic than using than using the formula stated in the page above.

2. ## C3 COURSEWORK - comparing methods of solving functions

Therefore, the root is 0.87939 to 5 decimal places. Lastly, we will use x=g(x) method to find the root in the interval [0, 1]: I need to rearrange the equation x³+3x²–3=0 in to the form of x=g(x), Let X1= 1 x³+3x²–3=0 n x 1 1 0.816496581 2 0.816496581 0.904741021 3

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