An Algorithm to Optimize the Calculation of the Fourth Order Runge-Kutta Method Applied to the Numerical Integration of Kinetics Coupled Differential Equations

作     者：Sadao Isotani Walter Maigon Pontuschka Seiji Isotani 

作者机构：Institute of PhysicsUniversity of Sao PauloSao PauloBrazil Institute of Mathematics and Computer ScienceUniversity of Sao PauloSão CarlosBrazil 

出 版 物：《Applied Mathematics》 (应用数学（英文）)

年 卷 期：2012年第3卷第11期

页      面：1583-1592页

学科分类：07[理学] 0701[理学-数学] 070101[理学-基础数学] 

主　　题：Coupled Differential Equations Runge-Kutta Kinetics Equations 

摘      要：The kinetic electron trapping process in a shallow defect state and its subsequent thermal- or photo-stimulated promotion to a conduction band, followed by recombination in another defect, was described by Adirovitch using coupled rate differential equations. The solution for these equations has been frequently computed using the Runge-Kutta method. In this research, we empirically demonstrated that using the Runge-Kutta Fourth Order method may lead to incorrect and ramified results if the numbers of steps to achieve the solutions is not “large enough. Taking into account these results, we conducted numerical analysis and experiments to develop an algorithm that determines the smallest non-critical number of steps in an automatic way to optimize the application of the Runge-Kutta Fourth Order method. This algorithm was implemented and tested in a variety of situations and the results have shown that our solution is robust in dealing with different equations and parameters.