Next Article in Journal
The Lambert-F Distributions Class: An Alternative Family for Positive Data Analysis
Next Article in Special Issue
Multiple Meixner Polynomials on a Non-Uniform Lattice
Previous Article in Journal
Applying the Modified Delphi Method and DANP to Determine the Critical Selection Criteria for Local Middle and Top Management in Multinational Enterprises
Previous Article in Special Issue
On the Finite Orthogonality of q-Pseudo-Jacobi Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Multivariate Bernstein Polynomials

by
Mama Foupouagnigni
1,2,† and
Merlin Mouafo Wouodjié
3,*,†
1
African Institute for Mathematical Sciences, P.O. Box 608, Limbe, Cameroon
2
Higher Teacher Training College, University of Yaounde I, P.O. Box 47, Yaounde, Cameroon
3
Institute of Mathematics, University of Kassel, Heinrich-Plett Str. 40, 34132 Kassel, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(9), 1397; https://doi.org/10.3390/math8091397
Submission received: 20 June 2020 / Revised: 27 July 2020 / Accepted: 27 July 2020 / Published: 21 August 2020

Abstract

:
In this paper, we first revisit and bring together as a sort of survey, properties of Bernstein polynomials of one variable. Secondly, we extend the results from one variable to several ones, namely—uniform convergence, uniform convergence of the derivatives, order of convergence, monotonicity, fixed sign for the p-th derivative, and deduction of the upper and lower bounds of Bernstein polynomials from those of the corresponding functions.

1. Introduction

It is not often easy to find exact solutions to most equations governing real-life problems. Even when such solutions can be found, for their practical use, there is a need to find their approximation by polynomials, when possible.
Weierstrass’s theorem asserts the possibility of uniform approximation by polynomials of continuous functions over a closed interval. An analytic function can be expanded in a uniformly convergent power series and a continuous but non-analytic function can be expanded in a uniformly convergent series of general polynomials, with no possibility of rearranging its terms so as to produce a convergent power series [1]. There are many proofs of the Weierstrass theorem. In this paper, we shall present one called “S. Bernstein’s proof” since it also gives an explicit representation of the polynomials for uniform approximation. An example of Bernstein polynomials of one variable is (see Reference [1]):
B n ( f ; x ) = k = 0 n f k n n k x k ( 1 x ) n k .
The polynomials B n ( f ; . ) are known to approximate uniformly their corresponding functions f in [ 0 , 1 ] , provided that f is continuous on the interval [ 0 , 1 ] . Moreover, the n th derivative of B n ( f ; . ) converges uniformly to f ( n ) , provided that f C n [ 0 , 1 ] ( f C n [ 0 , 1 ] means f is continuous on the interval [ 0 , 1 ] and all its derivatives f ( i ) , i = 1 , , n , exist and are also continuous on [ 0 , 1 ] ). According to the literature, there are several ways (with different notations) to define explicitly the Bernstein polynomials of one variable, which lead to the same results as in Reference [1]—uniform convergence as well as uniform convergence of the derivatives (see References [2,3,4,5,6,7]).
Bernstein polynomials have been used for approximation of functions in many areas of Mathematics and other fields such as smoothing in Statistics, and constructing and generalizing Bézier curves (see References [8,9]). We also have some applications:
  • in the quadrature rule on equally spaced knots (see Reference [10]);
  • to numerical differentiation and quadrature and perform numerical experiments showing the effectiveness of the considered technique (see Reference [11]): they perform some numerical tests for approximation by polynomials of a given function, numerical differentiation and numerical quadrature.
In contrast to other modes of approximation, in particular to Tchebyschev or best uniform approximation, the Bernstein polynomials yield smooth approximants; but there is a price that must be paid for these beautiful approximations properties: the convergence of the Bernstein polynomials is very slow.
P. J. Davis [1].
Several authors like A. Jafarian, S. Measoomy Nia, Alireza K. Golmankhaneh and D. Baleanu [12], D. H. Tucker [13], E. H. Doha, A. H. Bhrawy, M. A. Saker [14], George M. Phillips [15], Henryk Gzyl and José Luis Palacios [16], M. H. T Alshbool, I. Hashim [17], P. J. Davis [1], Richard V. Kadison [18] and others, have worked on Bernstein polynomials for functions with one variable, and have derived some beautiful properties whose most important ones are: the uniform convergence, the uniform convergence of derivatives, the fixed sign for the p th derivative, the deduction of upper and lower bounds of Bernstein polynomials from those of the corresponding functions.
The main problem here is to extend these properties of Bernstein polynomials of one variable by generalizing them for functions of several real variables. There are more than five ways with different notations to define Bernstein polynomials in the multi-dimensional case for a function f defined on [ 0 , 1 ] m , where m is the dimension of the space. Some used notations are x = ( x 1 , , x m ) [ 0 , 1 ] m , k = ( k 1 , , k m ) N m , | x | = x 1 + + x m , x k = x 1 k 1 x m k m , | k | = k 1 + + k m , k ! = k 1 ! k m ! and k n = k 1 n , , k m n , n N { 0 } . These notations are used in the following definitions:
(a)
the Bernstein polynomials on ( m 1 ) -dimensional simplex in cube [ 0 , 1 ] m (see References [19,20,21,22]):
b n [ f ] ( x ) : = | k | = n f k n n ! k ! x k , where x 1 + + x m = 1 ;
(b)
the Bernstein polynomials on a m-dimensional simplex and very useful in stochastic analysis (see References [7,23,24,25,26]):
B n , m ( f , x ) : = | k | n f k n n ! k ! ( n | k | ) ! x k ( 1 | x | ) n | k | , where x 1 + + x m 1 ;
(c)
the Bernstein polynomials on simplices with a pointwise orthogonality relation (see Reference [27]):
B n ( f , x ) : = | α | = n f α ¯ n ( | α | ) ! α ! x α ¯ ( 1 | x | ) α m + 1 ,
where α = ( α 1 , , α m + 1 ) N m + 1 and α ¯ = ( α 1 , , α m ) ;
(d)
the Bernstein polynomials on ( m 1 ) -dimensional simplex in cube [ 0 , 1 ] m used in the field of mathematical biology with quadrartic mappings of “standard” ( m 1 ) -dimensional simplex to self (see Reference [21]):
B n ( f ) ( x 1 , , x m ) = k 1 = 0 n k m = 0 n f k n i = 1 m n k i x i k i ( 1 x i ) n k i ;
(e)
the Bernstein polynomials related to the convergence in the uniform and L p -norms for polynomials with integral coefficients (see Reference [28] for two variables, References [1,29] for the generalized one):
B n 1 , , n m ( f ; x 1 , , x m ) = k 1 = 0 n 1 k m = 0 n m f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i ,
where n 1 , , n m N .
Using the above mentioned definitions, many authors did some investigations aimed at generalizing some properties of univariate Bernstein polynomials to Bernstein polynomials with two or more variables. Let us take for example the forms given by (2), (3) and (6).
  • For the multivariate Bernstein polynomials given by Equation ( 2 ) : they have proven the pointwise and uniform convergence (see References [21,22,30]). They have also presented two variations which are preserved and do coincide with classical convexity in the univariate case (see Reference [31]). Moreover, they showed the monotonicity of successive Bernstein polynomials for classically convex functions (see References [31,32]).
  • For the multivariate Bernstein polynomials given by Equation ( 3 ) : they have proven two identities for these polynomials on simplex, which are considered on a pointwise (see References [23,24,25,25]). They also showed that these polynomials converge to their corresponding functions f along with partial derivatives provided that the latter derivatives exist and are continuous (see References [7,26,33]).
  • For the Bernstein polynomials given by Equation ( 6 ) : they established (but only for two variables) the pointwise and uniform convergence; the same result was established for the partial derivatives of the corresponding function (see References [28,29,34,35,36]).
