Abstract:The study of trigonometric series has started at the beginning of the nineteenth century. Joseph Fourier made the important observation that almost every function of a closed interval can be decomposed into the sum of sine and cosine functions. This technique to develop a function into a trigonometric series was published for the first time in 1822 by Joseph Fourier. The resulting series is nowadays called Fourier series. Since Fourier’s time, many different approaches to understand the concept of Fourier series have been discovered, each of which emphasizes different aspects of the topic. Some of the more powerful and elegant approaches are based on mathematical ideas and tools that were not available at the time Fourier completed his original work. Although the original motivation was to solve the heat equation for a metal plate, it later became obvious that the same technique could be applied to a wide variety of mathematical and physical problems and has many applications in electrical engineering, vibration analysis, acoustics, optics, signal treatment, image processing, etc. Despite the importance of Fourier series, the method used until now to compute them via computer algebra systems (CAS) is essentially based on the same principle as in Fourier’s time, i.e. by the evaluation of certain integrals. Unfortunately this technique is not completely successful for many functions. Although numeric values of the Fourier coefficients might be available, symbolic values are often not accessible. Modern CAS like Maple or Mathematica can compute such integrals in many cases for a given n∈Z\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${n \in \mathbb{Z}}$$\end{document}. However if one is interested in the Fourier coefficients for all n∈Z\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$${n \in \mathbb{Z}}$$\end{document}, then n is considered as a given symbolic variable and such integrals can be computed only in few cases. In this paper we introduce an algorithmic approach to compute those Fourier coefficients, involving differential equations of a particular form, and recurrence equations. This approach extrapolates the computation of the Fourier series for functions for which the computation of Fourier coefficients via the definition is out of reach for current CAS.A holonomic recurrence equation for an, i.e. a recurrence equation which is linear, homogeneous and has polynomial coefficients, can be written in operator notation as L(an) = 0. The operator L can be interpreted as a non-commutative polynomial via the commutator rule Nn − nN = N, N denoting the shift operator Nan = an+1. In the last section we show how our algorithm can be used to factorize such recurrence operators in certain cases.

A new and efficient method for the computation of Legendre coefficients

A Fast Precise Algorithm for Duhamel Integrals

A New Method for Numerical Abel-Inversion

A faster multipole Legendre-Chebyshev transform

Algorithmic Approach for Formal Fourier Series

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop theory

The Fourier series expansions of the Legendre incomplete elliptic integrals of the first and second kind

Numerical Methods for Fast Nonlinear Fourier Transformation, Part I: Exponential Runge-Kutta and Linear Multistep Methods

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

Fast and backward stable transforms between spherical harmonic expansions and bivariate Fourier series

A fast algorithm for the linear canonical transform

Rigorous numerical integration of algebraic functions

Fast Computation of the Discrete Fourier Transform Square Index Coefficients

Enhanced and Generalized One-Step Neville Algorithm: Fractional Powers and Access to the Convergence Rate

The Fourier-Legendre series of Bessel functions of the first kind and the summed series involving $\,_{2}F_{3}$ hypergeometric functions that arise from them

Ferrers Functions of Arbitrary Degree and Order and Related Functions

How to compute the constant term of a power of a Laurent polynomial efficiently

Abel inversion using fast Fourier transforms.

Efficient scaling and squaring method for the matrix exponential

Efficient numerical algorithms for multi-precision and multi-accuracy calculation of the error functions and Dawson integral with complex arguments

Large Sums of Fourier Coefficients of Cusp Forms