Testcase for UnivariateRealSolver.
Because Brent-Dekker is guaranteed to converge in less than the default
maximum iteration count due to bisection fallback, it is quite hard to
debug.
A polynomial spline function consists of a set of
interpolating polynomials and an ascending array of domain
knot points, determining the intervals over which the spline function
is defined by the constituent polynomials.