Miwa deformation of WLZZ matrix models

In this paper we propose a universal way to define multivariate orthogonal polynomials for the positive branch of WLZZ matrix models. This allows us to generalize from the Gaussian case, previously considered in the literature, the strict superintegrability property and the Miwa deformation formulas.