In this study,we construct a family of single root finding method of optimal order four and then generalize this family for estimating of all roots of non-linear equation simultaneously.Convergence analysis proves tha...In this study,we construct a family of single root finding method of optimal order four and then generalize this family for estimating of all roots of non-linear equation simultaneously.Convergence analysis proves that the local order of convergence is four in case of single root finding iterative method and six for simultaneous determination of all roots of non-linear equation.Some non-linear equations are taken from physics,chemistry and engineering to present the performance and efficiency of the newly constructed method.Some real world applications are taken from fluid mechanics,i.e.,fluid permeability in biogels and biomedical engineering which includes blood Rheology-Model which as an intermediate result give some nonlinear equations.These non-linear equations are then solved using newly developed simultaneous iterative schemes.Newly developed simultaneous iterative schemes reach to exact values on initial guessed values within given tolerance,using very less number of function evaluations in each step.Local convergence order of single root finding method is computed using CAS-Maple.Local computational order of convergence,CPU-time,absolute residuals errors are calculated to elaborate the efficiency,robustness and authentication of the iterative simultaneous method in its domain.展开更多
N-body simulations of the Sun, the planets, and small celestial bodies are frequently used to model the evolution of the Solar System. Large numbers of numerical integrators for performing such simulations have been d...N-body simulations of the Sun, the planets, and small celestial bodies are frequently used to model the evolution of the Solar System. Large numbers of numerical integrators for performing such simulations have been developed and used;see, for example, [1,2]. The primary objective of this paper is to analyse and compare the efficiency and the error growth for different numerical integrators. Throughout the paper, the error growth is examined in terms of the global errors in the positions and velocities, and the relative errors in the energy and angular momentum of the system. We performed numerical experiments for the different integrators applied to the Jovian problem over a long interval of duration, as long as one million years, with the local error tolerance ranging from 10-16 to 10-18.展开更多
摘要In this study,we construct a family of single root finding method of optimal order four and then generalize this family for estimating of all roots of non-linear equation simultaneously.Convergence analysis proves that the local order of convergence is four in case of single root finding iterative method and six for simultaneous determination of all roots of non-linear equation.Some non-linear equations are taken from physics,chemistry and engineering to present the performance and efficiency of the newly constructed method.Some real world applications are taken from fluid mechanics,i.e.,fluid permeability in biogels and biomedical engineering which includes blood Rheology-Model which as an intermediate result give some nonlinear equations.These non-linear equations are then solved using newly developed simultaneous iterative schemes.Newly developed simultaneous iterative schemes reach to exact values on initial guessed values within given tolerance,using very less number of function evaluations in each step.Local convergence order of single root finding method is computed using CAS-Maple.Local computational order of convergence,CPU-time,absolute residuals errors are calculated to elaborate the efficiency,robustness and authentication of the iterative simultaneous method in its domain.
摘要N-body simulations of the Sun, the planets, and small celestial bodies are frequently used to model the evolution of the Solar System. Large numbers of numerical integrators for performing such simulations have been developed and used;see, for example, [1,2]. The primary objective of this paper is to analyse and compare the efficiency and the error growth for different numerical integrators. Throughout the paper, the error growth is examined in terms of the global errors in the positions and velocities, and the relative errors in the energy and angular momentum of the system. We performed numerical experiments for the different integrators applied to the Jovian problem over a long interval of duration, as long as one million years, with the local error tolerance ranging from 10-16 to 10-18.