Extragalactic sources of high-energy cosmic gamma-rays, like blazars and quasars, are usually highly variable and the temporal characteristics of their light curves are key to the understanding of their physics. In order to deal with blazar flares in different parts of electromagnetic spectrum - from gamma-rays and X-rays to visible light - we developed a new approach to correlation study based on Monte Carlo simulations and the Fourier convolution.