RESIDUE ANNIHILATION METHOD FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
Subject Areas : International Journal of Industrial Mathematics
Mehmet Pakdemirli
1
(Emeritus Professor, Mechanical Engineering Department, Manisa Celal Bayar University, Muradiye, Yunusemre, Manisa)
Keywords: Numerical Analysis, Residue Annihilation, Numerical Iteration Techniques, Convergence, Initial and Boundary Value Problems,
Abstract :
A new numerical algorithm is proposed for solving ordinary differential equations. The algorithm is names as Residue Annihilation Method (RAM). The method does not require transfer of the equation into a first order system of equations. For a k’th order nonlinear ordinary differential equation, a parametric solution containing k+2 parameters is assumed as an initial step. By imposing the compatibility conditions together with the annihilation of the residue and its first derivative, a nonlinear system with k+2 equations is obtained. Solving the system yields a recursive relation for the parameters. The assumed parameter values therefore vary at each integration step. Evaluating the parametric solution at each integration step yields the discrete numerical solution. A continuous approximate solution valid throughout the whole domain can also be expressed in terms of the Gamma Interval Functions. Sample ordinary differential equations up to third order derivatives are treated with the new method. The method can be applied to initial value problems directly and to boundary value problems when combined with shooting techniques. Depending on the assumed parametric solution, better convergence can be achieved to the real solution. The convergence rate for the algorithm is O(h2), h being the step size. By including higher order derivatives of the residue, convergence rate can be increased.
