Abstract:In this paper, a general technique for solving nonlinear, two-point boundary-value problems is presented; it is assumed that the differential system has ordern and is subject top initial conditions andq final conditions, wherep+q=n. First, the differential equations and the boundary conditions are linearized about a nominal functionx(t) satisfying thep initial conditions. Next, the linearized system is imbedded into a more general system by means of a scaling factor α, 0≤α≤1, applied to each forcing term. Then, themethod of particular solutions is employed in order to obtain the perturbation Δx(t)=αA(t) leading from the nominal functionx(t) to the varied function$$\tilde x$$(t); this method differs from the adjoint method and the complementary function method in that it employs only one differential system, namely, the nonhomogeneous, linearized system.The scaling factor (or stepsize) α is determined by a one-dimensional search starting from α=1 so as to ensure the decrease of the performance indexP (the cumulative error in the differential equations and the boundary conditions). It is shown that the performance index has a descent property; therefore, if α is sufficiently small, it is guaranteed that$$\tilde P$$<P. Convergence to the desired solution is achieved when the inequalityP≤ɛ is met, where ɛ is a small, preselected number.In the present technique, the entire functionx(t) is updated according to$$\tilde x$$(t)=x(t)+αA(t). This updating procedure is called Scheme (a). For comparison purposes, an alternate procedure, called Scheme (b), is considered: the initial pointx(0) is updated according to$$\tilde x$$(0)=x(0)+αA(0), and the new nominal function$$\tilde x$$(t) is obtained by forward integration of the nonlinear differential system. In this connection, five numerical examples are presented; they illustrate (i) the simplicity as well as the rapidity of convergence of the algorithm, (ii) the importance of stepsize control, and (iii) the desirability of updating the functionx(t) according to Scheme (a) rather than Scheme (b).

Root finding techniques that work

A novel Newton-Raphson style root finding algorithm

A robust iterative family for multiple roots of nonlinear equations: Enhancing accuracy and handling critical points

A Parallel Algorithm For Finding All The Roots Of A Polynomial

On the stability analysis of numerical schemes for solving non-linear polynomials arises in engineering problems

Annealing approach to root-finding

Nonlinear Equations Solving with Intelligent Optimization Algorithms: A Survey.

Root-finding: from Newton to Halley and beyond

Efficient Multiplicative Calculus-Based Iterative Scheme for Nonlinear Engineering Applications

Annealing approach to root finding

Optimal fourth- and eighth-order iterative methods for solving nonlinear equations with basins of attraction

An Optimal Fourth Order Derivative-Free Numerical Algorithm for Multiple Roots

Solving a Class of Nonlinear Eigenvalue Problems by Newton's Method

Survey of Real Root Finding of Univariate Polynomial Equation in CAGD/CG

Fifth-order Iterative Method for Finding Multiple Roots of Nonlinear Equations

From Halley to Secant: Redefining root finding with memory‐based methods including convergence and stability

General technique for solving nonlinear, two-point boundary-value problems via the method of particular solutions

Robust Numerical Methods for Nonlinear Regression

On a Stable Multiplicative Calculus-Based Hybrid Parallel Scheme for Nonlinear Equations

A reliable technique for nonlinear evolution equations with applications in mathematical physics

An effective approach for solving a class of nonlinear singular boundary value problems arising in different physical phenomena