On the stability of cosmological solutions in Brans-Dicke theory

We use scaled variables to deduce a set of coupled equations valid for isotropic and anisotropic Bianchi type I, V, IX models in some scalar-tensor theories. In particular, for the Brans-Dicke theory the equations decouple, then one is able to integrate the system for FRW and Bianchi I and V models. Additionally, we investigate, through the method of Lyapunov, the stability of our solutions; particular results are presented for BD theory. Furthermore, an extended Lyapunov theorem is demonstrated in the appendix.