Return to Article Details On a generalization of the Stancu-Schurer operator of higher order

On a generalization of the Stancu-Schurer operator of higher orderThanks: †Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, P.O. Box 66-1, Cluj-Napoca, RomaniaThanks: ‡Faculty of Mathematics and Computer Science, Babeş-Bolyai University, Str. M. Kogălniceanu 1, 400084 Cluj-Napoca, Romania, email: grigoriciuc@ictp.acad.ro andra.malina@ictp.acad.ro.

Eduard Ştefan Grigoriciuc †,‡ and Andra Malina †,‡
Dedicated to Professor Teodora Cătinaş on the occasion of her anniversary.
Date: May 05, 2026; accepted: May 20, 2026; published online: June 30, 2026.
Abstract.

In this paper, we introduce a generalization of the higher-order Stancu–Schurer operator. Starting from a particular form of the operator recently introduced by the authors, we extend the convex combination of two terms appearing in its expression to a convex combination of m functions, where m∈ℕ with m≥1. For these generalized operators, we study classical properties such as linearity, positivity, monotonicity, moments, and certain convexity properties. We conclude the paper with some remarks on a nonlinear extension with data-dependent weights.

Key words and phrases: 
Stancu-Schurer operator, generalized Stancu operator of higher order, convex combination, linear positive operator
1991 Mathematics Subject Classification
41A35; 41A36; 47A58

1. Introduction

For f∈C⁡[0,1] and n∈ℕ, Bernstein introduced in [3] the linear and positive operators

(1) ℬn⁢f⁢(x)=∑k=0nbn,k⁢(x)⁢f⁢(kn),

with the basis polynomials bn,k given by

(2) bn,k⁢(x)=(nk)⁢xk⁢(1−x)n−k,

for k∈{0,…,n} and x∈[0,1]. Another important family of operators was introduced by Stancu in [21] using a probabilistic approach. We denote these positive and linear operators by

(3) Sn⁢f⁢(x)=∑k=0n−rbn−r,k⁢(x)⁢[(1−x)⁢f⁢(kn)+x⁢f⁢(k+rn)],

where the polynomials bn−r,k are given by (2), f∈C⁡[0,1], x∈[0,1] and r is a non-negative integer satisfying n>2⁢r, for every n∈ℕ.

Let m∈ℕ with m≥1. Also, let n,r1,…,rm,R∈ℕ be such that n>2⁢R, where R=r1+…+rm. For every function f∈C⁡[0,1], the authors introduced in [17] the generalized Stancu operator of order m, given by

(4) Snm⁢f⁢(x)=∑k=0n−Rbn−R,k⁢(x)⁢[(1−∑j=1mxj)⁢f⁢(kn)+∑j=1mxj⁢f⁢(k+rjn)],

where the polynomials bn−R,k are given by (2), for all x∈Dm, where

Dm={x∈[0,1]: 1−∑j=1mxj≥0}⊆[0,1].

It is clear that for m=1 the operator Sn1 reduces to the classical Stancu operator Sn given by (3). Moreover, if R=0, then the operator Snm reduces to the Bernstein operator given by (1).

Remark 1.

The definition of the domain Dm is important for obtaining the positivity of the considered operators. Although Dm=[0,1] if and only if m=1, it can nevertheless be easily shown that Dm→[0,1/2], when m→∞. Hence, the domain of definition is not negligible for an arbitrary value of the order m∈ℕ. Concrete examples of the domains Dm can be obtained, as follows:

D1=[0,1],D2=[0,5−12],D3=[0,α],

where α=13⁢(17+3⁢333+17−3⁢333−1). In general,

Dm=[0,αm],

where αm∈[0,1] is the root of the equation x+x2+…+xm−1=0. Moreover, if m→∞ and x∈[0,1), then

limm→∞(x+x2+…+xm)−1=x1−x−1=2⁢x−11−x

and the root of the considered equation is x=12, i.e. D∞=[0,1/2].

To generalize the previous operator, a small modification can be made to its definition. Similarly to what has been done in [19] (see also [12], [18]), we consider p∈ℕ with p≥0 and C⁡[0,1+p] the space of continuous functions on the interval [0,1+p]. Then we obtain the generalized Stancu-Schurer operator of order m, given by

