We use numerical integrations to study the asymptotical behaviour of a homogeneous but anisotropic Bianchi type IX model in General Relativity with a massive scalar field. As it is well known, for a Brans-Dicke theory, the asymptotical behaviour of the metric functions is ruled only by the BransDicke coupling constant ω 0 with respect to the value −3/2. In this paper we examine if such a condition still exists with a massive scalar field. We also show that, contrary to what occurs for a massless scalar field, the singularity oscillatory approach may exist in presence of a massive scalar field having a positive energy density.