  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, $\Lambda$CDM and phantom
phases of the universe. The natures and the stability around the critical
points have also been investigated.
