Mixed Virtual Element Methods for general second order elliptic problems on polygonal meshes

In the present paper we introduce a Virtual Element Method (VEM) for the approximate solution of general linear second order elliptic problems in mixed form, allowing for variable coefficients. We derive a theoretical convergence analysis of the method and develop a set of numerical tests on a benchmark problem with known solution.