Численное интегрирование методом Симпсона с оценкой погрешности по методу Рунге

Для вычисления интеграла и построения интегральной кривой необходимо ввести исходную функцию, все начальные условия и выбрать количество разбиений.
На графике красной кривой обозначена интегральная кривая, зелеными линиями — разбиения.