A Continuous-Order Integral Operator for Maclaurin-type Reconstruction
Derek C. Braun
TL;DR
The paper develops a continuous-order generalization of the Maclaurin expansion by defining $\mathcal{T}[f](x)=\int_{0}^{\infty} \frac{D^{r}f(0)\,x^{r}}{\Gamma(r+1)}\,dr$ to reconstruct analytic functions from fractional-derivative data. It establishes an Euler–Maclaurin correction framework, introducing a correction series $\mathcal{E}[f](x)$ with explicit terms $\mathcal{E}_0,\mathcal{E}_1,\mathcal{E}_2$, to quantify and reduce the sum–integral discrepancy. Numerical experiments across a suite of analytic functions show that the uncorrected operator captures global structure with a modest offset, while incorporating the first three correction terms reduces the mean absolute error by about two orders of magnitude, to $\mathcal{O}(10^{-3})$. For monomials and Caputo-consistent boundary data, the method yields exact reconstructions, indicating a coherent framework for extending Maclaurin-type expansions to spectral representations in derivative order. The work points to a flexible, order-spectrum viewpoint in fractional calculus with potential implications for analytic reconstruction and spectral methods.
Abstract
I introduce a continuous-order integral analog of the Maclaurin expansion that reconstructs analytic functions from fractional derivative data. This operator replaces the discrete sum of integer-order derivatives in the classical Maclaurin expansion with an integral over fractional derivative orders, weighted by a Maclaurin-type kernel. I show that, under smoothness and decay assumptions, the discrepancy between this construction and the classical Maclaurin series is governed by the Euler--Maclaurin summation formula, and I define a corresponding correction series. Numerical experiments on a representative set of analytic functions show that the uncorrected operator reliably tracks the global structure of $f$, with a systematic, mostly constant offset and additional deviation near the origin. Adding the first three correction terms substantially reduces these discrepancies across the tested domains, lowering the mean absolute error by two orders of magnitude, from $10^{-1}$ to $10^{-3}$. For monomials, the order data collapse to a single atomic contribution; using Caputo-consistent boundary data, the operator reconstructs them exactly. Altogether, these results show that this continuous-order operator, together with its correction series, provides a coherent framework for extending the classical Maclaurin expansion to spectral representations in derivative order.
