Successive substitution method

Errors are proportional to At for small At. When the trapezoid rule is used with the finite difference method for solving partial differential equations, it is called the Crank-Nicolson method. The implicit methods are stable for any step size but do require the solution of a set of nonlinear equations, which must be solved iteratively. The set of equations can be solved using the successive substitution method or Newton-Raphson method. See Ref. 36 for an application to dynamic distillation problems. 

The equation in cell B1 is copied into cells Cl though El. Then turn on the iteration scheme in the spreadsheet and watch the solution converge. Whether or not convergence is achieved can depend on how you write the equations, so some experimentation may be necessary. Theorems for convergence of the successive substitution method are useful in this regard. 

Convergence proofs are available under certain conditions (Finlayson, 1980), and once the iterate value gets close to the solution, the convergence is very rapid. This method is generally better than the successive substitution method, except for special cases, but sometimes a good initial guess is required. 

The convergence of the successive substitution method is slow, i.e., it may require many iterations for the sequence to converge. The Newton-Raphson method has a faster rate of convergence, which is given as follows ... 

Subroutine GETRHO is used to calculate the density of the mixture from the Sanchez-Lacombe EOS by a successive-substitution method. The line numbers start from 1 in this subroutine. 

Fig. 3.16 Residual history of the successive substitution method. The convergence depends highly on the relaxation factor D. For D=0.1, oscillations occurred.

The underlying principle of the successive substitution method is to arrange the given equation to the form... 

The iteration scheme for the successive substitution method to search for a root is defined as... 

Divide the equation by a and add x,- to the LHS and RHS (a similar procedure is used in Appendix A for the successive substitution method for solving nonlinear equations) to yield the equation... 

The classical Newton method has quadratic convergence properties whereas the successive substitution method has linear rate of convergence. However, because of the overshoot, the Newton iteration may fail to converge when the initial estimate is not a good estimate of the solution of the system of nonlinear equations. 

In the successive substitution method, sometimes because of a poor initial... 

The Wegstein method converges, even under conditions in which the method of a = (-c) does not. Moreover, it accelerates the convergence when the successive substitution method is... 

Example 1.1 Solution of the Colebrook Equation by Successive Substitution, Linear Interpolation, and Newton Raph on Methods. Develop MATLAB functions to solve nonlinear equations by the successive substitution method, the linear interpolation, and the Newton-Raphson root-finding techniques. Use these functions to calculate the friction factor from the Colebrook equation [Eq. (1.4)] for flow of a fluid in a pipe with e/Z> =10 and Njf, = 10. Compare these methods with each other. 

Successive Substitution method to find one root of a nonlinear equation. 

Listing 3.1. Code segment for the successive substitution method applied to Eq. 

The successive substitution method has successfully calculated flie two real roots of flic fourfli order polynomial in this case. However, flic convergence is slow, and follows what is known as linear convergence. This means that as the true solution is approached, the error at each iteration is some linear fraction of the error at flic previous iteration. Let s now study the convergence rate and the condition for convergence of the successive substitution method. Let s let X, =, where x is the true solution value and Sx =x - x is the error in... 

---