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 Computer Algebra Package for Polynomial Sequence Recognition

A Maple Package for Computing Groebner Bases for Linear Recurrence Relations

On Calculating Partial Sums of Multiple Numerical Series by Methods of Computer Algebra

Ansatz in a Nutshell: A comprehensive step-by-step guide to polynomial, $C$-finite, holonomic, and $C^2$-finite sequences

Algorithmic Approach for Formal Fourier Series

Computer Algebra, Power Series and Summation

Polynomial Quotient Rings and Kronecker Substitution for Deriving Combinatorial Identities

Discovery of Exact Equations for Integer Sequences

A Combinatorial Approach to Frobenius Numbers of Some Special Sequences (Complete Version)

The NumericalCertification package in Macaulay2

Automated Search for Conjectures on Mathematical Constants using Analysis of Integer Sequences

A package on orthogonal polynomials and special functions

Exhibiting Randomness in Arithmetic using Mathematica and C

Algorithmic work with orthogonal polynomials and special functions

Automated Symbolic and Numerical Testing of DLMF Formulae using Computer Algebra Systems

An algebraic combinatorial approach to Sylvester's denumerant

A Polynomial Time Algorithm for Calculating Fourier-Dedekind Sums

Algorithm-assisted discovery of an intrinsic order among mathematical constants

The sum of a hybrid arithmetic function over a sparse sequence

A Polynomial Ring Connecting Central Binomial Coefficients and Gould's Sequence

The MatrixSchubert package for Macaulay2