##### Abstract

In this paper, a rigorous method to compute Floquet normal forms of fundamental matrix solutions of nonautonomous linear differential equations with periodic coefficients is introduced. The Floquet normal form of a fundamental matrix solution F(t) is a canonical decomposition of the form F(t) = Q(t)eRt, where Q(t) is a real periodic matrix and R is a constant matrix. To rigorously compute the Floquet normal form, the idea is to use the regularity of Q(t) and to simultaneously solve for R and Q(t) with the contraction mapping theorem in a Banach space of rapidly decaying coefficients. The explicit knowledge of R and Q can then be used to construct, in a rigorous computer-assisted way, stable and unstable bundles of periodic orbits of vector fields. The new proposed method does not require rigorous numerical integration of the ODE.