Fractions of Recurrence Operators for Generalized Fourier Series in Classical Orthogonal Polynomials
Pith reviewed 2026-05-07 10:36 UTC · model grok-4.3
The pith
The recurrence for coefficients in classical orthogonal polynomial series is the numerator of a fraction of recurrence operators.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The linear recurrence equation for the coefficients is interpreted as the numerator of a fraction of linear recurrence operators. This provides a simple and unified view of previous algorithms computing these recurrences, using a noncommutative Euclidean algorithm as the engine. The effectiveness is shown through various examples.
What carries the argument
Fractions of recurrence operators in a noncommutative ring, with the target recurrence as the numerator, and the noncommutative Euclidean algorithm to compute them.
If this is right
- Algorithms for computing recurrences from differential equations can be unified under this operator fraction framework.
- The noncommutative Euclidean algorithm serves as an efficient engine for these computations.
- The method applies directly to series in any classical orthogonal polynomial basis.
- Examples demonstrate practical computation of such recurrences for specific differential equations.
Where Pith is reading between the lines
- This framework could extend the unification to other types of orthogonal expansions or non-polynomial coefficient equations.
- Connections may exist to broader theories of noncommutative algebra in solving differential equations symbolically.
- Implementation in computer algebra systems could automate recurrence finding for generalized Fourier series.
Load-bearing premise
Every recurrence arising from a linear differential equation with polynomial coefficients can be represented exactly as the numerator of a fraction of recurrence operators without additional conditions on order or leading coefficients.
What would settle it
A counterexample consisting of a linear differential equation with polynomial coefficients and a classical orthogonal polynomial basis where the coefficient recurrence cannot be obtained as the numerator of any fraction of recurrence operators.
Figures
read the original abstract
We consider series expansions in bases of classical orthogonal polynomials. When such a series solves a linear differential equation with polynomial coefficients, its coefficients satisfy a linear recurrence equation. We interpret this equation as the numerator of a fraction of linear recurrence operators. This interpretation lets us give a simple and unified view of previous algorithms computing these recurrences, with a noncommutative Euclidean algorithm as the algorithmic engine. Finally, we demonstrate the effectiveness of our approach on various examples.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that when a series expansion in a basis of classical orthogonal polynomials solves a linear differential equation with polynomial coefficients, the coefficients satisfy a linear recurrence that can be interpreted as the numerator of a fraction of recurrence operators in the noncommutative ring. This algebraic view unifies prior algorithms for computing such recurrences by taking the noncommutative Euclidean algorithm as the central engine, and the approach is demonstrated on various examples.
Significance. If the interpretation applies generally, the work supplies a clean conceptual unification of existing methods for deriving coefficient recurrences in orthogonal-polynomial expansions of DE solutions. Framing the problem via fractions in the noncommutative operator ring and identifying the Euclidean algorithm as the algorithmic engine is a genuine strength; the explicit demonstrations on examples further support practical utility in special-functions computations.
major comments (1)
- [Abstract and algorithmic core] The central claim (stated in the abstract) that every recurrence arising from a linear DE with polynomial coefficients admits an exact representation as the numerator of a recurrence-operator fraction requires an explicit argument that the noncommutative Euclidean algorithm always terminates with the precise numerator; without this, the unification may hold only under unstated restrictions on operator order or leading-coefficient degrees.
minor comments (1)
- [Examples] The examples section would benefit from a concise table listing the input DE, the resulting recurrence, and the operators obtained, to make the unification with prior algorithms immediately visible.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of our work and for the constructive comment on the algorithmic core. We have revised the manuscript to supply the requested explicit argument.
read point-by-point responses
-
Referee: [Abstract and algorithmic core] The central claim (stated in the abstract) that every recurrence arising from a linear DE with polynomial coefficients admits an exact representation as the numerator of a recurrence-operator fraction requires an explicit argument that the noncommutative Euclidean algorithm always terminates with the precise numerator; without this, the unification may hold only under unstated restrictions on operator order or leading-coefficient degrees.
Authors: We agree that an explicit justification is required. The ring of recurrence operators with polynomial coefficients is a noncommutative Euclidean domain (Ore polynomial ring), so the Euclidean algorithm terminates for any input pair and the reduced numerator is exact. In the revised version we have inserted a new subsection (Section 3.2) that recalls the relevant ring axioms, states the division algorithm with remainder degree strictly lower than the divisor, and proves that the algorithm applied to the operator fraction arising from the DE yields precisely the numerator recurrence without restrictions on operator order or leading-coefficient degree. This makes the unification claim fully rigorous for all classical orthogonal polynomial bases. revision: yes
Circularity Check
No significant circularity; algorithmic construction from DE to recurrence is self-contained
full rationale
The paper starts from a linear differential equation with polynomial coefficients, interprets the resulting coefficient recurrence as the numerator of a fraction in the noncommutative ring of recurrence operators, and uses the noncommutative Euclidean algorithm to compute it. This is presented as a direct algorithmic construction that unifies prior methods. No step reduces the claimed result to a fitted parameter, self-definition, or load-bearing self-citation by construction. The derivation chain begins externally from the DE and produces the recurrence without circular reduction. The reader's assessment of score 2.0 aligns with possible minor self-citation but does not indicate load-bearing circularity in the core claim.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math The ring of linear recurrence operators with coefficients in the polynomial ring is a Euclidean domain (or admits a Euclidean algorithm) in the non-commutative sense.
- domain assumption If a series in classical orthogonal polynomials satisfies a linear differential equation with polynomial coefficients, then its coefficient sequence satisfies a linear recurrence relation.
Reference graph
Works this paper leans on
-
[1]
W. A. Al-Salam. The Bessel polynomials.Duke Math. J., 24:529–545, 1957
work page 1957
-
[2]
W. A. Al-Salam and T. S. Chihara. Another characterization of the classical orthogonal polynomials.SIAM J. Math. Anal., 3:65–70, 1972
work page 1972
-
[3]
I. Area, E. Godoy, A. Ronveaux, and A. Zarzo. Corrigendum: “Minimal recurrence rela- tions for connection coefficients between classical orthogonal polynomials: discrete case”.J. Comput. Appl. Math., 89(2):309–325, 1998. 50 A. BENOIT, N. BRISEBARRE, AND B. SAL VY
work page 1998
-
[4]
Askey.Orthogonal polynomials and special functions
R. Askey.Orthogonal polynomials and special functions. Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1975
work page 1975
- [5]
-
[6]
A. Benoit and B. Salvy. Chebyshev expansions for solutions of linear differential equations. In J. May, editor,ISSAC ’09: Proceedings of the twenty-second international symposium on Symbolic and algebraic computation, pages 23–30, 2009
work page 2009
-
[7]
A. Bostan, F. Chyzak, P. Lairez, and B. Salvy. Generalized Hermite reduction, creative telescoping and definite integration of D-finite functions. InISSAC’18—Proceedings of the 2018 ACM International Symposium on Symbolic and Algebraic Computation, pages 95–102. ACM Press, 2018
work page 2018
-
[8]
J. P. Boyd.Chebyshev and Fourier Spectral Methods. Dover Publications Inc., 2nd edition, 2000
work page 2000
-
[9]
S. Chen, M. Kauers, and M. F. Singer. Desingularization of Ore operators.J. Symbolic Comput., 74:617–626, 2016
work page 2016
-
[10]
T. S. Chihara.An introduction to orthogonal polynomials. Mathematics and its Applications, Vol. 13. Gordon and Breach Science Publishers, New York-London-Paris, 1978
work page 1978
-
[11]
C. W. Clenshaw. The numerical solution of linear differential equations in Chebyshev series. Proc. Cambridge Philos. Soc., 53:134–149, 1957
work page 1957
-
[12]
P. M. Cohn.Free ideal rings and localization in general rings, volume 3 ofNew Mathematical Monographs. Cambridge University Press, Cambridge, 2006. [13]NIST Digital Library of Mathematical Functions.https://dlmf.nist.gov/, Release 1.2.5 of 2025-12-15. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Mi...
work page 2006
-
[13]
D. Elliott. The expansion of functions in ultraspherical polynomials.J. Austral. Math. Soc., 1:428–438, 1959/1960
work page 1959
-
[14]
D. Elliott. The evaluation and estimation of the coefficients in the Chebyshev series expansion of a function.Math. Comp., 18:274–284, 1964
work page 1964
-
[15]
A. Erdélyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi.Higher transcendental func- tions. Vol. II. Robert E. Krieger Publishing Co., Inc., Melbourne, FL, 1981. Based on notes left by Harry Bateman, Reprint of the 1953 original
work page 1981
- [16]
- [17]
-
[18]
J. L. Fields and J. Wimp. Expansions of hypergeometric functions in hypergeometric func- tions.Math. Comp., 15:390–395, 1961
work page 1961
- [19]
-
[20]
K.O.Geddes.SymboliccomputationofrecurrenceequationsfortheChebyshevseriessolution of linear ODE’s. In C. M. Andersen, editor,Proceedings of the 1977 MACSYMA User’s Conference, pages 405–423, 1977. NASA CP-2012
work page 1977
- [21]
- [22]
-
[23]
D. Gottlieb and S. A. Orszag.Numerical analysis of spectral methods: theory and applica- tions. Society for Industrial and Applied Mathematics, Philadelphia, Pa., 1977. CBMS-NSF Regional Conference Series in Applied Mathematics, No. 26
work page 1977
-
[24]
J. Govaerts, M. N. Hounkonnou, and A. Z. Msezane, editors.Contemporary problems in mathematical physics. World Scientific Publishing Co., Inc., River Edge, NJ, 2002
work page 2002
-
[25]
E. Hille. Contributions to the theory of Hermitian series.Duke Math. J., 5:875–936, 1939
work page 1939
-
[26]
E. Hille. Contributions to the theory of Hermitian series. II. The representation problem. Trans. Amer. Math. Soc., 47:80–94, 1940. RECURRENCE OPERATORS FOR GENERALIZED FOURIER SERIES 51
work page 1940
-
[27]
E. Hille. Contributions to the theory of Hermitian series. III. Mean values.Internat. J. Math. Math. Sci., 3(3):407–421, 1980
work page 1980
-
[28]
W. Hodges.A shorter model theory. Cambridge University Press, Cambridge, 1997
work page 1997
-
[29]
M. N. Hounkonnou, S. Belmehdi, and A. Ronveaux. Linearization of arbitrary products of classical orthogonal polynomials.Appl. Math. (Warsaw), 27(2):187–196, 2000
work page 2000
-
[30]
W. T. Howell. On products of Laguerre polynomials.The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 24(161):396–405, 1937
work page 1937
-
[31]
E. L. Ince.Ordinary differential equations. Dover Publications, New York, 1956
work page 1956
-
[32]
M. E. H. Ismail.Classical and quantum orthogonal polynomials in one variable, volume 98 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 2005
work page 2005
-
[33]
R. Koekoek, P. A. Lesky, and R. F. Swarttouw.Hypergeometric orthogonal polynomials and theirq-analogues. Springer Monographs in Mathematics. Springer-Verlag, Berlin, 2010
work page 2010
-
[34]
W. Koepf and E. N. Chiadjeu. Algorithmic approach for formal Fourier series.Math. Comput. Sci., 9(3):365–389, 2015
work page 2015
- [35]
-
[36]
Lang.Algebra, volume 211 ofGraduate Texts in Mathematics
S. Lang.Algebra, volume 211 ofGraduate Texts in Mathematics. Springer-Verlag, New York, third edition, 2002
work page 2002
-
[37]
S. Lewanowicz. Construction of a recurrence relation of the lowest order for coefficients of the Gegenbauer series.Zastos. Mat., 15(3):345–396, 1976
work page 1976
-
[38]
S. Lewanowicz. Construction of the lowest-order recurrence relation for the Jacobi coefficients. Zastos. Mat., 17(4):655–675, 1983
work page 1983
-
[39]
S. Lewanowicz. Recurrence relations for the coefficients in Jacobi series solutions of linear differential equations.SIAM J. Math. Anal., 17(5):1037–1052, 1986
work page 1986
-
[40]
S. Lewanowicz. A new approach to the problem of constructing recurrence relations for the Jacobi coefficients.Zastos. Mat., 21(2):303–326, 1991
work page 1991
-
[41]
S. Lewanowicz. Quick construction of recurrence relations for the Jacobi coefficients.J. Com- put. Appl. Math., 43(3):355–372, 1992
work page 1992
-
[42]
S. Lewanowicz. Results on the associated classical orthogonal polynomials.J. Comput. Appl. Math., 65(1-3):215–231, 1995
work page 1995
-
[43]
S. Lewanowicz. Second-order recurrence relation for the linearization coefficients of the clas- sical orthogonal polynomials.J. Comput. Appl. Math., 69(1):159 – 170, 1996
work page 1996
-
[44]
S. Lewanowicz. Recurrences for the coefficients of series expansions with respect to classical orthogonal polynomials.Appl. Math. (Warsaw), 29(1):97–116, 2002
work page 2002
-
[45]
S. Lewanowicz. Construction of recurrences for the coefficients of expansions inq-classical orthogonal polynomials.J. Comput. Appl. Math., 153(1-2):295–309, 2003
work page 2003
-
[46]
S. Lewanowicz, E. Godoy, I. Area, A. Ronveaux, and A. Zarzo. Recurrence relations for the coefficients of the Fourier series expansions with respect toq-classical orthogonal polynomials. Numer. Algorithms, 23(1):31–50, 2000
work page 2000
-
[47]
S. Lewanowicz and P. Woźny. Recurrence relations for the coefficients in series expansions with respect to semi-classical orthogonal polynomials.Numer. Algorithms, 35(1):61–79, 2004
work page 2004
-
[48]
Y. L. Luke. Expansion of the confluent hypergeometric function in series of Bessel functions. Math. Tables Aids Comput., 13:261–271, 1959
work page 1959
-
[49]
Y. L. Luke.The special functions and their approximations, Vol. I. Mathematics in Science and Engineering, Vol. 53. Academic Press, New York, 1969
work page 1969
-
[50]
Y. L. Luke and R. L. Coleman. Expansion of hypergeometric functions in series of other hypergeometric functions.Math. Comp., 15:233–237, 1961
work page 1961
-
[51]
P. Maroni and Z. da Rocha. Connection coefficients between orthogonal polynomials and the canonical sequence: an approach based on symbolic computation.Numer. Algorithms, 47(3):291–314, 2008
work page 2008
-
[52]
P. Maroni and Z. da Rocha. Connection coefficients for orthogonal polynomials: symbolic computations, verifications and demonstrations in the Mathematica language.Numer. Algo- rithms, 63(3):507–520, 2013
work page 2013
-
[53]
J. C. Mason and D. C. Handscomb.Chebyshev Polynomials. Chapman & Hall/CRC, 2003
work page 2003
-
[54]
Mora.Solving polynomial equation systems
T. Mora.Solving polynomial equation systems. Vol. IV. Buchberger theory and beyond, vol- ume 158 ofEncyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 2016. 52 A. BENOIT, N. BRISEBARRE, AND B. SAL VY
work page 2016
-
[55]
A. F. Nikiforov and V. B. Uvarov.Special functions of mathematical physics. Birkhäuser Verlag, Basel, 1988
work page 1988
-
[56]
F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors.NIST Handbook of Mathematical Functions. Cambridge University Press, 2010
work page 2010
-
[57]
S. Olver and A. Townsend. A fast and well-conditioned spectral method.SIAM Rev., 55(3):462–489, 2013
work page 2013
-
[58]
O. Ore. Linear equations in non-commutative fields.Ann. of Math. (2), 32:463–477, 1931
work page 1931
-
[59]
O. Ore. Theory of non-commutative polynomials.Ann. of Math. (2), 34(3):480–508, 1933
work page 1933
-
[60]
Paszkowski.Zastosowania numeryczne wielomianów i szeregów Czebyszewa
S. Paszkowski.Zastosowania numeryczne wielomianów i szeregów Czebyszewa. Państwowe Wydawnictwo Naukowe, Warsaw, 1975. Podstawowe Algorytmy Numeryczne. [Fundamental Numerical Algorithms]
work page 1975
- [61]
-
[62]
A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev.Integrals and Series. Volume 2: Special functions. Gordon and Breach, 1986
work page 1986
-
[63]
A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev.Integrals and Series. Volume 3: More special functions. Gordon and Breach, 1989
work page 1989
-
[64]
E. D. Rainville.Special functions. Chelsea Publishing Co., Bronx, N.Y., first edition, 1971
work page 1971
- [65]
-
[66]
L. Rebillard.Étude théorique et algorithmique des séries de Chebyshev, solutions d’équations différentielles holonomes. PhD thesis, Institut national polytechnique de Grenoble, 1998. https://theses.hal.science/tel-00008571/
work page 1998
-
[67]
L. Rebillard and H. Zakrajšek. Recurrence relations for the coefficients in hypergeometric series expansions. In I. Kotsireas and E. Zima, editors,Computer Algebra 2006. Latest Ad- vances in Symbolic Algorithms, pages 158–180. World Sci. Publ., Hackensack, NJ, 2007
work page 2006
- [68]
-
[69]
A. Ronveaux, M. N. Hounkonnou, and S. Belmehdi. Generalized linearization problems.J. Phys. A, 28(15):4423–4430, 1995
work page 1995
-
[70]
A. Ronveaux, A. Zarzo, I. Area, and E. Godoy. Expansion of hypergeometric polynomials: several approaches. In Govaerts et al. [25], pages 446–459
-
[71]
A. Ronveaux, A. Zarzo, and E. Godoy. Recurrence relations for connection coefficients be- tween two families of orthogonal polynomials.J. Comput. Appl. Math., 62(1):67–73, 1995
work page 1995
-
[72]
B. Salvy and P. Zimmermann. Gfun: a Maple package for the manipulation of generating and holonomic functions in one variable.ACM Trans. Math. Softw., 20(2):163–177, 1994
work page 1994
-
[73]
P. Samuel. About Euclidean rings.J. Algebra, 19:282–301, 1971
work page 1971
-
[74]
R. P. Stanley.Enumerative combinatorics, volume 2. Cambridge University Press, 1999
work page 1999
-
[75]
G. Szegő.Orthogonal polynomials. American Mathematical Society, Providence, R.I., fourth edition, 1975. American Mathematical Society, Colloquium Publications, Vol. XXIII
work page 1975
-
[76]
R. Szwarc. Linearization and connection coefficients of orthogonal polynomials.Monatsh. Math., 113(4):319–329, 1992
work page 1992
-
[77]
H. C. Thacher, Jr. Conversion of a power to a series of Chebyshev polynomials.Comm. ACM, 7(3):181–182, 1964
work page 1964
-
[78]
M. van der Put and M. F. Singer.Galois theory of difference equations, volume 1666 of Lecture Notes in Mathematics. Springer-Verlag, Berlin, 1997
work page 1997
-
[79]
A. Verma. Certain expansions of the basic hypergeometric functions.Math. Comp., 20:151– 157, 1966
work page 1966
-
[80]
Wimp.Computation with Recurrence Relations
J. Wimp.Computation with Recurrence Relations. Pitman, Boston, 1984
work page 1984
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.