Рассматривается задача описания функции Ляпунова для аппроксимации области притяжения невозмущённого решения системы дифференциальных уравнений. Предлагаются алгоритмы для вычисления значений функции Ляпунова, и её производной в силу системы без использования их аналитических выражений. На конкретных примерах показывается, что использование функции Ляпунова, определенной таким образом, позволяет аппроксимировать область притяжения практически с любой степенью точности.