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 50th 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 m1. 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 fC[0,1] and n, Bernstein introduced in [3] the linear and positive operators

(1) nf(x)=k=0nbn,k(x)f(kn),

with the basis polynomials bn,k given by

(2) bn,k(x)=(nk)xk(1x)nk,

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) Snf(x)=k=0nrbnr,k(x)[(1x)f(kn)+xf(k+rn)],

where the polynomials bnr,k are given by (2), fC[0,1], x[0,1] and r is a non-negative integer satisfying n>2r, for every n.

Let m with m1. Also, let n,r1,,rm,R be such that n>2R, where R=r1++rm. For every function fC[0,1], the authors introduced in [17] the generalized Stancu operator of order m, given by

(4) Snmf(x)=k=0nRbnR,k(x)[(1j=1mxj)f(kn)+j=1mxjf(k+rjn)],

where the polynomials bnR,k are given by (2), for all xDm, where

Dm={x[0,1]: 1j=1mxj0}[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,512],D3=[0,α],

where α=13(17+3333+1733331). In general,

Dm=[0,αm],

where αm[0,1] is the root of the equation x+x2++xm1=0. Moreover, if m and x[0,1), then

limm(x+x2++xm)1=x1x1=2x11x

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 p0 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,pmf(x)=k=0n+pRbn+pR,k(x)[(1j=1mxj)f(kn)+j=1mxjf(k+rjn)],

for all xDm, where the polynomials bn+pR,k are given by (2) and fC[0,1+p].

An important step that can be undertaken is to replace the Bernstein polynomials bn+pR,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 m1. Also, let n,r1,,rm,R be such that n>2R, where R=r1++rm. For simplicity, let us denote λ=n+pR. Then, for every function fC[0,1+p], we denote by

(5) λmf(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)=1j=1mxj,

for all xDm, 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)=1tλ,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

