A numerical code to study the variability of Blazar emission

We present a numerical code, written in C, which can be used to simulate or to analyze the emission of Blazars over the entire electromagnetic spectrum. Our code can reproduce the following features: synchrotron emission, inverse Compton emission (Thomson Klein-Nishina regime) external Compton emission, accretion disk variability using a Cellular Automata algorithm, temporal evolution of the emitting plasma energy distribution, flaring phenomena, light curves in the rest and in the observer frame (taking account for time crossing effects). In this paper we will show mainly the accretion disk simulation, and the implications in the External Compton scenario.