Our subject is a set of equations
where the uj(n)(j = 1, 2, …, k) are k “unknown” functions of the integer variable n, the zi(n) (i = 1, 2, … h) are h “known” functions of n, and the Aij(n) are hk “known” operators
which are polynomials in E, each of fixed order pij but with coefficients which may vary with n. E is the usual operator defined by
Our first task is to determine whether the equations (1) are self-consistent. Secondly, if they are self-consistent, we ask what follows from them for a given subset of the unknowns, e.g. for (uj+1, …, uk) in other words we wish to eliminate (u1, …, uj). In particular we wish to eliminate all the variables but one, say uk. We shall in fact find that either uk is arbitrary or else that it has only to satisfy a single linear recurrence relation : and the order of that relation is of interest to us. Thirdly, we ask that reduction to standard form is possible, with or without a transformation of the unknowns themselves.