Sn,pm⁢f⁢(x)=∑k=0n+p−Rbn+p−R,k⁢(x)⁢[(1−∑j=1mxj)⁢f⁢(kn)+∑j=1mxj⁢f⁢(k+rjn)],

for all x∈Dm, where the polynomials bn+p−R,k are given by (2) and f∈C⁡[0,1+p].

An important step that can be undertaken is to replace the Bernstein polynomials bn+p−R,k with the generalized polynomials introduced by Cheney and Sharma in [14], constructed on the basis of Jensen’s formulas. This approach was considered by the authors in [17] in the context of generalized higher-order Stancu–Schurer operators or by T. Bostanci and G. Başcanbaz-Tunca in [4] for a Stancu-type extension of the Cheney-Sharma operator. Such a study falls within the research direction established by T. Cătinaş in the works [5], [6], [8], [9] and [10], where Cheney–Sharma type operators are investigated on various domains. This direction was also followed by the authors in the works [16] and [17]. Another direction that was considered by T. Cătinaş in [7] consists of a Stancu-type extension for another linear and positive operator, namely the Campiti-Metafune operator. Recently, some of the previously mentioned operators have also been studied in a complex setting (see [11], [12], [13], [15]).

2. The generalized Stancu-Schurer operator of higher order

Let m∈ℕ with m≥1. Also, let n,r1,…,rm,R∈ℕ be such that n>2⁢R, where R=r1+…+rm. For simplicity, let us denote λ=n+p−R. Then, for every function f∈C⁡[0,1+p], we denote by

(5) ℳλm⁢f⁢(x)=∑k=0λbλ,k⁢(x)⁢[Aλ,k⁢(x)⁢f⁢(kn)+∑j=1mBλ,k,j⁢(x)⁢f⁢(k+rjn)],

where the polynomials bλ,k are given by (2) and the functions Aλ,k, Bλ,k,j have the properties

(6) Aλ,k⁢(x)∈[0,1],Bλ,k,j⁢(x)≥0

for every j∈{1,…,m}, k∈{0,…,λ} and x∈[0,1]. Moreover, for a fixed k∈{0,…,λ}, we impose the condition

(7) Aλ,k⁢(x)+∑j=1mBλ,k,j⁢(x)=1,

for all x∈[0,1]. We call the operator ℳλm given by (5) the generalized Stancu-Schurer operator of higher order (in particular, of order m). It is clear that ℳλm is a generalization of the operator Sn,pm studied by the authors in [17]. Indeed, if

Bλ,k,j⁢(x)=xjandAλ,k⁢(x)=1−∑j=1mxj,

for all x∈Dm, then the assumptions for Aλ,k and Bλ,k,j presented above are fulfilled and ℳλm reduces to Sn,pm.

Remark 2.

Let us consider

(8) tλ,k,m⁢(x)=∑j=1mBλ,k,j⁢(x)

such that tλ,k,m⁢(x)∈[0,1], for all x∈[0,1]. Then

Aλ,k⁢(x)=1−tλ,k,m⁢(x)∈[0,1].

If tλ,k,m⁢(x)>0, we can define

(9) ωλ,k,j⁢(x)=Bλ,k,j⁢(x)tλ,k,m⁢(x)

with the property that ∑j=1mωλ,k,j⁢(x)=1. Then the operator ℳλm becomes

ℳλm⁢f⁢(x) =∑k=0λbλ,k(x)[(1−tλ,k,m(x))f(kn)+
+tλ,k,m(x)∑j=1mωλ,k,j(x)f(k+rjn)],

where the polynomials bλ,k are given by (2) and the functions tλ,k,m and ωλ,k,j are given by relation (8), respectively by (9). We notice here that the function tλ,k,m plays the role of a mixing parameter indicating how much data is assigned to shifted nodes, while the function ωλ,k,j can be seen as a probability function that measures how the data is distributed among the shifts.

Remark 3.

According to the previous form of the operator ℳλm, it is not difficult to observe that for m=1 we obtain the operator

ℳλ1⁢f⁢(x)=∑k=0λbλ,k⁢(x)⁢[(1−tλ,k,1⁢(x))⁢f⁢(kn)+tλ,k,1⁢(x)⁢f⁢(k+rn)]

that is a generalization of the Stancu operator defined by (3). The operator becomes a two point convex combination between the classical node k/n and the shifted node (k+r)/n with a mixing parameter tλ,k,1⁢(x)∈[0,1].

