• Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month
  1. 1
  2. 2
  3. 3
  4. 4
  5. 5
  6. 6
  7. 7
  8. 8
  9. 9
  10. 10
  11. 11
  12. 12
  13. 13

Solving Equations by numerical methods.

Extracts from this document...


Solving Equations by numerical methods. The fact that the vast majority of the numbers on the real number line are irrational (i.e. cannot be put in the form p/q, where p and q are integers) has an associated consequence that the vast majority of equations involving powers of x (i.e. polynomials) are insoluble by closed analytical techniques. This coursework's solution implies finding values of x say c1, c2-...c3 where f(c1)=0; f(c2)=0 and so on. Further transcendental equations and non- linear equations require numerical methods for their solution. It is to be noted that in no way are numerical methods inferior to analytical solutions, they are indeed the only practical solutions available, however on the down side no exact solution is possible and error bounds have to be placed on the solution given. In this coursework three main methods of solution from one-dimensional equations are given, each of the methods could be extended to the multi-dimensional case. The aim of this coursework is to show example of the methods in action, where these are successful and where they are not. A comparison between the methods will also be attempted. Method 1: Decimal Search SOLVING AN EQUATION WITH THIS METHOD Fig1. Shows the function f(x) = 12 ln(x) - x3/2. (Where ln(x) represents the natural logarithm of x.) To find the roots, the equation is: 12 ln(x) - x3/2 = 0 This equation can not be solved numerically because x appears inside the function ln(x) ...read more.


in the form x=g(x). A series of iterations attempts to improve on values of x; geometrically the method is relating on the intersection of the curve y=x with that of y=g(x). An initial value of x is chosen say x0. From this x value the value of the curve g(x) is obtained and a line drawn to the y=x curve, dropping down from this intersection point the value x1 is obtained from where the procedure continues as before. Eventually a point is reached simultaneous with the curve y=x and the curve g(x). In applying this method to the equation f(x) = 2ln(x-1) -x + 7.5 = 0, y = x f(x) = 2ln(x-1) - x + 7.5 x = 2ln(x - 1) + 7.5 = g(x) x = e^(0.5(x - 7.5)) + 1 = g(x) We immediately see the problem of this method. There are two choices for the iteration: 1. x n+1= 2ln(xn-1) + 7.5 2. xn+1= e^(0.5(xn-7.5)) + 1 1. Which of the choices of iteration for x should be used? Without a detailed analysis this question cannot be answered. However one criterion that could and should be used by the user is that for a successful iteration on an increasing function of x approaching the intersection with the curve y=x is that g'(x) should be less than 1 and greater than -1. For a decreasing function in this vicinity things are less clear cut. ...read more.


A simple program (see appendix) can be programmed into the calculator that generates the x-values, or alternatively spreadsheet software such as Excel can be used. Overall Conclusion. Decimal search is largely reliable, if not slow. It required little programming, and would only fail in situations where the use was not careful. However even in taking care it is possible to miss roots in highly oscillatory local situations. Newton Raphson also can present problems; one often has to conduct an analysis to find the initial point, (not always possible with more complicated equations) on the range of values of x one can get away with in the initial trial for convergence to a root. The conclusion here must also be that further research is required into this method. Fixed point estimation iteration presents real problems for the user. By experience he must select from a number of options the most likely iterative method (normally it appears to be safest to choose an iteration with a (1/n)th root in it. Also one cannot be sure which root if any your chosen iterate will converge too, one can however for increasing functions use the criteria of the gradient of the iterative function should be less than 1 in the vicinity of the root. In other examples such as xn+1=K (1-xn) chaotic functions arise. Also in the literature using complex numbers Julia sets with beautiful and chaotic boundaries arise. The world of fractals may also be mentioned here as an area to explore if more time was available. Daniel Antelo Pure2 Coursework St Charles Sixth Form College ?? ?? ?? ?? ...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)

    root is between -1 and 0 When I do the iterations starting at the point -1, it turns out like below. The green line is my original equation showing my real root and this clearly shows that the root that is being implied by the y = g(x)

  2. Marked by a teacher

    The Gradient Function

    5 star(s)

    2 4 2 4 8 8 3 9 18 12 4 16 32 16 The gradient here is equal to 4x - I can see this by investigating the values of x, as compared to the gradient. I predict that the gradient for 5 will be 20.

  1. Marked by a teacher

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

    3 star(s)

    But it will more precisely as it punish the insignificant value in the equation. The bigger the F stat the better. Compare the F-value in this equation with the first equation, which is 1.557e+004 and 1.404e+004. This means that the second equation works better than the first one, which gives opposite view from the R2.

  2. The open box problem

    For the ratio 3:1 I will start with the dimensions 9x3. x x 3 9-2x x 9 3-2x If I put the dimensions 9 and 3 into the original equation it comes out to be V=x(9-2x)(3-2x).

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

    This will also be explained in more detail in the interpretation. Tn Values of T Difference Ratio of Differences T1 0.88072085396282900 0.03875772638719590 T2 0.91947858035002500 0.00999934197580865 0.25799609285420 T4 0.92947792232583400 0.00252894688806504 0.25291133098391 T8 0.93200686921389900 0.00063439585807401 0.25085376884265 T16 0.93264126507197300 0.00015874172335739 0.25022503116480 T32 0.93280000679533000 0.00003969449396257 0.25005709351664 T64 0.93283970128929300 0.00000992419222789 0.25001432786251 T128 0.93284962548152100 0.00000248104805697 0.25000000000000

  2. Solving Equations Using Numerical Methods

    then the gradient is 0 and the next tangent is a horizontal line. As this is the case, it can't find the next root as the horizontal line continues forever and an overflow is created meaning that it will never reach the x-axis so it won't create a tangent back to the curve.

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

    -0.345008342 10 -0.27117 -0.34500834 -0.345008188 11 -0.27117 -0.34500819 -0.34500823 12 -0.27117 -0.34500823 -0.345008218 13 -0.27117 -0.34500822 -0.345008221 14 -0.27117 -0.34500822 -0.345008221 15 -0.27117 Root is -0.3450082 to 7d.p. Error bounds is -0.3450082�0.00000005. Root bounds is -0.34500825<x<-0.34500815. Check X Y -0.34500825 -0.345008213 -0.34500815 -0.34500824 X g'(x)

  2. C3 COURSEWORK - comparing methods of solving functions

    Therefore, we might miss the other two after finding the first one. This reason is that the two points are too close to each other. There are 3 root in the graph, 2 in the same interval [-3, -2] and 1 in the interval [0, 1].

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