Blends in Maple

A blend of two Taylor series for the same smooth real- or complex-valued function of a single variable can be useful for approximation. We use an explicit formula for a two-point Hermite interpolational polynomial to construct such blends. We show a robust Maple implementation that can stably and efficiently evaluate blends using linear-cost Horner form, evaluate their derivatives to arbitrary order at the same time, or integrate a blend exactly. The implementation is suited for use with evalhf. We provide a top-level user interface and efficient module exports for programmatic use. This work was presented at the Maple Conference 2020. See www.maplesoft.com/mapleconference