  We study the problem of the instability of inhomogeneous radiation universes
in quadratic lagrangian theories of gravity written as a system of evolution
equations with constraints. We construct formal series expansions and show that
the resulting solutions have a smaller number of arbitrary functions than that
required in a general solution. These results continue to hold for more general
polynomial extensions of general relativity.
