  The approximation properties of the Aldaz-Kounchev-Render (AKR) operators are
discussed and classes of functions for which these operators approximate better
than the classical Bernstein operators are described. The new results are then
extended to the bivariate case on the square $[0,1]^2$ and compared with other
existing results known in literature. Several numerical examples, illustrating
the relevance and supporting the theoretical findings, are presented