In particular, if tλ,k,1⁢(x)=x, for all x∈[0,1], then ℳλ1 reduces to the classical Stancu-Schurer operator (see [12], [19]). Moreover, if p=0, then ℳn−r1 is the Stancu operator given by (3).

3. Properties of the operator ℳλm

In this section, we consider some properties of the operator ℳλm. We mention that all the assumptions on bλ,k, Aλ,k and Bλ,k,j, respectively for the parameters considered are the same as in the previous section. In the interest of brevity, we do not repeat them here; for details, we refer the reader to Section 2.

Proposition 4.

The operator ℳλm is linear and positive on C⁡[0,1+p].

Proof.

In order to show that the operator ℳλm is linear, let us consider f,g∈C⁡[0,1+p] and α,β∈ℝ. Then

ℳλm⁢(α⁢f+β⁢g)⁢(x) =∑k=0λbλ,k(x)[Aλ,k(x)(αf+βg)(kn)
+∑j=1mBλ,k,j(x)(αf+βg)(k+rjn)]
=∑k=0λbλ,k(x)[Aλ,k(x)(αf(kn)+βg(kn))
+∑j=1mBλ,k,j(x)(αf(k+rjn)+βg(k+rjn))]
=∑k=0λbλ,k(x)[αAλ,k(x)f(kn)+∑j=1mαBλ,k,j(x)f(k+rjn)
+βAλ,k(x)g(kn)+∑j=1mβBλ,k,j(x)g(k+rjn)]
=α⁢∑k=0λbλ,k⁢(x)⁢[Aλ,k⁢(x)⁢f⁢(kn)+∑j=1mBλ,k,j⁢(x)⁢f⁢(k+rjn)]
+β∑k=0λbλ,k(x)[Aλ,k(x)g(kn)+∑j=1mBλ,k,j(x)g(k+rjn)]
=α⁢ℳλm⁢f⁢(x)+β⁢ℳλm⁢g⁢(x)
=(α⁢ℳλm⁢f+β⁢ℳλm⁢g)⁢(x),

for all x∈Dm. Moreover, since bλ,k, Aλ,k and Bλ,k,j are positive functions on [0,1], it follows that for any positive function f∈C⁡[0,1+p] we obtain that

ℳλm⁢f⁢(x)≥0,x∈Dm,

and this completes the proof. ∎

A direct consequence of the previous result is the monotonicity of the operator ℳλm.

Corollary 5.

Let f,g∈C⁡[0,1+p] be such that f⁡(x)≤g⁡(x), for all x∈[0,1+p]. Then

ℳλm⁢f⁢(x)≤ℳλm⁢g⁢(x),

for all x∈Dm.

Proof.

Let f,g∈C⁡[0,1+p] be two functions such that f⁡(x)≤g⁡(x), for all x∈[0,1+p]. Then we consider h∈C⁡[0,1+p] such that h⁡(x)=g⁡(x)−f⁡(x), for all x∈[0,1+p]. It is clear that h⁡(x)≥0 on [0,1+p] and then ℳλm⁢h⁢(x)≥0, for all x∈Dm according to Proposition 4. Moreover, based on the linearity of the operator ℳλm we deduce that

0≤ℳλm⁢h⁢(x)=ℳλm⁢(g−f)⁢(x)=ℳλm⁢g⁢(x)−ℳλm⁢f⁢(x)

and then

ℳλm⁢f⁢(x)≤ℳλm⁢g⁢(x),

for all x∈Dm and this completes the proof. ∎

Another important result concerning the operator ℳλm is the study of its moments, in particular the preservation (or approximation) of the test functions e0, e1 and e2. This result is presented in a broader setting and extends in a natural way the result obtained by the authors in [17, Theorems 3.1 and 3.2].

Theorem 6.

Let ej⁢(x)=xj, for every j∈{0,1,2} and x∈Dm. Then

(10) ℳλm⁢(e0,x)=1

and

(11) ℳλm⁢(e1,x)=λn⁢x+1n⁢∑k=0λbλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x),

for all x∈Dm. Moreover,

(12) ℳλm⁢(e2,x) =λ⁡(λ−1)n2⁢x2+λn2⁢x+2n2⁢∑k=0λk⁢bλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
+1n2∑k=0λbλ,k(x)∑j=1mrj2Bλ,k,j(x),

for all x∈Dm.

Proof.

