Motivated by a connection with the Lovasz Local Lemma, we design an algorithm to compute the multivariate independence polynomial of a graph in the "Shearer region". More generally, we can compute the independent polynomial in a natural complex extension of the Shearer region. A contribution of technical interest is a multivariate version of the correlation decay method.