λmf(x) =k=0λbλ,k(x)[(1tλ,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

λ1f(x)=k=0λbλ,k(x)[(1tλ,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 nr1 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,gC[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)]
=αλmf(x)+βλmg(x)
=(αλmf+βλmg)(x),

for all xDm. Moreover, since bλ,k, Aλ,k and Bλ,k,j are positive functions on [0,1], it follows that for any positive function fC[0,1+p] we obtain that

λmf(x)0,xDm,

and this completes the proof. ∎

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

Corollary 5.

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

λmf(x)λmg(x),

for all xDm.

Proof.

Let f,gC[0,1+p] be two functions such that f(x)g(x), for all x[0,1+p]. Then we consider hC[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 λmh(x)0, for all xDm according to Proposition 4. Moreover, based on the linearity of the operator λm we deduce that

0λmh(x)=λm(gf)(x)=λmg(x)λmf(x)

and then

λmf(x)λmg(x),

for all xDm 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 xDm. Then

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

and

(11) λm(e1,x)=λnx+1nk=0λbλ,k(x)j=1mrjBλ,k,j(x),

for all xDm. Moreover,

(12) λm(e2,x) =λ(λ1)n2x2+λn2x+2n2k=0λkbλ,k(x)j=1mrjBλ,k,j(x)
+1n2k=0λbλ,k(x)j=1mrj2Bλ,k,j(x),

for all xDm.

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 xDm. 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=1mrjnBλ,k,j(x)
=λnx+1nk=0λbλ,k(x)j=1mrjBλ,k,j(x),

for all xDm. 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+2n2k=0λkbλ,k(x)j=1mrjBλ,k,j(x)
+1n2k=0λbλ,k(x)j=1mrj2Bλ,k,j(x)
=λ2n2k=0λk2λ2bλ,k(x)+2n2k=0λkbλ,k(x)j=1mrjBλ,k,j(x)
+1n2k=0λbλ,k(x)j=1mrj2Bλ,k,j(x)
=λ2n2(x2+x(1x)λ)+2n2k=0λkbλ,k(x)j=1mrjBλ,k,j(x)
+1n2k=0λbλ,k(x)j=1mrj2Bλ,k,j(x)
=λ(λ1)n2x2+λn2x+2n2k=0λkbλ,k(x)j=1mrjBλ,k,j(x)
+1n2k=0λbλ,k(x)j=1mrj2Bλ,k,j(x),

for all xDm 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 xDm and any m with m1, the following relations hold:

  1. (1)

    limnλm(e1,x)=x;

  2. (2)

    λm(e1,x)[0,λ+rmaxn][0,1+p], where rmax=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 xDm. On the other hand, according to the assumption imposed on Bλ,k,j in Section 2, it is not difficult to observe that

0j=1mrjBλ,k,j(x)rmaxj=1mBλ,k,j(x)=rmax,

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

0λnxλm(e1,x)λx+rmaxnλ+rmaxn,

for all xDm. Moreover, it is clear that

rmax=max{r1,,rm}R=r1++rm,

and then

λ+rmaxn=n+pR+rmaxnn+pn=1+pn1+p,

for all n,p with p0 and n>2R. 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 fC[0,1+p]. Then

λmfC(Dm)fC[0,1+p].
Proof.

Let xDm and fC[0,1+p]. Using properties (6) and (7), we have that

|λmf(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)|]
fC[0,1+p]k=0λbλ,k(x)(Aλ,k(x)+j=1mBλ,k,j(x))
=fC[0,1+p]k=0λbλ,k(x)
=fC[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

λmf(x)f(λm(e1,x)),

for all xDm.

Proof.

In view of relation (5), we know that

λmf(x)=k=0λbλ,k(x)[Aλ,k(x)f(kn)+j=1mBλ,k,j(x)f(k+rjn)],

for all xDm. 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=1i=1mαk,i[0,1].

Then the operator λm can be written as

λmf(x)=k=0λbλ,k(x)[αk,0f(y0)+i=1mαk,if(yi)]=k=0λbλ,k(x)i=0mαk,if(yi),

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

λmf(x) =k=0λbλ,k(x)i=0mαk,if(yi)
k=0λbλ,k(x)f(i=0mαk,iyi)
f(k=0λbλ,k(x)i=0mαk,iyi)
=f(k=0λbλ,k(x)(αk,0y0+i=1mαk,iyi))
=f(k=0λbλ,k(x)(Aλ,k(x)kn+j=1mBλ,k,j(x)k+rjn))
=f(λm(e1,x)),

for all xDm. Hence,

λmf(x)f(λm(e1,x)),

for all xDm 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 pR, then λm(e1,x)x, for all xDm.

Proof.

In view of Theorem 6 we know that

λm(e1,x)=λnx+1nk=0λbλ,k(x)j=1mrjBλ,k,j(x),

for all xDm. Then

λm(e1,x)x =λnnx+1nk=0λbλ,k(x)j=1mrjBλ,k,j(x)
=pRnx+1nk=0λbλ,k(x)j=1mrjBλ,k,j(x)
pRnx,

since bλ,k and Bλ,k,j are positive functions on [0,1] and n,r1,,rm. Hence, if pR, then

λm(e1,x)x0,

for all xDm 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 pR and f is an increasing convex function on [0,1+p], then

(13) λmf(x)f(x),

for all xDm.

Proof.

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

(14) λmf(x)f(λm(e1,x)),

for all xDm. Moreover, since f is increasing on [0,1+p], it follows that for any two points u,v[0,1+p] with uv 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 xDm. Hence,

(15) f(x)=f(u)f(v)=f(λm(e1,x)),

for all xDm. In view of relations (14) and (15) we deduce that

λmf(x)f(x),

for all xDm 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)=cjtλ,k,m(x),

where

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

for all x[0,1] and cj0 such that j=1mcj=1.

Example 14.

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

Bλ,k,j(x)=cjtλ,k,m(x),

where

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

for all x[0,1] and cj0 such that j=1mcj=1 and q1.

Example 15.

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

Bλ,k,j(x)=cjtλ,k,m(x),

where

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

for all x[0,1] and cj0 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.5x2+0.3sin(8πx)+0.2x,
f2(x) =0.3sin(4πx)cos(6πx)+0.5x2+0.2,
f3(x) =sin(5x)ex,

for x[0,1]. In Figure 1 we display the functions fi together with their corresponding operators λmfi, 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 λmfi, i=1,2,3, considering all the three cases given in Examples 1315 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+pR 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 λ4fi 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 λ4fi, 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 λ4fi, 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) λBf(x)=k=0λbλ,k(x)[(1nf(x))f(kn)+nf(x)f(k+rn)],

where nf 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+pR reduces to λ=n+pr, where n,p,r with p0 and n>2r. 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 fC[0,1+p]. Then

λBf(0)=f(0)

and

λBf(1)=f(n+pn).

In particular, if p=0, then λBf(1)=f(1). The proof of this result is very simple, according to the well-known properties of the Bernstein operator nf.

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+pnx,

for all x[0,1]. Moreover,

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

for all x[0,1].

Proof.

In view of Theorem 6 for m=1 and Bλ,k,j(x)=nf(x), for all x[0,1], where nf 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 nf (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]