Numerical Algorithm for Design of Stability Polynomials for the First Order Methods

Abstract

The algorithm for coefficients determination for stability polynomials of degree up to m = 35 is developed. The coefficients correspond to explicit Runge-Kutta methods of the first accuracy order. Dependence between values of a polynomial at the points of extremum and both size and form of a stability domain is shown. Numerical results are given.