A numerical analysis of chaotic behaviour in Bianchi IX models

In this paper we investigate the chaotic behaviour of the Bianchi IX cosmological models using techniques developed in the study of dynamical systems and chaotic behaviour. We numerically calculate the Lyapunov exponent, , and show that instead of converging to a constant value, it decreases steadily. We study this effect further by studying the Lyapunov exponent using short-time averages. We show that the usual method of calculating is invalid in the case of a cosmological model.