Using property (7) and the well-known expressions for the moments of the Bernstein polynomials (1), we obtain

ℳλm⁢(e0,x)=∑k=0λbλ,k⁢(x)⁢(Aλ,k⁢(x)+∑j=1mBλ,k,j⁢(x))=∑k=0λbλ,k⁢(x)=1,

for all x∈Dm. Moreover, simple computations show that

ℳλm⁢(e1,x) =∑k=0λbλ,k⁢(x)⁢(Aλ,k⁢(x)⁢kn+∑j=1mBλ,k,j⁢(x)⁢k+rjn)
=∑k=0λbλ,k⁢(x)⁢kn+∑k=0λbλ,k⁢(x)⁢∑j=1mrjn⁢Bλ,k,j⁢(x)
=λn⁢x+1n⁢∑k=0λbλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x),

for all x∈Dm. Finally, we obtain that

ℳλm⁢(e2,x) =∑k=0λbλ,k⁢(x)⁢(Aλ,k⁢(x)⁢k2n2+∑j=1mBλ,k,j⁢(x)⁢(k+rjn)2)
=∑k=0λbλ,k⁢(x)⁢k2n2+2n2⁢∑k=0λk⁢bλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
+1n2∑k=0λbλ,k(x)∑j=1mrj2Bλ,k,j(x)
=λ2n2⁢∑k=0λk2λ2⁢bλ,k⁢(x)+2n2⁢∑k=0λk⁢bλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
+1n2∑k=0λbλ,k(x)∑j=1mrj2Bλ,k,j(x)
=λ2n2⁢(x2+x⁡(1−x)λ)+2n2⁢∑k=0λk⁢bλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
+1n2∑k=0λbλ,k(x)∑j=1mrj2Bλ,k,j(x)
=λ⁡(λ−1)n2⁢x2+λn2⁢x+2n2⁢∑k=0λk⁢bλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
+1n2∑k=0λbλ,k(x)∑j=1mrj2Bλ,k,j(x),

for all x∈Dm and this completes the proof. ∎

Remark 7.

In particular, if Bλ,k,j⁢(x)=xj, for all x∈[0,1] and j∈{1,…,m}, then Theorem 6 reduces to [17, Theorems 3.1 and 3.2]. Note that the remaining cases can be treated similarly, given the very general form of the proof of the previous result.

Lemma 8.

For all x∈Dm and any m∈ℕ with m≥1, the following relations hold:

  1. (1)

    limn→∞ℳλm⁢(e1,x)=x;

  2. (2)

    ℳλm⁢(e1,x)∈[0,λ+rm⁢a⁢xn]⊆[0,1+p], where rm⁢a⁢x=max⁡{r1,…,rm}.

In particular, if m=1, then ℳλm⁢(e1,x)∈[0,n+pn]. Moreover, if p=0, then ℳλm⁢(e1,x)∈[0,1].

Proof.

In view of relation (11) we deduce immediately that

limn→∞ℳλm⁢(e1,x)=x,

for all x∈Dm. On the other hand, according to the assumption imposed on Bλ,k,j in Section 2, it is not difficult to observe that

0≤∑j=1mrj⁢Bλ,k,j⁢(x)≤rm⁢a⁢x⁢∑j=1mBλ,k,j⁢(x)=rm⁢a⁢x,

for all x∈[0,1], where rm⁢a⁢x=max⁡{r1,…,rm}. Since ℳλm⁢(e1,x) is given by relation (11), we obtain that

0≤λn⁢x≤ℳλm⁢(e1,x)≤λ⁢x+rm⁢a⁢xn≤λ+rm⁢a⁢xn,

for all x∈Dm. Moreover, it is clear that

rm⁢a⁢x=max⁡{r1,…,rm}≤R=r1+…+rm,

and then

λ+rm⁢a⁢xn=n+p−R+rm⁢a⁢xn≤n+pn=1+pn≤1+p,

for all n,p∈ℕ with p≥0 and n>2⁢R. Finally, the particular cases presented are obtained by simple substitutions of the parameters in the derived relations and this completes the proof. ∎

Another important result concerning the operator ℳλm is presented in the following proposition. It establishes that the operator is bounded in the space of continuous functions and, moreover, it preserves the uniform norm in the sense that it does not exceed the norm of the function on which it acts. For the particular case presented in Remark 7, one may consult [17, Proposition 3.3].

