Bounds for Calderón-Zygmund operators with matrix A2 weights

Our method uses a Bellman function technique introduced by S. Treil to obtain the right estimates for the norm of dyadic Haar shift operators. We then apply the representation theorem of T. Hytönen to extend the result to general Calderón-Zygmund operators.