On Implementation of Numerical Methods for Solving Ordinary Differential Equations in Computer Algebra Systems

Abstract

This paper presents an original package for investigating numerical solutions of ordinary differential equations, which is built in the Sage computer algebra system. This project is focused on a closer integration of numerical and symbolic methods while primarily aiming to create a convenient tool for working with numerical solutions in Sage. The package defines two new classes: initial problems and approximate solutions. The first class defines tools for symbolic computations related to initial problems, while the second class defines tools for interpolating values of symbolic expressions on an approximate solution and estimating the error with the use of the Richardson method. An implementation of the Runge–Kutta method is briefly described, with its main feature being the possibility of working with arbitrary Butcher tableaux and arbitrary numeric fields.