  The two-matrix model can be solved by introducing bi-orthogonal polynomials.
In the case the potentials in the measure are polynomials, finite sequences of
bi-orthogonal polynomials (called
  "windows") satisfy polynomial ODEs as well as deformation equations (PDEs)
and finite difference equations (Delta-E) which are all Frobenius compatible
and define discrete and continuous isomonodromic deformations for the irregular
ODE, as shown in previous works of ours.
  In the one matrix model an explicit and concise expression for the
coefficients of these systems is known and it allows to relate the partition
function with the isomonodromic tau-function of the overdetermined system.
Here, we provide the generalization of those expressions to the case of
bi-orthogonal polynomials, which enables us to compute the determinant of the
fundamental solution of the overdetermined system of ODE+PDEs+Delta-E.