According to the last definition (see Equation ( 6 ) ) of the generalized Bernstein polynomial B n 1 , , n m ( f ; x 1 , , x m ) , there is just a few new results (pointwise and uniform convergence) given by some authors like Clemens Heitzinger (see Reference [29]). That is why we focus, in this paper, on the multivariate Bernstein polynomials given by the form defined by Equation (6).
Many authors have made tremendous efforts to prove some properties for multivariate Bernstein polynomials on simplex (see Reference [23]), which are considered on a pointwise convergence (good approximations of Bernstein polynomials for every continuous function on simplex and the higher dimensional q-analogue of Bernstein polynomials on simplex). Some have tried to study the approximation by combination of Bernstein polynomials or operators (the rate of convergence of derivatives and the smoothness of the derivatives of functions) in order to approximate, under definite conditions, the function more closely than the Bernstein polynomials (see References [3,6,37,38,39,40,41]). Moreover, they defined and also tried to study the q-Bernstein polynomials (see References [2,4,42,43,44,45]), the generalized Bernstein polynomials (see References [10,26]) and the modified Bernstein polynomials (see Reference [25]).
To achieve our aim which is to generalize these beautiful properties for functions with one variable to functions with several variables, we will use the following properties of univariate Bernstein polynomials (with some references where they can be found):
-
Property 1 (see Reference [1]): if f is bounded on [ 0 , 1 ] and continuous at x 0 [ 0 , 1 ] , then
lim n B n ( f ; x 0 ) = f ( x 0 ) ;
-
Property 2 (see References [1,3,5,6,7,10,11,39,46,47,48]): if f C [ 0 , 1 ] , the Bernstein polynomial B n ( f ; · ) converges uniformly to f;
-
Property 3 (see References [1,7,11]): if f C p [ 0 , 1 ] (f is continuous on the interval [ 0 , 1 ] and all its derivatives f ( i ) , i = 1 , , p , exist and are also continuous on [ 0 , 1 ] ), then B n p ( f ; · ) converges uniformly to f ( p ) ;
-
Property 4 (see References [1,21]): if f C p [ 0 , 1 ] , for some 0 p n and A f ( p ) ( x ) B ( A a n d B a r e c o n s t a n t s ) , 0 x 1 , then
A n p n ( n 1 ) ( n p + 1 ) B n ( p ) ( f ; x ) B , 0 x 1 ;
(the consequence is: if f is decreasing (resp. nondecreasing) on [ 0 , 1 ] , then B n ( f ; · ) is decreasing (resp. nondecreasing) on [ 0 , 1 ] .)
-
Property 5 (see References [1,31]): if f is convex in [ 0 , 1 ] , then B n ( f ; · ) is convex there and for n = 2 , 3 ,
B n 1 ( f ; x ) B n ( f ; x ) , 0 < x < 1 ;
-
Property 6 (see References [1,11,39]): if f ( x ) is bounded on [ 0 , 1 ] and x 0 is a point of [ 0 , 1 ] at which f exists and is continuous, then
lim n n B n ( f ; x 0 ) f ( x 0 ) = 1 2 x 0 ( 1 x 0 ) f ( x 0 ) .
The following intermediate results for the case of univariate Bernstein polynomials are needed to prove all those properties:
(a)
see References [1,15]: for x [ 0 , 1 ] we have B n ( 1 ; x ) = 1 ;
(b)
see Reference [1]: for x [ 0 , 1 ] and α R we have
B n ( exp ( α x ) ; x ) = x + exp ( α / n ) + ( 1 x ) n ;
(c)
see References [1,15]: the Bernstein polynomial may be expressed in the form
B n ( f ; x ) = k = 0 n n k k f n ( 0 ) x k ,
where f n ( x ) = f ( x n ) and △ is the forward difference operator with step size h = 1 n defined as follows (see Reference [15]):
f ( x i ) : = f ( x i + 1 ) f ( x i ) , k + 1 f ( x i ) : = k f ( x i ) = k f ( x i + 1 ) k f ( x i ) for k 1 , 1 f ( x i ) : = f ( x i ) and 0 f ( x i ) : = f ( x i ) ,
where k is called the k th difference;
(d)
See References [1,15,21,22]: for a given δ > 0 , n N { 0 } and x [ 0 , 1 ] , we have
k n x δ n k x k ( 1 x ) n k 1 4 n δ 2 ;
(e)
see References [1,15]: for any integer k 0 , the k th derivative of B n + k ( f ; x ) may be expressed in terms of k-th differences of f as
B n + k ( k ) ( f ; x ) = ( n + k ) ! n ! r = 0 n k f r n + k n r x r ( 1 x ) n r ,
for all n 0 , where △ is applied with step size h = 1 n + k ;
(f)
see (Reference [1] Lemma 6.3.5): there is a constant C independent of a given positive integer n such that for all x in [ 0 , 1 ]
k n x n 1 / 4 n k x k ( 1 x ) n k C n 3 / 2 .
With all those properties, we will prove the following for our generalized Bernstein polynomials:
-
Property 1: if f is bounded on [ 0 , 1 ] m and continuous on ( a 1 , , a m ) [ 0 , 1 ] m , then the Bernstein polynomial B n 1 , , n m ( f ; a 1 , , a m ) converges to f ( a 1 , , a m ) , as n i for all i;
-
Property 2: if f C [ 0 , 1 ] m , then the Bernstein polynomial B n 1 , , n m ( f ; · ) converges uniformly to f;
-
Property 3: if f C p 1 p m [ 0 , 1 ] m , then B n 1 , , n m p 1 , , p m ( f ; · ) converges uniformly to f ( p 1 ) , , ( p m ) , where f C p 1 , , p m [ 0 , 1 ] m means f is continuous on [ 0 , 1 ] m and p i f x p i is continuous on [ 0 , 1 ] m , 0 p i n i , i = 1 , , m ;
-
Property 4: if f C p 1 p m [ 0 , 1 ] m , 0 p i n i , i = 1 , , m , and A f ( p 1 ) ( p m ) ( x 1 , , x m ) B , ( x 1 , , x m ) [ 0 , 1 ] m , then
A i = 1 m n i p i n i ( n i 1 ) ( n i p i + 1 ) B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) B
for all ( x 1 , , x m ) [ 0 , 1 ] m . As a consequence, if f ( x 1 , , x m ) is decreasing (resp. nondecreasing) on [ 0 , 1 ] with respect to the variable x i , i 1 , , m , then B n 1 , , n m ( f ; x 1 , , x m ) is decreasing (resp. nondecreasing) on [ 0 , 1 ] with respect to the variable x i there.
-
Property 5: if f ( x 1 , , x m ) is convex in [ 0 , 1 ] m , then, for n i = 2 , 3 , with i = 1 , , m , we have
B n 1 1 , , n m 1 ( f ; x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m )
for all ( x 1 , , x m ) [ 0 , 1 ] m ;
-
Property 6: if f ( x 1 , , x m ) is bounded in [ 0 , 1 ] m and let ( a 1 , , a m ) be a point of [ 0 , 1 ] m at which f ( 2 ) ( 2 ) ( x 1 , , x m ) exists and is continuous, then, for n 1 = n 2 = = n m = n
lim n n B n , , n ( f ; a 1 , , a m ) f ( a 1 , , a m ) = 1 2 i = 1 m a i ( 1 a i ) 2 f x i 2 ( a 1 , , a m ) .
The following intermediate results for the multivariate case, that we will also prove, are needed in order to prove all these cited properties:
(a)
for m 1 , i { 1 , , m } , j i and n i integers, n i 1 , and x i [ 0 , 1 ]
x 1 j 1 x m j m = k 1 = 0 j 1 k m = 0 j m i = 1 m σ j i , k i n i k i x i k i ( 1 x i ) n i k i , B n 1 , , n m 1 ; x 1 , , x m = 1 , B n 1 , , n m i = 1 m x i j i ; x 1 , , x m = i = 1 m B n i x i j i ; x i ,
where σ j i , k i are real numbers;
(b)
for i { 1 , , m } , x i [ 0 , 1 ] and α i R
B n 1 , , n m exp i = 1 m α i x i ; x 1 , , x m = i = 1 m B n i exp ( α i x i ) ; x i = i = 1 m x i exp ( α i / n i ) + ( 1 x i ) n i ,
k 1 = 0 n 1 k m = 0 n m i = 1 m n i k i k i n i x i j i x i k i ( 1 x i ) n i k i = i = 1 m k i = 0 n i n i k i k i n i x i j i x i k i ( 1 x i ) n i k i ;
(c)
the Bernstein polynomial may be expressed in the form
B n 1 , , n m ( f ; x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m 1 t 1 m t m f n 1 , , n m ( 0 , , 0 ) i = 1 m n i t i x i t i ,
where f n 1 , , n m ( t 1 , , t m ) = f t 1 n 1 , , t m n m and i t i , i { 1 , , m } , is the t i th difference(s) related to the i t h variable x i , and applied with step size h i = 1 n i .
(d)
for a given δ i > 0 , i = 1 , , m and ( x 1 , , x m ) [ 0 , 1 ] m , we have
k 1 n 1 x 1 δ 1 k m n m x m δ m i = 1 m n i k i x i k i ( 1 x i ) n i k i 1 4 m i = 1 m n i δ i 2 ,
summing over those values of k i = 0 , , n i , i = 1 m , for which k i n i x i δ i ;
(e)
for p 1 , , p m 0 integers with 0 p i n i , i = 1 , 2 , , m , the p 1 , , p m th derivative of B n 1 + p 1 , , n m + p m ( f ; x 1 , , x m ) may be expressed in terms of p 1 , , p m th differences of f as
B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) = i = 1 m ( n i + p i ) ! n i ! t 1 = 0 n 1 t m = 0 n m 1 p 1 m p m f t 1 n 1 + p 1 , , t m n m + p m i = 1 m n i t i x i t i ( 1 x i ) n i t i ,
where the Δ k is acting with the mesh h k = 1 n k + p k ;
(f)
there is a constant C independent of given positive integers n i such that for all ( x 1 , , x m ) [ 0 , 1 ] m , i = 1 , , m ,
k 1 n 1 x 1 n 1 1 / 4 k m n m x m n m 1 / 4 i = 1 m n i k i x i k i ( 1 x i ) n i k i C i = 1 m n i 3 / 2 .

2. Definition and Properties

We shall now derive the properties of the m-dimensional generalization of the Bernstein polynomials. First of all, we give m-dimensional generalization of the Bernstein polynomials.
Definition 1
(Generalization). (P. J. Davis [1], page 122)
Let f ( x 1 , , x m ) be defined on [ 0 , 1 ] m and n i 1 , i = 1 , , m , integers. The ( n 1 , , n m ) t h Bernstein polynomial for f ( x 1 , , x m ) is given by
B n 1 , , n m ( f ; x 1 , , x m ) = k 1 = 0 n 1 k m = 0 n m f k 1 n 1 , , k m n m i = 1 m b k i , n i ( x i ) ,
where b k i , n i ( x i ) = n i k i x i k i ( 1 x i ) n i k i , i = 1 , , m .
Notice that
B n 1 , , n m ( f ; 0 , , 0 ) = f ( 0 , , 0 ) ; B n 1 , , n m ( f ; 1 , , 1 ) = f ( 1 , , 1 ) .
We call B n 1 , , n m the Bernstein operator; it maps a function f, defined on [ 0 , 1 ] m , to B n 1 , , n m f , where the function B n 1 , , n m f evaluated at ( x 1 , , x m ) is denoted by B n 1 , , n m ( f ; x 1 , , x m ) . This operator is linear, since it follows from Equation ( 27 ) that
B n 1 , , n m ( λ f + μ g ) = λ B n 1 , , n m f + μ B n 1 , , n m g ,
for all functions f and g defined on [ 0 , 1 ] m , and all real λ and μ .
Using the results obtained with one variable, it is easy to obtain some pertinent properties for the generalization of Bernstein polynomials.
Lemma 1.
1 
Let m 1 , j 1 , , j m be integers. Then,
x 1 j 1 x m j m = k 1 = 0 j 1 k m = 0 j m i = 1 m σ j i , k i b k i , j i ( x i ) ,
where σ j i , k i , k i = 0 , , j i , 1 i m , are real numbers.
2 
For i { 1 , , m } , let n i 1 , j i be integer and x i [ 0 , 1 ] . Then,
B n 1 , , n m 1 ; x 1 , , x m = 1 , B n 1 , , n m i = 1 m x i j i ; x 1 , , x m = i = 1 m B n i x i j i ; x i ,
(a similar result for one variable can be found in References [1,15] for the first equation of the system ( 31 ) ).
3 
For i { 1 , , m } , let x i [ 0 , 1 ] and α i R . Then,
B n 1 , , n m exp i = 1 m α i x i ; x 1 , , x m = i = 1 m B n i exp ( α i x i ) ; x i = i = 1 m x i exp ( α i / n i ) + ( 1 x i ) n i
(a similar result for one variable can be found in Reference [1]).
4 
For x i [ 0 , 1 ] , 1 i m , we have
k 1 = 0 n 1 k m = 0 n m i = 1 m n i k i k i n i x i j i x i k i ( 1 x i ) n i k i = i = 1 m k i = 0 n i n i k i k i n i x i j i x i k i ( 1 x i ) n i k i .
Proof of Lemma 1.
See References [1,15] for the proofs. □

3. Approximation by Bernstein Polynomials with Several Variables