Proposition 9.

Let f∈C⁡[0,1+p]. Then

‖ℳλm⁢f‖C⁡(Dm)≤‖f‖C⁡[0,1+p].
Proof.

Let x∈Dm and f∈C⁡[0,1+p]. Using properties (6) and (7), we have that

|ℳλm⁢f⁢(x)| =|∑k=0λbλ,k⁢(x)⁢[Aλ,k⁢(x)⁢f⁢(kn)+∑j=1mBλ,k,j⁢(x)⁢f⁢(k+rjn)]|
≤∑k=0λbλ,k⁢(x)⁢|Aλ,k⁢(x)⁢f⁢(kn)+∑j=1mBλ,k,j⁢(x)⁢f⁢(k+rjn)|
≤∑k=0λbλ,k⁢(x)⁢[Aλ,k⁢(x)⁢|f⁡(kn)|+∑j=1mBλ,k,j⁢(x)⁢|f⁡(k+rjn)|]
≤‖f‖C⁡[0,1+p]⁢∑k=0λbλ,k⁢(x)⁢(Aλ,k⁢(x)+∑j=1mBλ,k,j⁢(x))
=‖f‖C⁡[0,1+p]⁢∑k=0λbλ,k⁢(x)
=‖f‖C⁡[0,1+p].

∎

Among the properties of interest in the study of such operators is the convexity property. Below, we present a convexity result which extends the one obtained by the authors in [17, Subsection 3.2]. For more details, one may consult [2], [16] and [20].

Proposition 10.

If f is a convex function on [0,1+p], then

ℳλm⁢f⁢(x)≥f⁡(ℳλm⁢(e1,x)),

for all x∈Dm.

Proof.

In view of relation (5), we know that

ℳλm⁢f⁢(x)=∑k=0λbλ,k⁢(x)⁢[Aλ,k⁢(x)⁢f⁢(kn)+∑j=1mBλ,k,j⁢(x)⁢f⁢(k+rjn)],

for all x∈Dm. For simplicity, let us denote by yi=k+rin, for all i∈{0,…,m} with r0=0 and by αk,0=Aλ,k⁢(x)≥0, respectively by αk,i=Bλ,k,i⁢(x), for all i∈{1,…,m}. According to the assumptions imposed in Section 2, we know that for a fixed k∈{0,…,λ} the following relation holds:

αk,0=1−∑i=1mαk,i∈[0,1].

Then the operator ℳλm can be written as

ℳλm⁢f⁢(x)=∑k=0λbλ,k⁢(x)⁢[αk,0⁢f⁢(y0)+∑i=1mαk,i⁢f⁢(yi)]=∑k=0λbλ,k⁢(x)⁢∑i=0mαk,i⁢f⁢(yi),

for all x∈Dm. Since f is convex on [0,1+p], it follows that

ℳλm⁢f⁢(x) =∑k=0λbλ,k⁢(x)⁢∑i=0mαk,i⁢f⁢(yi)
≥∑k=0λbλ,k⁢(x)⁢f⁢(∑i=0mαk,i⁢yi)
≥f⁡(∑k=0λbλ,k⁢(x)⁢∑i=0mαk,i⁢yi)
=f⁡(∑k=0λbλ,k⁢(x)⁢(αk,0⁢y0+∑i=1mαk,i⁢yi))
=f⁡(∑k=0λbλ,k⁢(x)⁢(Aλ,k⁢(x)⁢kn+∑j=1mBλ,k,j⁢(x)⁢k+rjn))
=f⁡(ℳλm⁢(e1,x)),

for all x∈Dm. Hence,

ℳλm⁢f⁢(x)≥f⁡(ℳλm⁢(e1,x)),

for all x∈Dm and this completes the proof. ∎

The previous result can be further refined by imposing a simple condition on the parameters under consideration, as follows:

Lemma 11.

If p≥R, then ℳλm⁢(e1,x)≥x, for all x∈Dm.

Proof.

In view of Theorem 6 we know that

ℳλm⁢(e1,x)=λn⁢x+1n⁢∑k=0λbλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x),

for all x∈Dm. Then

ℳλm⁢(e1,x)−x =λ−nn⁢x+1n⁢∑k=0λbλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
=p−Rn⁢x+1n⁢∑k=0λbλ,k⁢(x)⁢∑j=1mrj⁢Bλ,k,j⁢(x)
≥p−Rn⁢x,

