Для динамической системы, у которых правые части являются квадратичными полиномами был составлен алгоритм аппроксимации разностной схемой 2-го порядка в нашей системе fdm for sage в системе компьютерной алгебры Sage. В этой системе начальная задача задается отдельно от метода ее решения. Ее решение, скажем, по явной схеме Эйлера задается как erk(pr), а по обратимой схеме - как cremona_scheme(pr). Обе функции возвращают элемент, принадлежащий классу Numsol, поэтому с новой схемой можно работать в этой системе так же, как со схемами Рунге-Кутты. При помощи этой реализации был проведен ряд компьютерных экспериментов, в ходе которых были найдены весьма неожиданные свойства этих схем.