On the numerical computation of bivariate Lagrange polynomials