since bλ,k and Bλ,k,j are positive functions on [0,1] and n,r1,…,rm∈ℕ. Hence, if p≥R, then

ℳλm⁢(e1,x)−x≥0,

for all x∈Dm and this completes the proof. ∎

We mention here that in the particular case Bλ,k,j⁢(x)=xj, for all x∈[0,1] we obtain the result proved by authors in [17, Lemma 3.1]. It is important to note that the condition imposed here is stronger than the one considered by the authors in [17]. However, in the particular case Bλ,k,j⁢(x)=xj, for all x∈[0,1] this condition can be relaxed, as in [17].

Theorem 12.

If p≥R and f is an increasing convex function on [0,1+p], then

(13) ℳλm⁢f⁢(x)≥f⁡(x),

for all x∈Dm.

Proof.

In view of Proposition 10 we know that if f is convex on [0,1+p], then

(14) ℳλm⁢f⁢(x)≥f⁡(ℳλm⁢(e1,x)),

for all x∈Dm. Moreover, since f is increasing on [0,1+p], it follows that for any two points u,v∈[0,1+p] with u≤v we have that f⁡(u)≤f⁡(v). In particular, for

u=x∈[0,1]⊆[0,1+p]

and

v=ℳλm⁢(e1,x)∈[0,λ+rmaxn]⊆[0,1+p]

we know (see Lemmas 8 and 11) that

u=x≤ℳλm⁢(e1,x)=v,

for all x∈Dm. Hence,

(15) f⁡(x)=f⁡(u)≤f⁡(v)=f⁡(ℳλm⁢(e1,x)),

for all x∈Dm. In view of relations (14) and (15) we deduce that

ℳλm⁢f⁢(x)≥f⁡(x),

for all x∈Dm and this completes the proof. ∎

4. Particular examples for the functions Aλ,k and Bλ,k,j

In [17] the authors obtained all the results presented above for the particular case Bλ,k,j⁢(x)=xj, for all x∈[0,1] and j∈{1,…,m}. In this section, we present other particular examples for the functions Aλ,k and Bλ,k,j that can be of interest for our study.

Example 13.

Let Bλ,k,j:[0,1]→[0,1] be given by

Bλ,k,j⁢(x)=cj⁢tλ,k,m⁢(x),

where

tλ,k,m⁢(x)=4⁢x⁢(1−x)⁢k⁡(λ−k)λ2,

for all x∈[0,1] and cj≥0 such that ∑j=1mcj=1.

Example 14.

Let Bλ,k,j:[0,1]→[0,1] be given by

Bλ,k,j⁢(x)=cj⁢tλ,k,m⁢(x),

where

tλ,k,m⁢(x)=(k+1)⁢xq(λ+1)⁢[xq+(1−x)q],

for all x∈[0,1] and cj≥0 such that ∑j=1mcj=1 and q≥1.

Example 15.

Let Bλ,k,j:[0,1]→[0,1] be given by

Bλ,k,j⁢(x)=cj⁢tλ,k,m⁢(x),

where

tλ,k,m⁢(x)=x2⁢(1+sin⁡(π⁢kλ)),

for all x∈[0,1] and cj≥0 such that ∑j=1mcj=1.

Remark 16.

Based on the examples presented above, we consider the functions f1,f2,f3:[0,1]→ℝ, given by

f1⁢(x) =0.5⁢x2+0.3⁢sin⁡(8⁢π⁢x)+0.2⁢x,
f2⁢(x) =0.3⁢sin⁡(4⁢π⁢x)⁢cos⁡(6⁢π⁢x)+0.5⁢x2+0.2,
f3⁢(x) =sin⁡(5⁢x)⁢e−x,

for x∈[0,1]. In Figure 1 we display the functions fi together with their corresponding operators ℳλm⁢fi, i=1,2,3 constructed using the functions Aλ,k and Bλ,k,j given in

  • •

    Example 13 for i=1,

  • •

    Example 14 for i=2,

  • •

    Example 15 for i=3, respectively.

Finally, Table 1 presents the maximum approximation errors for the operators ℳλm⁢fi, i=1,2,3, considering all the three cases given in Examples 13–15 and Table 2 presents the mean approximation errors. We consider the following parameters: m=4, r1=2,r2=1,r3=3,r4=2, p=8, λ=n+p−R with R=r1+r2+r3+r4 and n∈{50,100,200,500,1000}.