We express in this theorem our multivariate Bernstein polynomials in terms of differences. A similar result exists for one variable (see References [1,15]).
Theorem 1.
Let f ( x 1 , , x m ) be defined on [ 0 , 1 ] m . We have for positive integers n i , 1 i m ,
B n 1 , , n m ( f ; x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m 1 t 1 m t m f n 1 , , n m ( 0 , , 0 ) i = 1 m n i t i x i t i ,
where f n 1 , , n m ( t 1 , , t m ) = f t 1 n 1 , , t m n m and i t i , i { 1 , , m } , is the t i t h difference(s) (see References [1,15]) related to the i t h variable x i , and applied with step size h i = 1 n i .
Proof of Theorem 1.
We will sometimes use the square brackets to draw attention on computations of certain terms.
B n 1 , , n m ( f ; x 1 , , x m )
= k 1 = 0 n 1 k m = 0 n m f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i
= k 2 = 0 n 2 k m = 0 n m k 1 = 0 n 1 f k 1 n 1 , k 2 n 2 , , k m n m n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i .
We take the first variable of f and fix the others. Hence f k 1 n 1 , , k m n m is seen as a function g 1 k 1 n 1 , and we apply the following result from P. J. Davis [1], about Bernstein polynomials with one variable:
B n ( g ; x ) = t = 0 n t g n ( 0 ) n t x t ,
where n N and g n t = g t n . Then we get
B n 1 , , n m ( f ; x 1 , , x m ) = k 2 = 0 n 2 k m = 0 n m t 1 = 0 n 1 1 t 1 g 1 ( 0 ) x 1 t 1 n 1 t 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 2 = 0 n 2 k m = 0 n m t 1 = 0 n 1 1 t 1 f 0 , k 2 n 2 , , k m n m x 1 t 1 n 1 t 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 3 = 0 n 3 k m = 0 n m t 1 = 0 n 1 k 2 = 0 n 2 1 t 1 f 0 , k 2 n 2 , k 3 n 3 , , k m n m n 2 k 2 x 2 k 2 ( 1 x 2 ) n 2 k 2 x 1 t 1 n 1 t 1 × i = 3 m n i k i x i k i ( 1 x i ) n i k i .
We repeat the action with the second variable of the function 1 t 1 f 0 , k 2 n 2 , , k m n m , to get
B n 1 , , n m ( f ; x 1 , , x m ) = k 3 = 0 n 3 k m = 0 n m t 1 = 0 n 1 t 2 = 0 n 2 1 t 1 2 t 2 f 0 , 0 , k 3 n 3 , , k m n m i = 1 2 n i t i x i t i i = 3 m n i k i x i k i ( 1 x i ) n i k i .
Then B n 1 , , n m ( f ; x 1 , , x m )
= k 4 = 0 n 4 k m = 0 n m t 1 = 0 n 1 t 2 = 0 n 2 k 3 = 0 n 3 1 t 1 2 t 2 f 0 , 0 , k 3 n 3 , , k m n m n 3 k 3 x 3 k 3 ( 1 x 3 ) n 3 k 3 i = 1 2 n i t i x i t i × i = 4 m n i k i x i k i ( 1 x i ) n i k i .
Similarly, with the third variable of 1 t 1 2 t 2 f 0 , 0 , k 3 n 3 , , k m n m , we get
B n 1 , , n m ( f ; x 1 , , x m ) = k 4 = 0 n 4 k m = 0 n m t 1 = 0 n 1 t 2 = 0 n 2 t 3 = 0 n 3 1 t 1 2 t 2 3 t 3 f 0 , 0 , 0 , k 4 n 4 , , k m n m i = 1 3 n i t i x i t i i = 4 m n i k i x i k i ( 1 x i ) n i k i .
By repeating, successively, the actions with the 4 th , 5 th , , m th variables but with the functions 1 t 1 2 t 2 3 t 3 f 0 , 0 , 0 , k 4 n 4 , , k m n m , 1 t 1 2 t 2 3 t 3 4 t 4 f 0 , 0 , 0 , 0 , k 5 n 5 , , k m n m , , 1 t 1 2 t 2 m 1 t m 1 f 0 , 0 , 0 , , 0 , k m n m respectively, we obtain the result. □
In order to show the pointwise and uniform convergence of our multivariate Bernstein polynomials, let us state and prove the following lemma. A similar result exists for one variable (see References [1,15,21,22]).
Lemma 2.
For a given δ i > 0 , i = 1 , , m and ( x 1 , , x m ) [ 0 , 1 ] m , we have
k 1 n 1 x 1 δ 1 k m n m x m δ m i = 1 m n i k i x i k i ( 1 x i ) n i k i 1 4 m i = 1 m n i δ i 2 ,
summing over those values of k i = 0 , , n i for which k i n i x i δ i .
Proof of Lemma 2. 
We will sometimes use the square brackets to draw attention to computations of certain terms.
k 1 n 1 x 1 δ 1 k m n m x m δ m i = 1 m n i k i x i k i ( 1 x i ) n i k i = k 2 n 2 x 2 δ 2 k m n m x m δ m k 1 n 1 x 1 δ 1 n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i = k 1 n 1 x 1 δ 1 n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 k 3 n 3 x 3 δ 3
k m n m x m δ m k 2 n 2 x 2 δ 2 n 2 k 2 x 2 k 2 ( 1 x 2 ) n 2 k 2 i = 3 m n i k i x i k i ( 1 x i ) n i k i = i = 1 2 k i n i x i δ i n i k i x i k i ( 1 x i ) n i k i k 3 n 3 x 3 δ 3 k m n m x m δ m i = 3 m n i k i x i k i ( 1 x i ) n i k i .
By repeating, successively, the actions with k 3 , , k m , we obtain
k 1 n 1 x 1 δ 1 k m n m x m δ m i = 1 m n i k i x i k i ( 1 x i ) n i k i = i = 1 m k i n i x i δ i n i k i x i k i ( 1 x i ) n i k i .
Let us use the following result which can be found in References [1,21,22]:
k n x δ n k x k ( 1 x ) n k 1 4 n δ 2 ,
where n , k N , n 0 , k n , δ > 0 and x [ 0 , 1 ] . We then get
k 1 n 1 x 1 δ 1 k m n m x m δ m i = 1 m n i k i x i k i ( 1 x i ) n i k i i = 1 m 1 4 n i δ i 2 1 4 m i = 1 m n i δ i 2 .
Remark 1.
Lemma 3 is still valid if the summation is not over all the k i , 1 i m but over some of them. In this case, the upper bound will contain only terms involved in the summation.
The following theorem states for the pointwise and uniform convergence of our Bernstein polynomials with several variables. It should be mentionned that similar work has been done but only for two variables (see (Theorem 1.1 in Reference [26]), [30], (Theorems 3, 5 and 8 in [21]), (Theorem 6.7 in [33]), (Lemma 1 in [28]), (Theorem in [32]) and (Theorems 1 and 3 in [7]).
Theorem 2.
Let f ( x 1 , , x m ) be bounded in the m-dimensional cube 0 x i 1 , i = 1 , , m . Then B n 1 , , n m ( f ; x 1 , , x m ) converges towards f ( x 1 , , x m ) at any point of continuity of this function, as n i for all i. If f C [ 0 , 1 ] m , the limit holds uniformly in [ 0 , 1 ] m .
Proof of Theorem 2.
The function f ( x 1 , , x m ) is assumed bounded in [ 0 , 1 ] m . Hence for some M > 0 , | f ( x 1 , , x m ) | M , ( x 1 , , x m ) [ 0 , 1 ] m , and therefore for any two points ( α 1 , , α m ) , ( β 1 , , β m ) [ 0 , 1 ] m
f ( α 1 , , α m ) f ( β 1 , , β m ) 2 M .
Let ( x 1 , , x m ) [ 0 , 1 ] m be a point of continuity of f. Given ε > 0 , we can find δ 1 , , δ m depending on ( x 1 , , x m ) and ε such that f ( x 1 , , x m ) f ( y 1 , , y m ) < ε 2 , whenever | x i y i | < δ i , i = 1 , , m .
Since B n 1 , , n m ( 1 ; x 1 , , x m ) = 1 , we have
f ( x 1 , , x m ) = f ( x 1 , , x m ) B n 1 , , n m ( 1 ; x 1 , , x m ) = f ( x 1 , , x m ) k 1 = 0 n 1 k m = 0 n m i = 1 m n i k i x i k i ( 1 x i ) n i k i = k 1 = 0 n 1 k m = 0 n m f ( x 1 , , x m ) i = 1 m n i k i x i k i ( 1 x i ) n i k i .
Let us consider the set E = i = 1 m { 0 , , n i } , and for j = 1 , , m , we define the sets
Ω j = k j { 0 , , n j } : k j n j x j < δ j and F = E ( Ω 1 × × Ω m )
Then, F = k = 1 m F k , with F k = i = 1 m Ω i k [ α i k ] F : α i k { 0 , 1 } , i = 1 m α i k = k , where
Ω i k [ α i k ] = Ω i if α i k = 0 Ω i c if α i k = 1 and Ω i c = k i { 0 , , n i } : k i n i x i δ i .
For A k = i = 1 m Ω i k [ α i k ] F k , k = 1 , , m , let us define also I A k = i { 1 , , m } : α i k = 1 (that means card I A k = k 1 ). We have
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) = k 1 = 0 n 1 k m = 0 n m f ( x 1 , , x m ) i = 1 m n i k i x i k i ( 1 x i ) n i k i k 1 = 0 n 1 k m = 0 n m f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i , = k 1 = 0 n 1 k m = 0 n m f ( x 1 , , x m ) f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i , k 1 = 0 n 1 k m = 0 n m f ( x 1 , , x m ) f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i , Ω 1 Ω m f ( x 1 , , x m ) f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i + F f ( x 1 , , x m ) f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i .
Using the fact that f is continuous and bounded, we get
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) ε 2 Ω 1 Ω m i = 1 m n i k i x i k i ( 1 x i ) n i k i + 2 M F i = 1 m n i k i x i k i ( 1 x i ) n i k i , ε 2 k 1 = 0 n 1 k m = 0 n m i = 1 m n i k i x i k i ( 1 x i ) n i k i + 2 M l = 1 m A l F l i = 1 m n i k i x i k i ( 1 x i ) n i k i , since F = l = 1 m F l ε 2 B n 1 , , n m ( 1 ; x 1 , , x m ) + 2 M l = 1 m A l F l i = 1 m n i k i x i k i ( 1 x i ) n i k i , ε 2 + 2 M l = 1 m A l F l i = 1 m n i k i x i k i ( 1 x i ) n i k i ( using Equation ( 31 ) ) .
For l { 1 , , m } and A l F l ,
A l i = 1 m n i k i x i k i ( 1 x i ) n i k i = Ω 1 l [ α 1 l ] Ω m l [ α m l ] i = 1 m n i k i x i k i ( 1 x i ) n i k i = i = 1 m Ω i l [ α i l ] n i k i x i k i ( 1 x i ) n i k i ,
using a similar proof as in Lemma 1 . For i { 1 , , m } ,
-
if α i l = 0 , this means i I A l then,
Ω i l [ α i l ] n i k i x i k i ( 1 x i ) n i k i = Ω i n i k i x i k i ( 1 x i ) n i k i k i = 0 n i n i k i x i k i ( 1 x i ) n i k i = 1
-
if α i l = 1 that means i I A l then
Ω i l [ α i l ] n i k i x i k i ( 1 x i ) n i k i = Ω i c n i k i x i k i ( 1 x i ) n i k i 1 4 n i δ i 2
(using Equation (45)).
Hence, for i , l { 1 , , m } ,
Ω i l [ α i l ] n i k i x i k i ( 1 x i ) n i k i 1 if i I A l , 1 4 n i δ i 2 if i I A l ,
and this implies
A l i = 1 m n i k i x i k i ( 1 x i ) n i k i i I A l 1 4 n i δ i 2 .
Therefore,
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) ε 2 + 2 M l = 1 m A l F l i I A l 1 4 n i δ i 2 = ε 2 + l = 1 m A l F l 2 M i I A l 4 n i δ i 2 , ε 2 + l = 1 m A l F l 2 M n l i I A l 4 δ i 2 ( because card I A l = l ) , ε 2 + 1 n l = 1 m A l F l 2 M i I l 4 δ i 2 ( with n = min n i : i = 1 , , m ) .
For n 1 , , n m sufficiently large (that means n sufficiently large), we have 1 n l = 1 m A l F l 2 M i I l δ i 2 < ε 2 . Therefore,
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) ε 2 + ε 2 = ε .
Here, it should be noticed that the minimum integer n ˜ , such that
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) ε n 1 , , n m n n ˜ ,
depends on δ 1 , , δ m which themselves depend on x 1 , , x m and ε . That is why we only have a pointwise convergence.
Suppose now that f C [ 0 , 1 ] m , then f is uniformly continuous on [ 0 , 1 ] m that is, given ε > 0 , we can find δ 1 , , δ m depending just on ε , but not on any element in [ 0 , 1 ] , such that f ( x 1 , , x m ) f ( y 1 , , y m ) < ε 2 , for all ( x 1 , , x m ) , ( y 1 , , y m ) [ 0 , 1 ] m satisfying | x i y i | < δ i , i = 1 , , m . By following the same process, we obtain that
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) ε 2 + 1 n l = 1 m 2 M i I l δ i 2
for all x 1 , , x m [ 0 , 1 ] and n = min n i : i = 1 , , m . Hence, there exist n ˜ depending only on ε and not on x 1 , , x m such that
f ( x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) < ε , x 1 , , x m [ 0 , 1 ]
and n 1 , , n m n n ˜ . Therefore, the convergence is uniform in [ 0 , 1 ] m . □
We express now, in Lemma 3 , the p 1 , , p m -th derivative of our multivariate Bernstein polynomials of order n 1 + p 1 , , n m + p m -th in terms of p 1 , , p m -th differences of the corresponding functions. A similar result exists for one variable (see References [1,15]).
Lemma 3.
For p 1 , , p m 0 integers with 0 p i n i , i = 1 , 2 , , m , the p 1 , , p m -th derivative of B n 1 + p 1 , , n m + p m ( f ; x 1 , , x m ) may be expressed in terms of p 1 , , p m -th differences of f as
B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) = i = 1 m ( n i + p i ) ! n i ! t 1 = 0 n 1 t m = 0 n m 1 p 1 m p m f t 1 n 1 + p 1 , , t m n m + p m i = 1 m n i t i x i t i ( 1 x i ) n i t i .
Proof of Lemma 3. 
We will sometimes use the square brackets to draw attention on computations of certain terms. We write
B n 1 + p 1 , , n m + p m ( f ; x 1 , , x m ) = k 1 = 0 n 1 + p 1 k m = 0 n m + p m f k 1 n 1 + p 1 , , k m n m + p m i = 1 m n i k i x i k i ( 1 x i ) n i k i , = k 2 = 0 n 2 + p 2 k m = 0 n m + p m k 1 = 0 n 1 + p 1 f k 1 n 1 + p 1 , , k m n m + p m n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 × i = 2 m n i k i x i k i ( 1 x i ) n i k i .
We take the first variable of f and fix the others. Hence, f k 1 n 1 + p 1 , , k m n m + p m is seen as a function g 1 k 1 n 1 + p 1 :
B n 1 + p 1 , , n m + p m ( f ; x 1 , , x m ) = k 2 = 0 n 2 + p 2 k m = 0 n m + p m k 1 = 0 n 1 + p 1 g 1 k 1 n 1 + p 1 n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 2 = 0 n 2 + p 2 k m = 0 n m + p m B n 1 + p 1 ( g 1 , x 1 ) i = 2 m n i k i x i k i ( 1 x i ) n i k i .
Hence
B n 1 + p 1 , , n m + p m p 1 ( f ; x 1 , , x m ) = k 2 = 0 n 2 + p 2 k m = 0 n m + p m B n 1 + p 1 p 1 ( g 1 , x 1 ) i = 2 m n i k i x i k i ( 1 x i ) n i k i .
Let us use the following result from P. J. Davis [1]:
B n + p ( p ) ( f ; x ) = ( n + p ) ! n ! t = 0 n p f n + p t n t x t ( 1 x ) n t ,
where n , p N and f n + p t = f t n + p . We then get
B n 1 + p 1 , , n m + p m p 1 ( f ; x 1 , , x m ) = k 2 = 0 n 2 + p 2 k m = 0 n m + p m ( n 1 + p 1 ) ! n ! t 1 = 0 n 1 1 p 1 g 1 t 1 n 1 + p 1 n 1 t 1 x 1 t 1 ( 1 x 1 ) n 1 t 1 i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 2 = 0 n 2 + p 2 k m = 0 n m + p m ( n 1 + p 1 ) ! n ! t 1 = 0 n 1 1 p 1 f t 1 n 1 + p 1 , k 2 n 2 + p 2 , , k m n m + p m n 1 t 1 x 1 t 1 ( 1 x 1 ) n 1 t 1 × i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 3 = 0 n 3 + p 3 k m = 0 n m + p m ( n 1 + p 1 ) ! n ! t 1 = 0 n 1 k 2 = 0 n 2 + p 2 1 p 1 f t 1 n 1 + p 1 , k 2 n 2 + p 2 , , k m n m + p m
× n 2 k 2 x 2 k 2 ( 1 x 2 ) n 2 k 2 n 1 t 1 x 1 t 1 ( 1 x 1 ) n 1 t 1 i = 3 m n i k i x i k i ( 1 x i ) n i k i .
We repeat the action with the second variable of the function 1 p 1 f t 1 n 1 + p 1 , k 2 n 2 + p 2 , , k m n m + p m , we get
B n 1 + p 1 , , n m + p m p 1 , p 2 ( f ; x 1 , , x m ) = k 3 = 0 n 3 + p 3 k m = 0 n m + p m i = 1 2 ( n i + p i ) ! n i ! t 1 = 0 n 1 t 2 = 0 n 2 1 p 1 2 p 2 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , k 3 n 3 + p 3 , , k m n m + p m i = 1 2 n i t i x i t i ( 1 x i ) n i t i i = 3 m n i k i x i k i ( 1 x i ) n i k i , = k 3 = 0 n 3 + p 3 k m = 0 n m + p m i = 1 2 ( n i + p i ) ! n i ! t 1 = 0 n 1 t 2 = 0 n 2 k 3 = 0 n 3 + p 3 1 p 1 2 p 2 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , k 3 n 3 + p 3 , , k m n m + p m n 3 k 3 x 3 k 3 ( 1 x 3 ) n 3 k 3 i = 1 2 n i t i x i t i ( 1 x i ) n i t i i = 4 m n i k i x i k i ( 1 x i ) n i k i .
After, we continue the same action with the third variable of 1 t 1 2 t 2 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , k 3 n 3 + p 3 , , k m n m + p m and get
B n 1 + p 1 , , n m + p m p 1 , p 2 , p 3 ( f ; x 1 , , x m ) = k 4 = 0 n 4 + p 4 k m = 0 n m + p m i = 1 3 ( n i + p i ) ! n i ! t 2 = 0 n 2 1 p 1 2 p 2 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , t 3 n 3 + p 3 , k 4 n 4 + p 4 , , k m n m + p m i = 1 3 n i t i x i t i ( 1 x i ) n i k i i = 4 m n i k i x i k i ( 1 x i ) n i k i ,
which can be written as
B n 1 + p 1 , , n m + p m p 1 , p 2 , p 3 ( f ; x 1 , , x m ) = k 5 = 0 n 5 + p 5 k m = 0 n m + p m i = 1 3 ( n i + p i ) ! n i ! t 1 = 0 n 1 t 2 = 0 n 2 t 3 = 0 n 3 k 4 = 0 n 4 + p 4 1 p 1 2 p 2 3 p 3 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , t 3 n 3 + p 3 , k 4 n 4 + p 4 , , k m n m + p m n 4 t 4 x 4 t 4 ( 1 x 4 ) n 4 k 4 i = 1 3 n i t i x i t i ( 1 x i ) n i k i × i = 5 m n i k i x i k i ( 1 x i ) n i k i .
By repeating successively the actions with the 4 th , 5 th , , m th variables but with the functions 1 t 1 2 t 2 3 t 3 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , t 3 n 3 + p 3 , k 4 n 4 + p 4 , , k m n m + p m , 1 t 1 2 t 2 3 t 3 4 t 4 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , t 3 n 3 + p 3 , t 4 n 4 + p 4 , k 5 n 5 + p 5 , , k m n m + p m , , 1 t 1 2 t 2 m 1 t m 1 f t 1 n 1 + p 1 , t 2 n 2 + p 2 , , t m 1 n m 1 + p m 1 , k m n m + p m respectively, we obtain the result.  □
Let us see how to deduce the monotonicity, the fixed sign of the p 1 , , p m -th derivative, the convexity, and the upper and lower bounds of our multivariate Bernstein polynomials from those of the corresponding function. Similar results exist for one variable (see References [1,21,31]).
Theorem 3.
Let 0 p i n i , i = 1 , , m , be fixed integers and the restriction of f to each of its variables i belong to C p i [ 0 , 1 ] , i = 1 , , m . If
A f ( p 1 ) ( p m ) ( x 1 , , x m ) B , ( x 1 , , x m ) [ 0 , 1 ] m ,
t h e n A i = 1 m n i p i n i ( n i 1 ) ( n i p i + 1 ) B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) B
with ( x 1 , , x m ) [ 0 , 1 ] m . For all p i = 0 , i = 1 , , m , the multiplier of B n 1 , , n m p 1 , , p m is to be interpreted as 1 .
-
If f ( p 1 ) ( p m ) ( x 1 , , x m ) 0 , ( x 1 , , x m ) [ 0 , 1 ] m , then
B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) 0 , ( x 1 , , x m ) [ 0 , 1 ] m
-
If f ( x 1 , , x m ) is convex on [ 0 , 1 ] m , then B n 1 , , n m ( f ; x 1 , , x m ) is convex on [ 0 , 1 ] m .
Proof of Theorem 3.
We begin with Equation (62) and replace n i by n i p i , i = 1 , , m . We have
B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) = i = 1 m n i ! ( n i p i ) ! t 1 = 0 n 1 p 1 t m = 0 n m p m 1 p 1 m p m f t 1 n 1 , , t m n m × i = 1 m n i t i x i t i ( 1 x i ) n i p i t i .
Let us use the following result from P. J. Davis [1]: for g C p [ a , b ] and a = e 0 < e 1 < < e n = b with e i = e 0 + i h R , h = e i + 1 e i , a subdivision of [ a , b ] ,
ε 0 ] e 0 , e p [ , p g ( x 0 ) h p = g ( p ) ( ε 0 ) .
Then we get for any variable of f with step size h i = 1 n i , i = 1 , , m ,
1 p 1 m p m f t 1 n 1 , , t m n m = f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) n 1 p 1 n m p m ,
where t i n i < ε t i < t i + p i n i , i = 1 , , m . Thus,
B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) = i = 1 m n i ! ( n i p i ) ! t 1 = 0 n 1 p 1 t m = 0 n m p m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) n 1 p 1 n m p m i = 1 m n i p i t i x i t i ( 1 x i ) n i p i t i = i = 1 m n i ( n i 1 ) ( n i p i + 1 ) n i p i t 1 = 0 n 1 p 1 t m = 0 n m p m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) × i = 1 m n i p i t i x i t i ( 1 x i ) n i p i t i .
A f ( p 1 ) ( p m ) ( x 1 , , x m ) B , ( x 1 , , x m ) [ 0 , 1 ] m
A f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) B
A f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) t 1 = 0 n 1 p 1 t m = 0 n m p m i = 1 m n i p i t i x i t i ( 1 x i ) n i p i t i B
since t 1 = 0 n 1 p 1 t m = 0 n m p m i = 1 m n i p i t i x i t i ( 1 x i ) n i p i t i = 1
A t 1 = 0 n 1 p 1 t m = 0 n m p m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) i = 1 m n i p i t i x i t i ( 1 x i ) n i p i t i B
A i = 1 m n i p i n i ( n i 1 ) ( n i p i + 1 ) B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) B for all
( x 1 , , x m ) [ 0 , 1 ] m ( using ( 77 ) ) .
When all p i = 0 , the multiplier of B n 1 , , n m p 1 , , p m is
i = 1 m n i p i n i ( n i 1 ) ( n i p i + 1 ) = i = 1 m ( n i p i ) ! n p i n i ! = 1 ( p i = 0 i = 1 , , m ) .
*
If f ( p 1 ) ( p m ) ( x 1 , , x m ) 0 , ( x 1 , , x m ) [ 0 , 1 ] m , we set A = 0 and obtain
0 i = 1 m n i p i n i ( n i 1 ) ( n i p i + 1 ) B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m )
for all ( x 1 , , x m ) [ 0 , 1 ] m . This means B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) 0 , ( x 1 , , x m ) [ 0 , 1 ] m .
*
If f is convex, 1 2 m 2 f t 1 n 1 , , t m n m 0 and hence we have f ( 2 ) ( 2 ) ( ε t 1 , , ε t m ) n 1 2 n m 2 0 since f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) n 1 p 1 n m p m = 1 p 1 m p m f t 1 n 1 , , t m n m , this means f ( 2 ) ( 2 ) ( ε t 1 , , ε t m ) 0 . Thus, from (77) with p i = 2 i = 1 , , m ,
B n 1 , , n m 2 , , 2 ( f ; x 1 , , x m ) 0 , ( x 1 , , x m ) [ 0 , 1 ] m .
This implies that B n 1 , , n m is convex in every closed interval of ( 0 , 1 ) m . Since B n 1 , , n m is continuous, it is convex in [ 0 , 1 ] m .  □
Corollary 1.
If f ( x 1 , , x m ) is decreasing (resp. nondecreasing) on [ 0 , 1 ] with respect to the variable x i , i { 1 , , m } , then B n 1 , , n m ( f ; x 1 , , x m ) is decreasing (resp. nondecreasing) with respect to the variable x i , i { 1 , , m } on [ 1 , 1 ] .
Proof of Corollary 1.
We just use Theorem 3 , and the rest follows easily. □
The following theorem states for the pointwise and uniform convergence of the p 1 , , p m th derivative of our multivariate Bernstein polynomials, compare ([35] Theorem of Kingsley) and ([28] Theorem 1) (just for two variables), and ([25] Theorem 4) and ([7] Theorem 2).
Theorem 4.
Let 0 p i n i , i = 0 , , m , be fixed integers, f uniformly continuous on [ 0 , 1 ] m and p i f x i p i also uniformly continuous on [ 0 , 1 ] m , i = 1 , , m . Then B n 1 , , n m p 1 , , p m ( f ; x 1 , , x m ) converges towards f ( p 1 ) ( p m ) ( x 1 , , x m ) uniformly on [ 0 , 1 ] m .
Proof of Theorem 4. 
We have already shown that the the above result holds for p i = 0 , i = 0 , , m . We have to look at the case where p i are not all zero.
We begin with the expression for B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) given in Equation (62) and write using Equation ( 75 ) with h i = 1 n i + p i , i = 1 , , m
1 p 1 m p m f t 1 n 1 + p 1 , , t m n m + p m = f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) i = 1 m ( n i + p i ) p i ,
where t i n i + p i < ε t i < t i + p i n i + p i , i = 1 , , m . We get
B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) = i = 1 m ( n i + p i ) ! n i ! t 1 = 0 n 1 t m = 0 n m 1 p 1 m p m f t 1 n 1 + p 1 , , t m n m + p m i = 1 m n i t i x i t i ( 1 x i ) n i t i = i = 1 m ( n i + p i ) ! n i ! ( n i + p i ) p i t 1 = 0 n 1 t m = 0 n m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) i = 1 m n i t i x i t i ( 1 x i ) n i t i ,
which gives us
i = 1 m n i ! ( n i + p i ) p i ( n i + p i ) ! B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) i = 1 m n i t i x i t i ( 1 x i ) n i t i .
We then approximate f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) writing
f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) = f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) f ( p 1 ) ( p m ) t 1 n 1 , , t m n m + f ( p 1 ) ( p m ) t 1 n 1 , , t m n m .
We thus obtain
i = 1 m n i ! ( n i + p i ) p i ( n i + p i ) ! B n 1 + p 1 , , n m + p m p 1 , , p m ( f ; x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m f ( p 1 ) ( p m ) t 1 n 1 , , t m n m + f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) f ( p 1 ) ( p m ) t 1 n 1 , , t m n m i = 1 m n i t i x i t i ( 1 x i ) n i t i = T 1 ( x 1 , , x m ) + T 2 ( x 1 , , x m ) where
T 1 ( x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m f ( p 1 ) ( p m ) t 1 n 1 , , t m n m × i = 1 m n i t i x i t i ( 1 x i ) n i t i , T 2 ( x 1 , , x m ) = t 1 = 0 n 1 t m = 0 n m f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) f ( p 1 ) ( p m ) t 1 n 1 , , t m n m i = 1 m n i t i x i t i ( 1 x i ) n i t i .
Since t i n i + p i < t i n i < t i + p i n i + p i , i = 1 , , m , t i = 0 , , n i , it follows from the bounds on ε t i , i = 1 , , m , that
p i n i + p i = t i n i + p i t i + p i n i + p i < ε t i t i n i < t i + p i n i + p i t i n i + p i = p i n i + p i
for i = 1 , , m . This means ε t i t i n i < p i n i + p i , i = 1 , , m . From the uniform continuity of f ( p 1 ) ( p m ) ( x 1 , , x m ) , given ε > 0 , we can find n i 0 such that for all n i n i 0 and all t i ,
f ( p 1 ) ( p m ) ( ε t 1 , , ε t m ) f ( p 1 ) ( p m ) t 1 n 1 , , t m n m < ε .
Thus, T 2 ( x 1 , , x m ) converges uniformly to zero on [ 0 , 1 ] m . We have
i = 1 m n i ! ( n i + p i ) p i ( n i + p i ) ! 1 as n i , i = 1 , , m ,
and we see from Theorem 2 ) with f ( p 1 ) ( p m ) in place of f that T 1 ( x 1 , , x m ) converges uniformly to f ( p 1 ) ( p m ) ( x 1 , , x m ) . This completes the proof. □
The following theorem, in the case where f is a convex function, is the generalization of the inequality related to the decreasing of the sequence of Bernstein’s polynomials with one variable (see References [1,31]).
Theorem 5.
Let f ( x 1 , , x m ) be convex in [ 0 , 1 ] m . Then, for n i = 2 , 3 , , i = 1 , , m , we have
B n 1 1 , , n m 1 ( f ; x 1 , , x m ) B n 1 , , n m ( f ; x 1 , , x m ) , ( x 1 , , x m ) [ 0 , 1 ] m .
Proof of Theorem 5.
The case of one variable is already solved (see P. J. Davis [1]). Because of long calculations and space limitations, we will just show how this theorem works for three variables (we did not choose the case of two variables because with three variables, our strategy is well illustrated and we can easily use it for any large number of variables). Set t i = x i 1 x i , i = 1 , 2 , 3 . Then,
x i k i ( 1 x i ) k i = t i k i 1 + t i = ( 1 x i ) 1 , t h a t i s , x i k i ( 1 x i ) k i = t i k i x i k i ( 1 x i ) k i 1 = t i k i ( 1 + t i ) ,
and
( 1 x 1 ) n 1 ( 1 x 2 ) n 2 ( 1 x 3 ) n 3 B n 1 1 , n 2 1 , n 3 1 ( f ; x 1 , x 2 , x 3 ) B n 1 , n 2 , n 3 ( f ; x 1 , x 2 , x 3 ) = k 1 = 0 n 1 1 k 2 = 0 n 2 1 k 3 = 0 n 3 1 f k 1 n 1 1 , k 2 n 2 1 , k 3 n 3 1 i = 1 3 n i 1 k i x i k i ( 1 x i ) k i 1 k 1 = 0 n 1 k 2 = 0 n 2 k 3 = 0 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 i = 1 3 n i k i x i k i ( 1 x i ) k i = k 1 = 0 n 1 1 k 2 = 0 n 2 1 k 3 = 0 n 3 1 f k 1 n 1 1 , k 2 n 2 1 , k 3 n 3 1 i = 1 3 n i 1 k i t i k i ( 1 + t i ) k 1 = 0 n 1 k 2 = 0 n 2 k 3 = 0 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 i = 1 3 n i k i t i k i = i 1 , i 2 , i 3 { 0 , 1 } k 1 = i 1 n 1 1 + i 1 k 2 = i 2 n 2 1 + i 2 k 3 = i 3 n 3 1 + i 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 t 1 k 1 t 2 k 2 t 3 k 3 k 1 = 0 n 1 k 2 = 0 n 2 k 3 = 0 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 n 1 k 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 = i 1 , i 2 , i 3 { 0 , 1 } k 1 = i 1 n 1 1 + i 1 k 2 = i 2 n 2 1 + i 2 k 3 = i 3 n 3 1 + i 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 t 1 k 1 t 2 k 2 t 3 k 3 k 1 = 0 n 1 k 2 = 0 n 2 k 3 = 0 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 n 1 k 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 = F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) .
F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) is a polynomial of the variables t 1 , t 2 and t 3 . It can be written as
F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) = k 1 = 0 n 1 k 2 = 0 n 2 k 3 = 0 n 3 F n 1 , n 2 , n 3 k 1 , k 2 , k 3 t 1 k 1 t 2 k 2 t 3 k 3 .
To have our result, we just need to prove that F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) is positive.
For i e { 0 , 1 } with e { 1 , 2 , 3 } , we define
i e ¯ = 0 if i e = 1 ; 1 if i e = 0 .
Now, we will analyse all the possibilities of the summand F n 1 , n 2 , n 3 k 1 , k 2 , k 3 t 1 k 1 t 2 k 2 t 3 k 3 of our polynomial F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) . For simplicity and depending of the upper and lower bounds of ( 99 ) , we will write F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) as
F n 1 , n 2 , n 3 ( t 1 , t 2 , t 3 ) = F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) ,
where
*
F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) is the sum with index 1 k 1 n 1 1 , 1 k 2 n 2 1 and 1 k 3 n 3 1 ;
*
F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) is the sum where all the index k i , i { 1 , 2 , 3 } start at 1, but with just one index reaching the upper bound n i , i { 1 , 2 , 3 } ;
*
F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) is the sum where all the index k i are less or equal to n i 1 , i { 1 , 2 , 3 } , but just one index k i , i { 1 , 2 , 3 } reaching the lower bound 0;
*
F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) is the rest of the sum in Equation ( 99 ) after removing the terms F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) , F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) and F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) .
-
Case 1: First, we take the sum with index 1 k 1 n 1 1 , 1 k 2 n 2 1 and 1 k 3 n 3 1 in Equation ( 99 ) . We obtain
F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 , i 2 , i 3 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 n 1 k 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 .
Next, we add and subtract the following expression in F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) to obtain:
k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 , i 2 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 k 3 + i 1 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 n 2 , k 3 n 3 n 1 1 k 1 i 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 .
Next, we group now the terms in F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) given by the previous expression in the following three blocks: F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) = B 1 + B 2 + B 3 with
B 1 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 , i 2 , i 3 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 + f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 i 3 ¯ n 3 1 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 ¯ f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 , B 2 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 , i 2 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 k 3 + f k 1 i 1 n 1 1 , k 2 i 2 ¯ n 2 1 , k 3 n 3 n 1 1 k 1 i 1 n 2 1 k 2 i 2 ¯ n 3 k 3
f k 1 i 1 n 1 1 , k 2 n 2 , k 3 n 3 n 1 1 k 1 i 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 , B 3 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 n 2 , k 3 n 3 n 1 1 k 1 i 1 n 2 k 2 n 3 k 3 + f k 1 i 1 ¯ n 1 1 , k 2 n 2 , k 3 n 3 n 1 1 k 1 i 1 ¯ n 2 k 2 n 3 k 3
f k 1 n 1 , k 2 n 2 , k 3 n 3 n 1 k 1 n 2 k 2 n 3 k 3 t 1 k 1 t 2 k 2 t 3 k 3 .
We are going to show that B i 0 for i { 1 , 2 , 3 } . It is obvious that
n 1 k = n k n k n ( n k + 1 ) and n 1 k 1 = n k k n ( n k ) with k , n positive integers .
Now, using Equation ( 107 ) in each B i we get
B 1 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 , i 2 { 0 , 1 } n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 k 3 n 3 k 3 n 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 1 + k 3 n 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 1 n 3 1 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 + i 1 , i 2 { 0 , 1 } n 1 1 k 1 i 1 n 2 1 k 2 i 2 n 3 k 3 k 3 n 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 1 n 3 1 + n 3 k 3 n 3 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 1 f k 1 i 1 n 1 1 , k 2 i 2 n 2 1 , k 3 n 3 t 1 k 1 t 2 k 2 t 3 k 3 ,
B 2 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 1 { 0 , 1 } n 1 1 k 1 i 1 n 2 k 2 n 3 k 3 n 2 k 2 n 2 f k 1 i 1 n 1 1 , k 2 n 2 1 , k 3 n 3 + k 2 n 2 f k 1 i 1 n 1 1 , k 2 1 n 2 1 , k 3 n 3 f k 1 i 1 n 1 1 , k 2 n 2 , k 3 n 3 + i 1 { 0 , 1 } n 1 1 k 1 i 1 n 2 k 2 n 3 k 3 k 2 n 2 f k 1 i 1 n 1 1 , k 2 1 n 2 1 , k 3 n 3 + n 2 k 2 n 2 f k 1 i 1 n 1 1 , k 2 n 2 1 , k 3 n 3 f k 1 i 1 n 1 1 , k 2 n 2 , k 3 n 3 t 1 k 1 t 2 k 2 t 3 k 3 ,
B 3 = k 1 = 1 n 1 1 k 2 = 1 n 2 1 k 3 = 1 n 3 1 n 1 k 1 n 2 k 2 n 3 k 3 n 1 k 1 n 1 f k 1 n 1 1 , k 2 n 2 , k 3 n 3 + k 1 n 1 f k 1 1 n 1 1 , k 2 n 2 , k 3 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 + n 1 k 1 n 2 k 2 n 3 k 3 k 1 n 1 f k 1 1 n 1 1 , k 2 n 2 , k 3 n 3 + n 1 k 1 n 1 f k 1 n 1 1 , k 2 n 2 , k 3 n 3 f k 1 n 1 , k 2 n 2 , k 3 n 3 t 1 k 1 t 2 k 2 t 3 k 3 ,
which implies that B i 0 , i { 1 , 2 , 3 } by the convexity of f and the fact that t 1 , t 2 and t 3 are also positive.
-
Case 2: We take the sum in Equation ( 99 ) , where all the index start at 1, but just one index k i , i { 1 , 2 , 3 } reaches the upper bound n i , i { 1 , 2 , 3 } . This corresponds to F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) which can be written as
F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) = F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 2 , 2 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 2 , 3 ( t 1 , t 2 , t 3 )
where F n 1 , n 2 , n 3 2 , i ( t 1 , t 2 , t 3 ) is the sum where the index k i reaches the upper bound n i for i { 1 , 2 , 3 } . We have to show that F n 1 , n 2 , n 3 2 , i ( t 1 , t 2 , t 3 ) , i { 1 , 2 , 3 } , are positive. Let us start with F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) :
F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) = k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 , i 3 { 0 , 1 } f 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 f 1 , k 2 n 2 , k 3 n 3 n 2 k 2 n 3 k 3 t 1 n 1 t 2 k 2 t 3 k 3 .
We add and subtract the following expression in F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 )
k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 { 0 , 1 } f 1 , k 2 i 2 n 2 1 , k 3 n 3 n 2 1 k 2 i 2 n 3 k 3 t 1 n 1 t 2 k 2 t 3 k 3 .
We now group the terms in F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) in the following two blocks: F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) = B 1 + B 2 with
B 1 = k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 , i 3 { 0 , 1 } f 1 , k 2 i 2 n 2 1 , k 3 i 3 n 3 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 + f 1 , k 2 i 2 n 2 1 , k 3 i 3 ¯ n 3 1 n 2 1 k 2 i 2 n 3 1 k 3 i 3 ¯ f 1 , k 2 i 2 n 2 1 , k 3 n 3 n 2 1 k 2 i 2 n 3 k 3 t 1 n 1 t 2 k 2 t 3 k 3 ,
B 2 = k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 { 0 , 1 } f 1 , k 2 i 2 n 2 1 , k 3 n 3 n 2 1 k 2 i 2 n 3 k 3 + f 1 , k 2 i 2 ¯ n 2 1 , k 3 n 3 n 2 1 k 2 i 2 ¯ n 3 k 3 f 1 , k 2 n 2 , k 3 n 3 n 2 k 2 n 3 k 3 t 1 n 1 t 2 k 2 t 3 k 3 .
We are going to show that B i 0 for i { 1 , 2 } . Using ( 107 ) in each B i , we get
B 1 = k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 { 0 , 1 } n 2 1 k 2 i 2 n 3 k 3 n 3 k 3 n 3 f 1 , k 2 i 2 n 2 1 , k 3 n 3 1 + k 3 n 3 f 1 , k 2 i 2 n 2 1 , k 3 1 n 3 1 f 1 , k 2 i 2 n 2 1 , k 3 n 3 + i 2 { 0 , 1 } n 2 1 k 2 i 2 n 3 k 3 k 3 n 3 f 1 , k 2 i 2 n 2 1 , k 3 1 n 3 1 + n 3 k 3 n 3 f 1 , k 2 i 2 n 2 1 , k 3 n 3 1 f 1 , k 2 i 2 n 2 1 , k 3 n 3 t 1 n 1 t 2 k 2 t 3 k 3 ,
B 2 = k 2 = 1 n 2 1 k 3 = 1 n 3 1 n 2 k 2 n 3 k 3 n 2 k 2 n 2 f 1 , k 2 n 2 1 , k 3 n 3 + k 2 n 2 f 1 , k 2 1 n 2 1 , k 3 n 3 f 1 , k 2 n 2 , k 3 n 3 t 1 n 1 t 2 k 2 t 3 k 3 ,
which implies that B i 0 , i { 1 , 2 } by the convexity of f and the fact that t 1 , t 2 and t 3 are also positive.
Using a similar reasoning as for F n 1 , n 2 , n 3 2 , 1 ( t 1 , t 2 , t 3 ) , one shows that F n 1 , n 2 , n 3 2 , 2 ( t 1 , t 2 , t 3 ) and F n 1 , n 2 , n 3 2 , 3 ( t 1 , t 2 , t 3 ) are also positive. The main difference appears in the expression that we add and subtract:
(i)
for F n 1 , n 2 , n 3 2 , 2 ( t 1 , t 2 , t 3 ) , the added and subtracted expression is
k 1 = 1 n 1 1 k 3 = 1 n 3 1 i 1 { 0 , 1 } f k 1 i 1 n 1 1 , 1 , k 3 n 3 n 1 1 k 1 i 1 n 3 k 3 t 1 k 1 t 2 n 2 t 3 k 3 ;
(ii)
for F n 1 , n 2 , n 3 2 , 3 ( t 1 , t 2 , t 3 ) , the added and subtracted expression is
k 1 = 1 n 1 1 k 2 = 1 n 2 1 i 1 { 0 , 1 } f k 1 i 1 n 1 1 , k 2 n 2 , 1 n 1 1 k 1 i 1 n 2 k 2 t 1 k 1 t 2 k 2 t 3 n 3 .
-
Case 3: We take the sum in Equation ( 99 ) , where all the index k i are less or equal to n i 1 , i { 1 , 2 , 3 } , but just one index k i , i { 1 , 2 , 3 } reaching the lower bound 0. This corresponds to F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) which can be written as
F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) = F n 1 , n 2 , n 3 3 , 1 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 3 , 2 ( t 1 , t 2 , t 3 ) + F n 1 , n 2 , n 3 3 , 3 ( t 1 , t 2 , t 3 ) ,
where F n 1 , n 2 , n 3 3 , i ( t 1 , t 2 , t 3 ) is the sum where the index k i reaches the lower bound 0 for i { 1 , 2 , 3 } .
We have
F n 1 , n 2 , n 3 3 , 1 ( t 1 , t 2 , t 3 ) = k 2 = 1 n 2 1 k 3 = 1 n 3 1 F n 1 , n 2 , n 3 0 , k 2 , k 3 t 2 k 2 t 3 k 3 ; F n 1 , n 2 , n 3 3 , 1 ( t 1 , t 2 , t 3 ) = k 1 = 1 n 1 1 k 3 = 1 n 3 1 F n 1 , n 2 , n 3 k 1 , 0 , k 3 t 1 k 1 t 3 k 3 ; F n 1 , n 2 , n 3 3 , 1 ( t 1 , t 2 , t 3 ) = k 1 = 1 n 1 1 k 2 = 1 n 2 1 F n 1 , n 2 , n 3 k 1 , k 2 , 0 t 1 k 1 t 2 k 2 .
Using a similar reasoning as in case 2, one shows that F n 1 , n 2 , n 3 3 , i ( t 1 , t 2 , t 3 ) , i { 1 , 2 , 3 } , are positive. The main difference appears in the expression that we add and subtract:
(i)
for F n 1 , n 2 , n 3 3 , 1 ( t 1 , t 2 , t 3 ) , the added and subtracted expression is
k 2 = 1 n 2 1 k 3 = 1 n 3 1 i 2 { 0 , 1 } f 0 , k 2 i 2 n 2 1 , k 3 n 3 n 2 1 k 2 i 2 n 3 k 3 t 2 k 2 t 3 k 3 ;
(ii)
for F n 1 , n 2 , n 3 3 , 2 ( t 1 , t 2 , t 3 ) , the added and subtracted expression is
k 1 = 1 n 1 1 k 3 = 1 n 3 1 i 1 { 0 , 1 } f k 1 i 1 n 2 1 , 0 , k 3 n 3 n 1 1 k 1 i 1 n 3 k 3 t 1 k 1 t 3 k 3 ;
(iii)
for F n 1 , n 2 , n 3 3 , 3 ( t 1 , t 2 , t 3 ) , the added and subtracted expression is
k 1 = 1 n 1 1 k 2 = 1 n 2 1 i 3 { 0 , 1 } f k 1 i 1 n 2 1 , k 2 n 2 , 0 n 1 1 k 1 i 1 n 2 k 2 t 1 k 1 t 2 k 2 .
-
Case 4: We take in Equation ( 99 ) the rest of the sum after removing F n 1 , n 2 , n 3 1 ( t 1 , t 2 , t 3 ) , F n 1 , n 2 , n 3 2 ( t 1 , t 2 , t 3 ) and F n 1 , n 2 , n 3 3 ( t 1 , t 2 , t 3 ) . This corresponds to F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) and has all its terms related with the following monomials: t e k e t j n j t p n p , t e n e t j n j t p n p , t j 0 t e 0 t p 0 , t e k e t j n j , t j k j in ( 99 ) where e , j , p { 1 , 2 , 3 } , e j , e p , j p , and k e { 1 , , n e 1 } . In this case we do not need to add and subtract some expressions in F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) which can be written as
F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) = k 1 = 1 n 1 1 n 1 k 1 n 1 k 1 n 1 f k 1 n 1 1 , 0 , 0 + k 1 n 1 f k 1 1 n 1 1 , 0 , 0 f k 1 n 1 , 0 , 0 t 1 k 1 + k 2 = 1 n 2 1 n 2 k 2 n 2 k 2 n 2 f 0 , k 2 n 2 1 , 0 + k 2 n 2 f 0 , k 2 1 n 2 1 , 0 f 0 , k 2 n 2 , 0 t 2 k 2 + k 1 = 1 n 1 1 n 3 k 3 n 3 k 3 n 3 f 0 , 0 , k 3 n 3 1 + k 3 n 3 f 0 , 0 , k 3 1 n 3 1 f 0 , 0 , k 3 n 3 t 3 k 3 + k 1 = 1 n 1 1 n 1 k 1 n 1 k 1 n 1 f k 1 n 1 1 , 1 , 1 + k 1 n 1 f k 1 1 n 1 1 , 1 , 1 f k 1 n 1 , 1 , 1 t 1 k 1 t 2 n 2 t 3 n 3 + k 2 = 1 n 2 1 n 2 k 2 n 2 k 2 n 2 f 1 , k 2 n 2 1 , 1 + k 2 n 2 f 1 , k 2 1 n 2 1 , 1 f 1 , k 2 n 2 , 1 t 1 n 1 t 2 k 2 t 3 n 3 + k 1 = 1 n 1 1 n 3 k 3 n 3 k 3 n 3 f 1 , 1 , k 3 n 3 1 + k 3 n 3 f 1 , 1 , k 3 1 n 3 1 f 1 , 1 , k 3 n 3 t 1 n 1 t 2 n 2 t 3 k 3 + k 1 = 1 n 1 1 n 1 k 1 n 1 k 1 n 1 f k 1 n 1 1 , 1 , 0 + k 1 n 1 f k 1 1 n 1 1 , 1 , 0 f k 1 n 1 , 1 , 0 t 1 k 1 t 2 n 2 + k 1 = 1 n 1 1 n 1 k 1 n 1 k 1 n 1 f k 1 n 1 1 , 0 , 1 + k 1 n 1 f k 1 1 n 1 1 , 0 , 1 f k 1 n 1 , 0 , 1 t 1 k 1 t 3 n 3 + k 2 = 1 n 2 1 n 2 k 2 n 2 k 2 n 2 f 1 , k 2 n 2 1 , 0 + k 2 n 2 f 1 , k 2 1 n 2 1 , 0 f 1 , k 2 n 2 , 0 t 1 n 1 t 2 k 2 + k 2 = 1 n 2 1 n 2 k 2 n 2 k 2 n 2 f 0 , k 2 n 2 1 , 1 + k 2 n 2 f 0 , k 2 1 n 2 1 , 1 f 0 , k 2 n 2 , 1 t 2 k 2 t 3 n 3 + k 1 = 1 n 1 1 n 3 k 3 n 3 k 3 n 3 f 1 , 0 , k 3 n 3 1 + k 3 n 3 f 1 , 0 , k 3 1 n 3 1 f 1 , 0 , k 3 n 3 t 1 n 1 t 3 k 3 + k 1 = 1 n 1 1 n 3 k 3 n 3 k 3 n 3 f 0 , 1 , k 3 n 3 1 + k 3 n 3 f 0 , 1 , k 3 1 n 3 1 f 0 , 1 , k 3 n 3 t 2 n 2 t 3 k 3 .
Then F n 1 , n 2 , n 3 4 ( t 1 , t 2 , t 3 ) is positive since f is convex and all the t i , i { 1 , 2 , 3 } , are positive. □

