Differential equation solver

Solves y′ = f(x, y) with fourth-order Runge-Kutta or the adaptive Dormand-Prince method, supporting common functions and an optional comparison with the analytic solution.