We introduce a new family of symmetric multivariate polynomials, whose coefficients are meromorphic functions of two parameters $(q,t)$ and polynomial in a further two parameters $(u,v)$. We evaluate these polynomials explicitly as a matrix product. At $u=v=0$ they reduce to Macdonald polynomials, while at $q=0$, $u=v=s$ they recover a family of inhomogeneous symmetric functions originally introduced by Borodin.