We consider the problem of representing a solution to the Cauchy problem for an ordinary differential equation as a Fourier series in polynomials
l r,k α (
x) (
k = 0, 1,...) that are Sobolev-orthonormal with respect to the inner product
$$leftlangle {f,g} rightrangle = sumlimits_{v = 0}^{r - 1} {{f^{(v)}}(0){g^{(v)}}} (0) + intlimits_0^infty {{f^{(r)}}(t)} {g^{(r)}}(t){t^alpha }{e^{ - t}}dt$$
, and generated by the classical orthogonal Laguerre polynomials
L k α (
x) (
k = 0, 1,...). The polynomials
l r,k α (
x) are represented as expressions containing the Laguerre polynomials
L n α?r (
x). An explicit form of the polynomials
l r,k+r α (
x) is established as an expansion in the powers
x r+l ,
l = 0,...,
k. These results can be used to study the asymptotic properties of the polynomials
l r,k α (
x) as
k→∞and the approximation properties of the partial sums of Fourier series in these polynomials.