Univariate Exclusion Solver

The methods of this family are based on a test which answers either

If we are in the first case, we remove the interval. Otherwise, the interval is splitted into two subintervals, to which we apply the test recursively, until a given precision is reached. Consequently, the remaining sub-domains might contain a root, but not guaranty is provided at this step.

The test implemented in thi exclusion solver uses the function

\[ \Delta_f (x_0, t) = |f (x_0) | - \sum_{i = 0}^{\deg (f)} \frac{|f^{(k)} (x_0) |}{k!} t^k . \]

We check that if $\Delta_f (x_0, r) > 0$, the interval $[x_0 - r, x_0 + r]$ does not contain a root of $f$, and thus can be removed. The corresponding function is

Seq<double> solve(f,Exclusion(),a,b);

It returns a sequence of type Seq<C>, where C is the type of the interval bounds a,b.

Implementation

See also:
synaps/usolve/Exclusion.h

SYNAPS DOCUMENTATION
logo