On the bivariate Padovan polynomials matrix