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

Numerical solution of equations

Extracts from this document...

Introduction

Pure Mathematic 2 Coursework Numerical solution of equations By Michael Pang I am going to show 3 of the numerical methods for solving the equation which cannot be solved algebraically. They are interval estimation, fixed point estimation and Newton-Raphson method. Those of these numerical methods are used when algebraic ones are not available. When you found any equations which cannot be solved algebraically, probably you will draw the graph and see where the roots are. However, the other problem is found that you cannot get all roots accuracy. Actually, the numerical methods cannot find the exactly root but the answers are more accuracy than sketching the graphs. Therefore, we usually provide the answer to 5 or 6 decimal places depending on what the questioner needs. Finally, we check the answer by setting the lower bound and upper bound to see whether it has sign change or not. Even if the three numerical can solve the equation non-algebraically, they have its advantages and disadvantages. And now I am going to show how these methods works and their problems by using the equation F(x) = x�-9x+3. Interval estimation Assume that the roots of the equation F(x) = x�-9x+3, and I am looking for the roots which F(x) = 0. The roots of the equation are the values of x which the graph of y = x�-9x+3 crosses the x axis. By using the computer, I recognise that there are 3 roots on the graph and the value of these 3 roots. However, if we cannot use the computer, we have to use one of the numerical methods - interval estimation. ...read more.

Middle

They are - The curve touches the x-axis - The roots are close together - There is a discontinuity The curve touches the x-axis In the above graph, the curve touches the x-axis. Therefore, the sign change could not be found and none of the methods of interval estimation can be used. The roots are close together In the graph, there 3 roots close to together between the interval (0,1). In decimal search, F(0.1) = 0, so that 1 of the root is 0.1 and so further interval cannot be found. In addition, interval bisection and linear interpolation would be unaware of the existence of other roots since the 3 roots are too close together. There is a discontinuity in F(x) This equation y = 1/(x-1.1) has no root, but all change of sign methods will converge on a false root at x = 1.1 Avoid problem To avoid those problems occurring, it is important to start by sketching the graph and therefore the problem will be found. Fixed point estimation It is much different from the method of interval bisection. It involves an iterative process rather than finding interval within the roots must lie. Generally, iterative process means that a sequence of numbers by continued repetition of the same procedure. But this method has some limited which will be introduced when the process are showing. How does it method work? Before starting to use this method, the equation F(x) = 0 has to be rearranged into the form of x = g(x). ...read more.

Conclusion

The other methods to show they are right, I should check their upper bound and lower bound where I set it in Fixed point estimation. And now I show it again to make it clear. - The largest root : F(2.816915) = 0.0000142, F(2.816905) = -0.00013403 - The medium root : F(0.337615) = -0.00005233 F(0.337605) = 0.00003425 - The smallest root : F(-3.154525) = -0.00004153 F(-3.154515) = 0.00016701 Advantage Newton-Raphson method is fast and convenient, the processes are short and that why I didn't use a lot of presentation to show my understanding. It also gives an extremely rapid rate of convergence. I only use 4 to 5 steps for each of the root, so I find that it was much easier to finish than the first 2 methods. Disadvantage Even this method is the fastest and the most convenient, it still has several problems. Firstly, in my case, if I didn't choose the right initial value, I couldn't find the root by only 4 to 5 steps or even faster in other case. However, if the first value I h chose is not close to the root, or near the turning point, a turning point of y = F(x), the iteration may diverge, or converge to another root. And if the initial point is at a stationary point, F'(X) = 0 so the method cannot proceed. Secondly, all numerical methods for solving equations cannot be proceeded when the function is discontinuous function as it doesn't have any real roots. Thirdly, the function is not defined over the whole of. The tangent at the graph may meet the axis at a point outside the domain of the function. ...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

See related essaysSee related essays

Related AS and A Level Core & Pure Mathematics essays

  1. Marked by a teacher

    C3 Coursework - different methods of solving equations.

    5 star(s)

    Because this only works out the value to 3 decimal places, I will use the Newton Raphson Formula to work it out to 5 decimal places. I will summarise it in a table like previously So I have to work out F' of my equation.

  2. Marked by a teacher

    The Gradient Function

    5 star(s)

    2.00000000000000 4.1 2.024845673 0.248457 4 2.00000000000000 4.01 2.002498439 0.249844 4 2.00000000000000 4.001 2.000249984 0.249984 2 1.41421356237310 2.1 1.449137675 0.349241 2 1.41421356237310 2.01 1.417744688 0.353113 2 1.41421356237310 2.001 1.414567072 0.353509 1 1.00000000000000 1.1 1.048808848 0.488088 1 1.00000000000000 1.01 1.004987562 0.498756 1 1.00000000000000 1.001 1.000499875 0.499875 x x 0.5 Gradient 1 1

  1. Marked by a teacher

    Estimate a consumption function for the UK economy explaining the economic theory and statistical ...

    3 star(s)

    (0.108) (0.097) Figure 4 Table 5 EQ( 5) Modelling LC by OLS (using project12min.xls) The estimation sample is: 1950 to 2002 Coefficient Std.Error t-value t-prob Part.R^2 LY_1 0.999415 0.0004450 2246. 0.000 1.0000 sigma 0.0411118 RSS 0.0878894078 log-likelihood 94.4484 DW 0.46 no. of observations 53 no. of parameters 1 mean(LC)

  2. Am going to use numerical methods to solve equations that can't be solved algebraically

    -2 x2= -1.756756757 x3= -1.607996214 x4= -1.63240644 x5= -1.639214373 x6= -1.639214372 x7= -1.639214336 x8= -1.639214311 x9= -1.639214268 x10= -1.639214194 x11= -1.639214172 x12= -1.63921403 x13= -1.639213786 x14= -1.639212651 x15= -1.639212651 x16= -1.639211421 x17= -1.639209325 x18= -1.63919953 x19= -1.63919953 x20= -1.639188919 x21 = -1.639170721 .

  1. C3 Coursework: Numerical Methods

    The fastest procedure in terms of speed of convergence was the Newton Raphson Method. To find an estimate of the root to 5 decimal places using the Newton Raphson method it took only 5 iterations. As with the other methods, this method was made easier using software.

  2. MEI numerical Methods

    must equal negative and f(b) must be positive. In the question we are provided with an interval which is 0 < x < ?/2. Lets assume K = 1 F(0) = -1 F(?/2) = a very large positive number. Our current interval is 0 < x < ?/2, in order to get a better approximation we must find the midpoint of the interval (h), and do the f(h).

  1. Solutions of equations

    Answer Root B: 0.37282 (5.d.p.) Finding Root C X0 = 1 X1 = 2.12979 X2 = 1.62013 X3 = 1.40784 X6 = 1.34816 Answer Root C: 1.34846(5.dp.) Newton-Raphson Failure For the Newton-Raphson method to fail the X0 starting position must be between a turning point and the root.

  2. In this coursework, I am going to solve equations by using the Numerical Methods. ...

    So the root are in [-2,-1] interval. Now I will calculate the values between -2 and -1 so that I can get a value to 1 decimal place. F(x) -2 -1.9 -1.8 -1.7 -1.6 -1.5 -1.4 -1.3 y -20.4 -14.843 -10.109 -6.1048 -2.7468 0.04375 2.33875 4.20431 sign - - - - - + + + Now we

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