Large Dynamical System Analysis for Anisotropic Universe in Brans-Dicke Theory

In this work, we have studied the Brans-Dicke (BD) cosmology in anisotropic models. We present three dimensional dynamical system describing the evolution of anisotropic models containing perfect fluid and BD scalar field with self-interacting potential. The relevant equations have been transformed into the dynamical system. The critical points and the corresponding eigen values have been found in radiation, dust, dark energy, $\mathrm{\Lambda }$CDM and phantom phases of the universe. The natures and the stability around the critical points have also been investigated.