Abstract: We describe an efficient numerical approach to calculate the longitudinal and transverse Kubo conductivities of large systems using Bastin's formulation. We expand the Green's functions in terms of Chebyshev polynomials and compute the conductivity tensor for any temperature and chemical potential in a single step. To illustrate the power and generality of the approach, we calculate the conductivity tensor for the quantum Hall effect in disordered graphene and analyze the effect of the disorder in a Chern insulator in Haldane's model on a honeycomb lattice.