AbstractA simple dynamical model for a BL Lacertae active galaxy is presented. The model consists of a logarithmic potential with an additional term representing internal perturbations. The time independent and the evolving model are investigated. In both cases we search for regular and chaotic motion and study the velocity distribution near the centre of the system. Numerical calculations suggest that responsible for the chaotic phenomena is the internal perturbation, the flattening parameter and the dense nucleus. The radius of the nucleus also affects the maximum velocity in the central regions of the galaxy. Our numerical outcomes are supported by theoretical arguments and analytical calculations. A linking of our numerical outcomes to observational data is also presented. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)