Refer to caption

f1⁢(x)

Refer to caption

f2⁢(x)

Refer to caption

f3⁢(x)

Figure 1. Functions fi,i=1,2,3, and the corresponding ℳλ4⁢fi operators.
fi\n Example 50 100 200 500 1000
f1 1 0.2549 0.1724 0.1014 0.0448 0.0231
2 0.3118 0.1845 0.1041 0.0445 0.0229
3 0.2742 0.1793 0.1059 0.0469 0.0243
f2 1 0.1580 0.1142 0.0726 0.0339 0.0179
2 0.1594 0.1197 0.0754 0.0349 0.0184
3 0.1699 0.1252 0.0788 0.0367 0.0194
f3 1 0.0472 0.0236 0.0118 0.0047 0.0024
2 0.0411 0.0207 0.0106 0.0043 0.0022
3 0.0631 0.0321 0.0162 0.0065 0.0033
Table 1. Maximum approximation errors for ℳλ4⁢fi, i=1,2,3.
fi\n Example 50 100 200 500 1000
f1 1 0.1195 0.0767 0.0438 0.0190 0.0098
2 0.1465 0.0912 0.0512 0.0219 0.0112
3 0.1367 0.0864 0.0488 0.0211 0.0108
f2 1 0.0724 0.0513 0.0314 0.0143 0.0075
2 0.0841 0.0588 0.0353 0.0159 0.0083
3 0.0798 0.0561 0.0340 0.0154 0.0080
f3 1 0.0210 0.0105 0.0052 0.0021 0.0010
2 0.0272 0.0139 0.0070 0.0028 0.0014
3 0.0312 0.0157 0.0079 0.0032 0.0016
Table 2. Mean approximation errors for ℳλ4⁢fi, i=1,2,3.

5. A nonlinear extension with data-dependent weights

In this section, we present a nonlinear extension of the operator ℳλm introduced in Section 2. By weakening the assumptions on the functions A and B defined in relation (6), we will consider that the weights appearing in the operator’s expression are themselves Bernstein-type operators. Thus, we obtain the operator ℳλB defined by

(16) ℳλB⁢f⁢(x)=∑k=0λbλ,k⁢(x)⁢[(1−ℬn⁢f⁢(x))⁢f⁢(kn)+ℬn⁢f⁢(x)⁢f⁢(k+rn)],

where ℬn⁢f is the classical Bernstein operator given by (1) and bλ,k are the Bernstein polynomials given by (2). Note that in this section we consider the particular case m=1 and then λ=n+p−R reduces to λ=n+p−r, where n,p,r∈ℕ with p≥0 and n>2⁢r. It is clear that the operator is no longer linear, given the dependence of the weights on the function f. Moreover, the positivity of the operator is preserved only under certain conditions. Nevertheless, for the introduced operator, important properties in approximation theory can still be studied, as can be observed in this final section.

Remark 17.

Let ℳλB be given by (16) and f∈C⁡[0,1+p]. Then

ℳλB⁢f⁢(0)=f⁡(0)

and

ℳλB⁢f⁢(1)=f⁡(n+pn).

In particular, if p=0, then ℳλB⁢f⁢(1)=f⁡(1). The proof of this result is very simple, according to the well-known properties of the Bernstein operator ℬn⁢f.

Proposition 18.

Let ej⁢(x)=xj, for every j∈{0,1,2} and x∈[0,1]. Then

ℳλB⁢(e0,x)=1

and

ℳλB⁢(e1,x)=n+pn⁢x,

for all x∈[0,1]. Moreover,

ℳλB⁢(e2,x)=1n2⁢[λ⁢x⁢(1−x)+λ2⁢x2+(2⁢r⁢λ⁢x+r2)⁢(x2+x⁡(1−x)n)],

for all x∈[0,1].

Proof.

In view of Theorem 6 for m=1 and Bλ,k,j⁢(x)=ℬn⁢f⁢(x), for all x∈[0,1], where ℬn⁢f is the Bernstein operator given by (1). ∎

It is not difficult to observe that for p=0 we obtain better results for the operator ℳλB (in particular, we have that ℳλB⁢(e1,x)=x, for all x∈[0,1]) that are similar to those obtained for the Bernstein operator ℬn⁢f (see [21]).

Acknowledgements.