4. Order of Convergence

To provide information about the speed of convergence of multivariate Bernstein polynomials, we need the following intermediate result. This result also exists for one variable (see (Reference [1] Lemma 6.3.5)).
Lemma 4.
There is a constant C independent of n i such that, for all
( x 1 , , x m ) [ 0 , 1 ] m , i = 1 , , m ,
k 1 n 1 x 1 n 1 1 / 4 k m n m x m n m 1 / 4 i = 1 m n i k i x i k i ( 1 x i ) n i k i C i = 1 m n i 3 / 2 .
Proof of Lemma 4. 
We will sometimes use the square brackets to draw attention to computations of certain terms.
k 1 n 1 x 1 n 1 1 / 4 k m n m x m n m 1 / 4 i = 1 m n i k i x i k i ( 1 x i ) n i k i = k 2 n 2 x 2 n 2 1 / 4 k m n m x m n m 1 / 4 k 1 n 1 x 1 n 1 1 / 4 n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 × i = 2 m n i k i x i k i ( 1 x i ) n i k i , = k 1 n 1 x 1 n 1 1 / 4 n 1 k 1 x 1 k 1 ( 1 x 1 ) n 1 k 1 k 3 n 3 x 3 n 3 1 / 4 k m n m x m n m 1 / 4 k 2 n 2 x 2 n 2 1 / 4 n 2 k 2 x 2 k 2 ( 1 x 2 ) n 2 k 2 i = 3 m n i k i x i k i ( 1 x i ) n i k i = i = 1 2 k i n i x i n i 1 / 4 n i k i x i k i ( 1 x i ) n i k i k 3 n 3 x 3 n 3 1 / 4 k m n m x m n m 1 / 4 i = 3 m n i k i x i k i ( 1 x i ) n i k i .
By repeating successively the actions with k 3 , k 4 , , k m , we obtain
k 1 n 1 x 1 n 1 1 / 4 k m n m x m n m 1 / 4 i = 1 m n i k i x i k i ( 1 x i ) n i k i = i = 1 m k i n i x i n i 1 / 4 n i k i x i k i ( 1 x i ) n i k i .
Let us use the following result which can be found in References [1,21,22]: There is a constant C independent of a given positive integer n such that for all x in [ 0 , 1 ] ,
k n x n 1 / 4 n k x k ( 1 x ) n k C n 3 / 2 .
By applying this to any x i , i = { 1 , , m } , we have
k 1 n 1 x 1 n 1 1 / 4 k m n m x m n m 1 / 4 i = 1 m n i k i x i k i ( 1 x i ) n i k i
i = 1 m C i n i 3 / 2 with C i constants
= C i = 1 m n i 3 / 2 with C constant ( C = C 1 C m ) .
 □
