In mathematics a P-recursive equation can be solved for polynomial solutions.
Sergei A. Abramov in 1989 and Marko Petkovšek in 1992 described an algorithm which finds all polynomial solutions of those recurrence equations with polynomial coefficients.
[1][2] The algorithm computes a degree bound for the solution in a first step.
In a second step an ansatz for a polynomial of this degree is used and the unknown coefficients are computed by a system of linear equations.
This article describes this algorithm.
In 1995 Abramov, Bronstein and Petkovšek showed that the polynomial case can be solved more efficiently by considering power series solution of the recurrence equation in a specific power basis (i.e. not the ordinary basis
[3] Other algorithms which compute rational or hypergeometric solutions of a linear recurrence equation with polynomial coefficients also use algorithms which compute polynomial solutions.
be a field of characteristic zero and
a recurrence equation of order
, polynomial right-hand side
and unknown polynomial sequence
denotes the degree of a polynomial
denotes the leading coefficient of the polynomial.
denotes the falling factorial and
the set of nonnegative integers.
This is called a degree bound for the polynomial solution
This bound was shown by Abramov and Petkovšek.
[1][2][3][4] The algorithm consists of two steps.
In a first step the degree bound is computed.
In a second step an ansatz with a polynomial
of that degree with arbitrary coefficients in
is made and plugged into the recurrence equation.
Then the different powers are compared and a system of linear equations for the coefficients of
This is called the method undetermined coefficients.
[5] The algorithm returns the general polynomial solution of a recurrence equation.
Applying the formula for the degree bound on the recurrence equation
Hence one can use an ansatz with a quadratic polynomial
Plugging this ansatz into the original recurrence equation leads to
This is equivalent to the following system of linear equations
Therefore the only polynomial solution is