The authors thank the referee(s) for carefully reading the manuscript and providing helpful suggestions.

References

  • [1]
  • [2] G. Başcanbaz-Tunca, A. Erençin and F. Taşdelen, Some properties of Bernstein type Cheney and Sharma operators, Gen. Math., 24 (2016) nos.1-2, pp. 17-25.
  • [3] S.N. Bernstein, Démonstration du théoréme de Weierstrass fondée sur le calcul des probabilités, Commun. Kharkov Math. Soc., 13 (1912/1913) no. 1, pp. 1-2.
  • [4] T. Bostanci and G. Başcanbaz-Tunca, A Stancu type extension of Cheney and Sharma operator, J. Numer. Anal. Approx. Theory, 47 (2018) no. 2, pp. 124–134. https://doi.org/10.33993/jnaat472-1133
  • [5] T. Cătinaş, Extension of Some Cheney-Sharma Type Operators to a Triangle With One Curved Side, Miskolc Math. Notes 21 (2020), pp. 101–111. https://dx.doi.org/10.18514/MMN.2020.2686
  • [6] T. Cătinaş, Cheney–Sharma Type Operators on a Triangle with Straight Sides, Symmetry 14 (2022), no. 11, 2446. https://doi.org/10.3390/sym14112446
  • [7] T. Cătinaş, A Stancu type extension of the Campiti-Metafune operator, J. Numer. Anal. Approx. Theory, 54 (2025) no. 2, pp. 229–236. DOI: 10.33993/jnaat542-1638
  • [8] T. Cătinaş and I. Buda, An extension of the Cheney-Sharma operator of the first kind, J. Numer. Anal. Approx. Theory, 52 (2023) no. 2, pp. 172–181. https://doi.org/10.33993/jnaat522-1373
  • [9] T. Cătinaş and D. Otrocol, Iterates of multivariate Cheney-Sharma operators, J. Comput. Anal. Appl. 15 (2013), No. 7, pp. 1240–1246.
  • [10] T. Cătinaş and D. Otrocol, Iterates of Cheney-Sharma type operators on a triangle with curved side, J. Comput. Anal. Appl. 28 (2020), No. 4, pp. 737–744.
  • [11] N. Çetin, A new generalization of complex Stancu operators, Math. Methods Appl. Sci. 42 (2019), pp 5582–5594.
  • [12] N. Çetin, A new complex generalized Bernstein-Schurer operator, Carpathian J. Math., 37 (2021) no. 1, pp. 81-89.
  • [13] N. Çetin and N.M. Mutlu, Complex generalized Stancu-Schurer operators, Math. Slovaca, 74 (2024) no.5, pp. 1215-1232.
  • [14] E.W. Cheney and A. Sharma, On a generalization of Bernstein polynomials, Riv. Mat. Univ. Parma, 2 (1964), pp. 77-84.
  • [15] S.G. Gal, Approximation by Complex Bernstein and Convolution Type Operators, Series on Concrete and Applicable Mathematics, vol. 8, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2009.
  • [16] E.S. Grigoriciuc, A Stancu type extension of the Cheney-Sharma Chlodovsky operators, J. Numer. Anal. Approx. Theory, 53 (2024) no. 1, pp. 103-117. DOI: 10.33993/jnaat531-1406
  • [17] E.S. Grigoriciuc and A. Malina, A Stancu–Schurer type extension of higher order of the Cheney–Sharma operators, Math. Slovaca, 76 (2026) no. 1, pp. 225-244. DOI: 10.1515/ms-2025-1155
  • [18] D. Miclăuş, The generalization of some results for Schurer and Schurer-Stancu operators, Rev. Anal. Numér. Théor. Approx., 40 (2011) no. 1, pp. 52-63.
  • [19] F. Schurer, Linear positive operators in approximation theory, Math. Inst. Technol. Univ. Delft Rep., 1962.
  • [20] D. Söylemez and F. Taşdelen, On Cheney-Sharma Chlodovsky operators, Bull. Math. Anal. Appl., 11 (2019) no. 1, pp. 36-43.
  • [21] D.D. Stancu Quadrature formulas constructed by using certain linear positive operators, in Numerical Integration. ISNM 57: International Series of Numerical Mathematics, vol. 57, G. Hämmerlin, Ed., Basel, Birkhäuser, 1981, pp. 241-251.
  • [22]