Решение систем дифференциальных уравнений неявными методами на с++

Необходимо написать методы для решения систем дифференциальных уравнений неявными методами на ЯП с++.

1) Метод Адамса-Мултона (используя метод простой итерации)
2) Неявный метод Рунге-Кутты (любая реализация, например метод Эйлера первого порядка)
3) Метод Адамса. Метод прогноза–коррекции

ВАЖНО.
Методы должны принимать входные параметры:
1) уравнение или несколько уравнений (если это система)
2) отрезок
3) шаг
4) начальные условия
5) точность вычисления (знаков после запятой)

Пример входных параметров:
y1' = 1.5*e^(2*x)-x^2-x-0.5
y2' = 3*x^2*y+x^2*e^x
y1(0) = 1
y2(0) = 1
начало отрезка - 0
конец отрезка - 3
шаг - 0.1
точность - 9 (девять знаков после запятой)

Данные методы необходимы для калькулятора систем дифференциальных уравнений
При необходимости есть парсер для уравнений на с++ и пример реализации явных методов для этого же калькулятора.

Цена договорная, возможно поэтапно, при необходимости, желаемый срок исполнения 2-3 дня

Исполнителю предоставлю все вспомогательные файлы, при необходимости.