In this paper two types of tensor product finite macro-elements are contrasted, the former being the well known Lagrange type and the latter the Bézier (Bernstein) type elements. Although they have a different mathematical origin and seemingly are irrelevant, they both are based on complete polynomials thus sharing the same functional space, i.e. the classes { x n } and { y n } . Therefore, from the theoretical point of view it is anticipated that they should lead to numerically identical results in both static and dynamic analysis. For both types of elements details are provided concerning the main computer programming steps, while selective parts of a typical MATLAB® code are presented. Numerical application includes static (Laplace, Poisson), eigenvalue (acoustics) and transient (heat conduction) problems of rectangular, circular and elliptic shapes, which were treated as a single macroelement. In agreement to the theory, in all six examples the results obtained using Bézier and Lagrange polynomials were found to be identical and of exceptional accuracy. 
Bézier versus Lagrange polynomials-based finite element analysis of 2-D potential problems