Abstract
Our purpose in this research is the development of higher order Runge-Kutta methods for solving stiff systems. We have developed methods of order five, six, and seven. We studied their stability Region and applications for solving stiff systems. Then we developed the corresponding implicit forms of these methods and we analyzed their stability and implementation for solving stiff systems.