Remark 2.
Lemma 4 is still valid if the summation is not on all the k i , 1 i m but on some of them. In this case, the upper bound will contain only terms involved in the summation.
The following Theorem will help us to conclude on the order of convergence of our Bernstein polynomials with several variables, see References [22,36].
Theorem 6.
Let f ( x 1 , , x m ) be bounded in [ 0 , 1 ] m and let ( a 1 , , a m ) be a point of [ 0 , 1 ] m at which f ( 2 ) ( 2 ) ( a 1 , , a m ) exists and is continuous. Then,
lim n j , 1 j m B n 1 , , n m ( f ; a 1 , , a m ) f ( a 1 , , a m ) 1 2 j = 1 m a j ( 1 a j ) n j 2 f x j 2 ( a 1 , , a m ) = 0 ,
and
lim N N B N , , N ( f ; a 1 , , a m ) f ( a 1 , , a m ) = 1 2 i = 1 m a i ( 1 a i ) 2 f x i 2 ( a 1 , , a m ) .
Proof of Theorem 6. 
Since f ( 2 ) ( 2 ) ( x 1 , , x m ) exists and is continuous at ( a 1 , , a m ) , then using Taylor expansion of order two
f ( x 1 , , x m ) = f ( a 1 , , a m ) + j = 1 m f x j ( a 1 , , a m ) ( x j a j ) + 1 2 j = 1 m 2 f x j 2 ( a 1 , , a m ) ( x j a j ) 2 + i , j = 1 , i j m 2 f x i x j ( a 1 , , a m ) ( x i a i ) ( x j a j ) + s ( x 1 , , x m ) j = 1 m ( x j a j ) 2 ,
where lim ( x 1 , , x m ) ( a 1 , , a m ) s ( x 1 , , x m ) = 0 . Set ( x 1 , , x m ) = k 1 n 1 , , k m n m , multiplying both sides by i = 1 m n i k i a i k i ( 1 a i ) n i k i and sum from k i = 0 to k i = n i , for all i = 1 , , m , we obtain
B n 1 , , n m ( f ; a 1 , , a m ) = f ( a 1 , , a m ) + j = 1 m f x j ( a 1 , , a m ) k 1 = 0 n 1 k m = 0 n m k j n j a j i = 1 m n i k i a i k i ( 1 a i ) n i k i + 1 2 j = 1 m 2 f x j 2 ( a 1 , , a m ) k 1 = 0 n 1 k m = 0 n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i + 1 2 i , j = 1 , i j m 2 f x i x j ( a 1 , , a m ) k 1 = 0 n 1 k m = 0 n m k i n i a i k j n j a j l = 1 m n l k l a l k l ( 1 a l ) n l k l + k 1 = 0 n 1 k m = 0 n m s k 1 n 1 , , k m n m j = 1 m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i = f ( a 1 , , a m ) + j = 1 m f x j ( a 1 , , a m ) k j = 0 n j k j n j a j n j k j a j k j 1 a j n j k j × k 1 = 0 , 1 j n 1 k m = 0 , m j n m i = 1 , i j m n i k i a i k i ( 1 a i ) n i k i + 1 2 j = 1 m 2 f x j 2 ( a 1 , , a m ) × k j = 0 n j k j n j a j 2 n j k j a j k j ( 1 a j ) n j k j k 1 = 0 , 1 j n 1 k m = 0 , m j n m i = 1 , i j m n i k i a i k i ( 1 a i ) n i k i + 1 2 i , j = 1 , i j m 2 f x i x j a 1 , , a m k i = 0 n i k i n i a i n i k i a i k i ( 1 a i ) n i k i k j = 0 n j k j n j a j n j k j × a j k j ( 1 a j ) n j k j k 1 = 0 , 1 i , j n 1 k m = 0 , m i , j n m l = 1 , l i , j m n l k l a l k l ( 1 a l ) n l k l + j = 1 m k 1 = 0 n 1 k m = 0 n m s k 1 n 1 , , k m n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i .
Let us use the following results from P. J. Davis [1] which are valid for all x [ 0 , 1 ] :
k = 0 n n k k n x x k ( 1 x ) n k = 0 ; k = 0 n n k k n x 2 x k ( 1 x ) n k = x ( 1 x ) n ; k = 0 n n k k n x 3 x k ( 1 x ) n k = x ( 1 x ) 2 x + 1 n 2 ; k = 0 n n k k n x 4 x k ( 1 x ) n k = x ( 1 x ) ( 3 n 6 ) x ( 1 x ) + 1 n 3 .
The second equation in system ( 135 ) can also be found in (Reference [33] Corollary 5.3). Hence, with Equations ( 31 ) and ( 135 ) , we have
B n 1 , , n m ( f ; a 1 , , a m ) = f ( a 1 , , a m ) + j = 1 m a j ( 1 a j ) 2 n j 2 f x j 2 ( a 1 , , a m ) + j = 1 m k 1 = 0 n 1 k m = 0 n m s k 1 n 1 , , k m n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i .
Let
S = j = 1 m k 1 = 0 n 1 k m = 0 n m s k 1 n 1 , , k m n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i
and ε > 0 . We can find n i , i = 1 , , m sufficiently large such that | x i a i | < 1 n i 1 / 4 implies | s ( x 1 , , x m ) | ε since lim ( x 1 , , x m ) ( a 1 , , a m ) s ( x 1 , , x m ) = 0 . Let us consider the set
E = i = 1 m { 0 , , n i } , and for j = 1 , , m , define the sets
Ω j = k j { 0 , , n j } : k j n j a j < 1 n j 1 / 4 and F = E ( Ω 1 × × Ω m ) .
Then F = k = 1 m F k , with F k = i = 1 m Ω i k [ α i k ] F : α i k { 0 , 1 } , i = 1 m α i k = k , where Ω i k [ α i k ] = Ω i if α i k = 0 Ω i c if α i k = 1 and Ω i c = k i { 0 , , n i } : k i n i x i 1 n i 1 / 4 . For A k = i = 1 m Ω i k [ α i k ] F k , k = 1 , , m , let us also define I A k = i { 1 , , m } : α i k = 1 for k = 1 , , m (that means card I A k = k 1 ). We have
| S | j = 1 m Ω 1 Ω m s k 1 n 1 , , k m n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i + F s k 1 n 1 , , k m n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i
ε j = 1 m k 1 = 0 n 1 k m = 0 n m k j n j a j 2 i = 1 m n i k i a i k i ( 1 a i ) n i k i + M l = 1 m A l F l i = 1 m n i k i a i k i ( 1 a i ) n i k i ) since F = l = 1 m F l where M = sup ( x 1 , , x m ) [ 0 , 1 ] m s ( x 1 , , x m ) j = 1 m ( x j a j ) 2 ε j = 1 m k j = 0 n j k j n j a j 2 n j k j a j k j ( 1 a j ) n j k j k 1 = 0 , 1 j n 1 k m = 0 , m j n m i = 1 , i j m n i k i a i k i ( 1 a i ) n i k i + M l = 1 m A l F l i I A l C i n i 3 / 2 ,
with C i constants (we bounded the summation on F using Equation ( 129 ) and the same procedure as we have used to show ( 55 ) in Theorem 2. Hence,
| S | ε j = 1 m k j = 0 n j k j n j a j 2 n j k j a j k j ( 1 a j ) n j k j + M l = 1 m A l F l i I A l C n 3 / 2 using ( 31 ) , a n d n = min { n j , 1 j m } n i , i = 1 , , m and C = max ( C i : i = 1 , , m ) ε j = 1 m a j ( 1 a j ) n j + M l = 1 m A l F l C l n 3 l / 2 using Equation ( 135 ) and card I A l = l ε n j = 1 m a j ( 1 a j ) + 1 n l = 1 m A l F l M × C l n 1 / 2 since n 3 l / 2 > n 3 / 2 .
For n 1 , , n m sufficiently large (that means n sufficiently large), we have l = 1 m A l F l M × C l n 1 / 2 < ε . This implies
| S | ε n 1 + j = 1 m a j ( 1 a j ) .
Therefore,
B n 1 , , n m ( f ; a 1 , , a m ) f ( a 1 , , a m ) j = 1 m a j ( 1 a j ) 2 n j 2 f x j 2 ( a 1 , , a m ) = | S | ε n 1 + j = 1 m a j ( 1 a j ) .
Since ε is arbitrary, we have
lim n j , 1 j m B n 1 , , n m ( f ; a 1 , , a m ) f ( a 1 , , a m ) j = 1 m a j ( 1 a j ) 2 n j 2 f x j 2 ( a 1 , , a m ) = 0 .
 □
Remark 3.
By taking n i = N , 1 i m , in Theorem 6 , we have
lim N N B N , , N ( f ; a 1 , , a m ) f ( a 1 , , a m ) = 1 2 i = 1 m a i ( 1 a i ) 2 f x i 2 ( a 1 , , a m ) .
From the previous Remark, one can deduce that the order of convergence of Bernstein polynomials with several variables for a bounded function f ( x 1 , , x m ) , ( x 1 , , x m ) [ 0 , 1 ] m , in some neighborhood of ( a 1 , , a m ) [ 0 , 1 ] m included in [ 0 , 1 ] m , is one, provided that f ( 2 ) ( a 1 , , a m ) exists. That is less than Taylor’s order of convergence.

5. Conclusions

Our concern here was the multivariate Bernstein polynomials given by Equation ( 6 ) which we recall here
B n 1 , , n m ( f ; x 1 , , x m ) = k 1 = 0 n 1 k m = 0 n m f k 1 n 1 , , k m n m i = 1 m n i k i x i k i ( 1 x i ) n i k i ,
and which is a generalization in several variables of the univariate Bernstein polynomials given by ( 1 )
B n ( f ; x ) = k = 0 n f k n n k x k ( 1 x ) n k .
We extended the beautiful properties of the above mentioned univariate Bernstein polynomials by generalizing them for functions of several real variables: uniform convergence, uniform convergence of the derivatives, order of convergence, monotonicity, fixed sign for the p-th derivative, and deduction of the upper and lower bounds of the Bernstein polynomials from those of the corresponding functions. To achieve that, we first recalled not only these properties for univariate Bernstein polynomials, but also the intermediate properties used to prove them (see Equations ( 11 ) ( 16 ) ). Then, we extended these intermediate properties for the univariate case to obtain the generalised intermediate properties for the multivariate case given by ( 6 ) (see ( 20 ) ( 26 ) ). Finally, using those generalized intermediate properties in combination with some others tools related to functions with several variables, we got the generalization of these beautiful properties of the univariate Bernstein polynomials to multivariate Bernstein polynomials given by Equation ( 6 ) .
More precisely, our contributions in this paper related to this form of multivariate Bernstein polynomials can be summarized as follows:
  • deriving the expression of these multivariate Bernstein polynomials in terms of differences (similar to the case of one variable in References [1,15]);
  • establishing the pointwise and uniform convergence of these Bernstein polynomials with several variables (a generalization of the case of two variables treated in (Reference [28] Lemma 1));
  • deriving the expression of the p 1 , , p m -th derivative of these multivariate Bernstein polynomials of order n 1 + p 1 , , n m + p m -th in terms of p 1 , , p m -th differences of the corresponding functions (similar to the case of one variable in References [1,15]);
  • proving the monotonicity, the fixed sign of the p 1 , , p m -th derivative, the convexity, and the upper and lower bounds of these multivariate Bernstein polynomials from those of the corresponding functions (similar to the case of one variable in References [1,31]);
  • establishing the pointwise and uniform convergence of the p 1 , , p m -th derivative of these multivariate Bernstein polynomials (a generalization of the case of two variables treated in (Reference [35] Theorem of Kingsley) and (Reference [28] Theorem 1));
  • generalizing the inequality related to the decreasing of the sequence of Bernstein’s polynomials with one variable (see References [1,31]) when f is a convex function;
  • proving that the order of convergence of these Bernstein polynomials with several variables is one.
With these new results, we have proved that the conclusion made by P.J. Davis in Reference [1] for Bernstein polynomials of one variable is true for the case of several variables, that is, the Bernstein polynomials of several variables yield smooth approximants, but their convergence is very slow.

Author Contributions

The conceptualization and the supervision were done by the first author (M.F.). The formal analysis and the investigation were done by the second author (M.M.W.). All authors have read and agreed to the published version of the manuscript.

Funding

The work of these two authors was initially funded by the African Institute for Mathematical Sciences-Cameroon, during the academic year 2013/2014.

Acknowledgments

The authors would like to acknowledge the support of the two institutions where this work has been carried out: The Mathematics Institute of the University of Kassel in Germany and the African Institute for Mathematical Sciences in Cameroon.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relation-ships that could have appeared to influence the work reported in this paper.

References

  1. Davis, P.J. Interpolation and Approximation; Dover Publications: New York, NY, USA, 1975. [Google Scholar]
  2. Kim, T.; Jang, L.-C.; Yi, H. A note on the modified q-Bernstein polynomials. Discret. Dyn. Nat. Soc. 2010. [Google Scholar] [CrossRef] [Green Version]
  3. Lu, W.-M.; Zhang, L. Pointwise approximation theorems for combinations of Bernstein polynomials with inner singularities. Appl. Math. 2010, 2. [Google Scholar] [CrossRef] [Green Version]
  4. Simsek, Y.; Acikgoz, M. A new generating function of (q-)Bernstein-type polynomials and their interpolation function. Abstr. Appl. Anal. 2010. [Google Scholar] [CrossRef] [Green Version]
  5. Xie, L.S. Pointwise approximation theorems for combinations and derivatives of Bernstein polynomials. Acta Math. Sin. 2005, 21, 1241–1248. [Google Scholar] [CrossRef]
  6. Xie, L. Pointwise simultaneous approximation by combination of Bernstein operators. J. Approx. Theory 2005, 137, 1–21. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, A.; Veretennikova, E.V. On partial derivatives of multivariate Bernstein polynomials. Sib. Adv. Math. 2016, 26, 294–305. [Google Scholar]
  8. Gal, S.G.; Anastassiou, G.A. Shape-Preserving Approximation by Real and Complex Polynomials; Birkhäuser: Boston, MA, USA, 2010. [Google Scholar]
  9. Peña, J.M. Shape Preserving Representation in Computer-Aided Geometric Design; Nova Science Publishers, Inc.: New York, NY, USA, 1999. [Google Scholar]
  10. Occorsio, D. Some new properties of generalized Bernstein polynomials. Stud. Univ. Babes-Bolyai Math. 2011, 56, 147–160. [Google Scholar]
  11. Costabile, F.; Gualtieri, M.I.; Serra, S. Asymptotic expansion and extrapolation for Bernstein Polynomials with applications. BIT Numer. Math. 1996, 36, 676–687. [Google Scholar] [CrossRef]
  12. Jafarian, A.; Nia, S.M.; Golmankhaneh, A.K.; Baleanu, D. On Bernstein Polynomials Method to the System of Abel Integral Equations. Abstr. Appl. Anal. 2014. [Google Scholar] [CrossRef]
  13. Tucker, D.H. A note on Bernstein polynomial type approximations. Proc. Am. Math. Soc. 1967, 18, 492–494. [Google Scholar] [CrossRef]
  14. Doha, E.H.; Bhrawy, A.H.; Saker, M.A. Integrals of Bernstein polynomials: An application for the solution of high even-order differential equations. Appl. Math. Lett. 2011, 24, 559–565. [Google Scholar] [CrossRef] [Green Version]
  15. Phillips, G.M. Interpolation and Approximation by Polynomials; Springer: New York, NY, USA, 2003. [Google Scholar]
  16. Gzyl, H.; Palacios, J.L. On the Approximation Properties of Bernstein Polynomials via Probabilistic Tools. Boletín de la Asociación Matemática Venez. 2003, 10, 5–13. [Google Scholar]
  17. Alshbool, M.H.T.; Hashim, I. Multistage Bernstein polynomials for the solutions of the Fractional Order Stiff Systems. J. King Saud Univ. Sci. 2016, 28, 280–285. [Google Scholar] [CrossRef] [Green Version]
  18. Kadison, R.V.; Liu, Z. Bernstein Polynomials and Approximation. 2018. Available online: https://www.math.upenn.edu/kadison/bernstein.pdf (accessed on 30 April 2020).
  19. Dalmo, R. Matrix factorization of multivariate Bernstein polynomials. Int. J. Pure Appl. Math. 2015, 103, 749–780. [Google Scholar] [CrossRef] [Green Version]
  20. Dinghas, A. Über einige Identitäten vom Bernsteinschen Typus. Det Koneglige Norske Videnskabers Selskab. 1951, 24, 96–97. [Google Scholar]
  21. Duchoň, M. A generalized Bernstein Approximation theorem. Tatra Mt. Math. Publ. 2011, 49, 99–109. [Google Scholar] [CrossRef] [Green Version]
  22. Lorentz, G.G. Bernstein Polynomials; American Mathematical Soc.: New York, NY, USA, 2013. [Google Scholar]
  23. Bayad, A.; Kim, T.; Rim, S.-H. Bernstein Polynomials on Simplex. 2011. Available online: https://arxiv.org/pdf/1106.2482.pdf (accessed on 30 April 2020).
  24. Campiti, M.; Raşa, I. Extrapolation Properties of Multivariate Bernstein Polynomials. Mediterr. J. Math. 2019, 16, 109. [Google Scholar] [CrossRef]
  25. Derriennic, M.-M. On multivariate Approximation by Bernstein-Type polynomials. J. Approx. Theory 1985, 45, 155–166. [Google Scholar] [CrossRef] [Green Version]
  26. Ding, C.; Cao, F. K-Functionals and multivariate Bernstein polynomials. J. Approx. Theory 2008, 155, 125–135. [Google Scholar] [CrossRef] [Green Version]
  27. Jetter, K.; Stöckler, J. An identity for multivariate Bernstein polynomials. Comput. Aided Geom. Des. 2003, 20, 563–577. [Google Scholar] [CrossRef]
  28. Martinez, F.L. Some properties of two-dimensional Bernstein polynomials. J. Approx. Theory 1989, 59, 300–306. [Google Scholar] [CrossRef] [Green Version]
  29. Heitzinger, C. Simulation and Inverse Modeling of Semiconductor Manufacturing Processes; Technische Universität Wien: Vienna, Australia, 2002. [Google Scholar]
  30. Ditzian, Z. Direct estimate for Bernstein polynomials. J. Approx. Theory 1994, 79, 165–166. [Google Scholar] [CrossRef] [Green Version]
  31. Sauer, T. Multivariate Bernstein polynomials and convexity. Comput. Aided Geom. Des. 1991, 8, 465–478. [Google Scholar] [CrossRef]
  32. Popoviciu, T. Sur l’approximation des fonctions convexes d’ordre supérieur. Bull. Mathématique De La Société Roumaine Des Sci. 1941, 43, 85–141. [Google Scholar]
  33. Fellhauer, A. Approximation of Smooth Functions Using Bernstein Polynomials in Multiple Variables. 2016. Available online: https://arxiv.org/pdf/1609.01940.pdf (accessed on 30 April 2020).
  34. Acikgoz, M.; Araci, S. A new approach to modified q-Bernstein polynomials for functions of two variables with their generating and interpolation functions. Turk. J. Anal. Number Theory 2012. [Google Scholar] [CrossRef] [Green Version]
  35. Kingsley, E.H. Bernstein polynomials for functions of two variables of class C(k). Proc. Am Math. Soc. 1951, 2, 64–71. [Google Scholar]
  36. Voronowskaja, E. Détermination de la forme asymptotique d’approximation des fonctions par les polynomes de M. Bernstein. C. R. Acad. Sci. URSS 1932, 79, 79–85. [Google Scholar]
  37. Agrawal, P.N. Inverse theorem in simultaneous approximation by Micchelli combination of Bernstein polynomials. Demonstr. Math. 1994, 27, 501–510. [Google Scholar] [CrossRef]
  38. Bernstein, S.N. Complément à l’article de E. Voronowskaja: The asymptotic behavior of the approximation of functions by their Bernstein polynomials. C. R. Acad. Sci. URSS 1932, 86–92. [Google Scholar]
  39. Butzer, P.L. Linear combinations of Bernstein polynomials. Can. J. Math. 1953, 5, 559–567. [Google Scholar] [CrossRef]
  40. Guo, S.; Yue, S.; Li, C.; Yang, G.; Sun, Y. A Pointwise approximation theorem for linear combinations of Bernstein polynomials. Abstr. Appl. Anal. 1996, 1. [Google Scholar] [CrossRef] [Green Version]
  41. WooD, B. Uniform approximation by linear combinations of Bernstein-Type polynomials. J. Approx. Theory 1984, 41, 51–55. [Google Scholar] [CrossRef] [Green Version]
  42. Acikgoz, M.; Sen, E.; Araci, S. q-Bernstein-Type polynomials for functions of two variables with their generating and interpolation functions. Turk. J. Anal. Number Theory 2013, 1, 36–42. [Google Scholar] [CrossRef] [Green Version]
  43. Oruç, H.; Tuncer, N. On the q-Bernstein polynomials. In Algorithms for Approximation IV, The Proceedings of the 2001 International Symposium, Queensgate, Huddersfield, UK, 16–20 July 2001; University of Huddersfield: Huddersfield, UK, 2002. [Google Scholar]
  44. Phillips, G.M. On generalized Bernstein polynomials. Numer. Anal. 1996, 263–269. [Google Scholar] [CrossRef]
  45. Vishwakarma, D.K.; Kumar, A.A.; Kumar, A. Multivariate q-Bernstein-Schurer-Kantorovich operators. J. Math. Syst. Sci. 2016, 6, 234–241. [Google Scholar] [CrossRef] [Green Version]
  46. Bernstein, S.N. Démonstration du théorème de Weierstrass fondée sur le calcul de probabilités. Commun. Kharkov Math. Soc. 1913, 13, 1–2. [Google Scholar]
  47. Wouodjié, M.M. On Multivariate Bernstein Polynomials. Master’s Thesis, African Institute for Mathematical Sciences, Limbé, Cameroon, 2014. Available online: http://www.mathematik.uni-kassel.de/~merlin/Merlin_Master_Thesis_AIMS_Cameroon.pdf (accessed on 30 April 2020).
  48. Reimer, M. Multivariate Polynomial Approximation; Birkhäuserc: Basel, Switzerland, 2012. [Google Scholar]

Share and Cite

MDPI and ACS Style

Foupouagnigni, M.; Mouafo Wouodjié, M. On Multivariate Bernstein Polynomials. Mathematics 2020, 8, 1397. https://doi.org/10.3390/math8091397

AMA Style

Foupouagnigni M, Mouafo Wouodjié M. On Multivariate Bernstein Polynomials. Mathematics. 2020; 8(9):1397. https://doi.org/10.3390/math8091397

Chicago/Turabian Style

Foupouagnigni, Mama, and Merlin Mouafo Wouodjié. 2020. "On Multivariate Bernstein Polynomials" Mathematics 8, no. 9: 1397. https://doi.org/10.3390/math8091397

APA Style

Foupouagnigni, M., & Mouafo Wouodjié, M. (2020). On Multivariate Bernstein Polynomials. Mathematics, 8(9), 1397. https://doi.org/10.3390/math8091397

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop