Alvaro Meseguer - Fundamentals of Numerical Mathematics for Physicists and Engineers

(1.19)

so that ( 1.17) now leads to the secant method :

Secant Iteration:

(1.20)

It is common practice to start the indexing at by taking the initial values and , so that the first iterate is the same as the one obtained with the chord method. It should be clear that both chord and secant methods involve just the new evaluation per iteration, the remaining terms being previously computed (and stored) in the past.

Starting from Newton's code as a reference, the reader may easily program chord and secant algorithms and compare their efficiency with Newton's method when solving the test cubic equation starting from the same initial interval , as was done with Newton's and bisection methods. Figure 1.2b shows the resulting convergence history of the two methods. From the two sets of data of Figure 1.2b we may conclude that while the chord method seems to converge linearly, the secant methods converges superlinearly with an approximate order . It can be formally proved that the exact convergence order of the chord method is , whereas for the secant method the order is the golden ratio (see exercises at the end of the chapter).

A small point‐like object of mass initially at rest is released from the highest point of the surface A of height , as shown in the following figure. The object slowly starts to slide under the effects of gravity until it reaches point B at which it loses contact with the surface. The goal of this practical is to determine the abscissa of that point.

First, assuming there are no friction forces, show that the speed of the object at , i.e. still in contact with the surface, is

Then show that the horizontal velocity is given by the expression in the figure on the left. Imposing that there is no horizontal acceleration at point B (i.e. ) show that (excluding points where )

at the point where contact is lost.

Consider the family of surfaces shown in the figure on the right. Releasing the object from point , find the abscissas for , and 0.9.

As mentioned before, nonlinear equations usually involve parameters and constants coming from known external data. For example, Eq. ( 1.5) for the quantum energy levels depend on the well's depth and width , as well as on the mass of the particle . Therefore, it is obvious that changing the value of , , or in ( 1.5) will change the value of the corresponding energy levels . In particular, we may also expect that if the parameters are changed just by a tiny amount, then the changes in the solutions should also be very small.