Approximation of smooth functions using Bernstein polynomials in multiple variables