use the finite-dimensional subspace spanned by the first N vectors in some orthogonal polynomial basis, such as the Legendre polynomials. Notes References . . . . Curve fitting Numerical differential equations