Rigorous Numerics for Nonlinear Differential Equations Using Chebyshev Series

A computational method based on Chebyshev series to rigorously compute solutions of initial and boundary value problems of analytic nonlinear vector fields is proposed. The idea is to recast solutions as fixed points of an operator defined on a Banach space of rapidly decaying Chebyshev coefficients and to use the so-called radii polynomials to show the existence of a unique fixed point nearby an approximate solution. As applications, solutions of initial value problems in the Lorenz equations and symmetric connecting orbits in the Gray-Scott equation are rigorously computed. The symmetric connecting orbits are obtained by solving a boundary value problem with one of the boundary values in the stable manifold.