Next Article in Journal
Solving an Integral Equation via Tricomplex-Valued Controlled Metric Spaces
Next Article in Special Issue
On Perfectness of Systems of Weights Satisfying Pearson’s Equation with Nonstandard Parameters
Previous Article in Journal
A Variational Formulation for Fins with Nonzero Contact Thermal Resistance at the Base
Previous Article in Special Issue
Conditional Expanding of Functions by q-Lidstone Series
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quadrature Methods for Singular Integral Equations of Mellin Type Based on the Zeros of Classical Jacobi Polynomials

Faculty of Mathematics, Chemnitz University of Technology, 09107 Chemnitz, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Axioms 2023, 12(1), 55; https://doi.org/10.3390/axioms12010055
Submission received: 14 November 2022 / Revised: 27 December 2022 / Accepted: 28 December 2022 / Published: 3 January 2023
(This article belongs to the Special Issue Orthogonal Polynomials, Special Functions and Applications)

Abstract

:
In this paper we formulate necessary conditions for the stability of certain quadrature methods for Mellin type singular integral equations on an interval. These methods are based on the zeros of classical Jacobi polynomials, not only on the Chebyshev nodes. The method is considered as an element of a special C * -algebra such that the stability of this method can be reformulated as an invertibility problem of this element. At the end, the mentioned necessary conditions are invertibility properties of certain linear operators in Hilbert spaces. Moreover, for the proofs we need deep results on the zero distribution of the Jacobi polynomials.

1. Introduction

The present paper is part of the efforts done during the last three decades to establish necessary and sufficient conditions for the stability of numerical methods for singular integral equations by using so called C * -algebra techniques. The integral equations under consideration contain strong singular integral operators of Cauchy and Mellin type. In general, they are of the form
a ( x ) u ( x ) + b ( x ) π i 1 1 u ( y ) d y y x + c ( x ) 1 1 H 1 + x 1 + y d y 1 + y + c + ( x ) 1 1 H + 1 x 1 y d y 1 y + 1 1 K ( x , y ) u ( y ) d y = f ( x ) , 1 < x < 1 ,
where the functions a , b , c ± : [ 1 , 1 ] C , f : ( 1 , 1 ) C , H ± : R + C as well as K : ( 1 , 1 ) × ( 1 , 1 ) C are given and u : ( 1 , 1 ) C is looked for. As usual, R and C denote the sets of real and complex numbers, respectively. Moreover, by R + = t R : t > 0 we refer to the set of positive real numbers. The minimal conditions on the given functions are the piecewise continuity of the coefficient functions a , b , and c ± as well as the continuity of the kernel functions H ± ( t ) and K ( x , y ) . Moreover, the right-hand side f should belong to a Hilbert space L α , β 2 . Equation (1) is considered in this space L α , β 2 and written shortly as
A u : = a I + b S + c M H + c + M H + + + K u = f .
For real numbers α , β > 1 , the Hilbert space L α , β 2 : = L α , β 2 ( 1 , 1 ) is defined by the inner product
f , g α , β = 1 1 f ( x ) g ( x ) ¯ v α , β ( x ) d x ,
where v α , β ( x ) = ( 1 x ) α ( 1 + x ) β is a classical Jacobi weight. Hence the norm in L α , β 2 is given by f α , β = f , f α , β . We call a function a : [ 1 , 1 ] C piecewise continuous if it is continuous at ± 1 , the one-sided limits a ( x ± 0 ) exist for all x ( 1 , 1 ) , and at least one of them coincides with the function value a ( x ) . The set of these piecewise continuous functions is denoted by PC : = PC [ 1 , 1 ] .
For a continuous function H : ( 0 , ) C , i.e., H C ( R + ) , and a continuous function K : ( 1 , 1 ) × ( 1 , 1 ) C the Mellin-type operator M H and the integral operator K are defined by
( M H ± u ) ( x ) = 1 1 H 1 x 1 y u ( y ) d y 1 y , x ( 1 , 1 ) ,
and
( K u ) ( x ) = 1 1 K ( x , y ) u ( y ) d y , x ( 1 , 1 ) ,
respectively. Moreover, S denotes the Cauchy singular integral operator defined in the sense of a principal value integral as
S u ( x ) = 1 π i 1 1 u ( y ) d y y x , x ( 1 , 1 ) .
Furthermore, the operator of multiplication by a bounded function a : [ 1 , 1 ] C is defined by
a I : L α , β 2 L α , β 2 , f a f , a PC .
Thus, I : L α , β 2 L α , β 2 itself denotes the identity operator. If B : L α , β 2 L α , β 2 is another operator, then we use the abbreviation a B for the product of the multiplication operator a I and the operator B , i.e., a B : = a I B .
The numerical methods for the approximate solution of Equation (1), which are of interest here, are collocation and collocation-quadrature methods, which we will describe later on in more detail. Since we know that solutions of (1) usually contain singularities at the endpoints of the integration interval, in these methods we look for an approximate solution u n ( x ) to u ( x ) of the form
u n ( x ) = v ρ 0 , τ 0 ( x ) p n ( x ) = ( 1 x ) ρ 0 ( 1 + x ) τ 0 p n ( x ) ,
where ρ 0 , τ 0 > 1 are real numbers and p n ( x ) is a polynomial of degree less than n . In case of a collocation method we choose a sequence of collocation points
1 < x n n < x n 1 , n < < x 1 n < 1 , k = 1 , , n ,
where n N —the set of positive integers, and try to determine u n ( x ) with the help of the conditions
A u n ( x j n ) = f n ( x j n ) , j = 1 , , n ,
where f n v ρ 0 , τ 0 P n is an approximation to f , P n denotes the space of algebraic polynomials of degree less than n , and v ρ 0 , τ 0 P n is considered as a subspace of L α , β 2 (i.e., equipped with the norm of L α , β 2 ).
To realize a so called collocation-quadrature method, in a first step we approximate the integral operators M H ± and K with the help of a quadrature method of interpolation type
1 1 g ( x ) ω ( x ) d x k = 1 n λ k n g ( x k n ) , λ k n = 1 1 k n ( x ) ω ( x ) d x ,
where
k n ( x ) = j = 1 , j k n x x k n x j n x k n , k = 1 , , n ,
are the fundamental Lagrange interpolation polynomials with respect to the nodes x k n . Thus, the Mellin type operator is approximated by
M n , H ± u ( x ) = k = 1 n λ k n ω ( x k n ) H 1 x 1 x k n u ( x k n ) 1 x k n
and the Fredholm integral operator K by
K n u ( x ) = k = 1 n λ k n ω ( x k n ) K ( x , x k n ) u ( x k n ) .
In the second step, we again use the nodes x k n as collocation points and try to determine u n ( x ) by solving
a u n + b S u n + c M n , H u n + c + M n , H + + u n + K n u n ( x j n ) = f n ( x j n ) , j = 1 , , n .
Note that, in the collocation-quadrature method (9), the quadrature rule (8) is not applied to S u n ( x ) .
Both the collocation method (7) and the collocation-quadrature method (9) can be written as an operator equation
A n u n = f n ,
where A n : v ρ 0 , τ 0 P n v ρ 0 , τ 0 P n is a linear operator (cf. (44) and (45)). The definition of the stability of the method (10) or, in other words, of the stability of the sequence A n = A n n = 1 of the operators A n , includes the unique solvability of (10) for all sufficiently large n and the uniform boundedness of the inverse operators A n 1 : v ρ 0 , τ 0 P n v ρ 0 , τ 0 P n (see Definition 1).
Now, the application of C * -algebra techniques is based on the idea to consider the sequence A n as an element of a suitable C * -algebra and to translate stability into invertibility modulo zero sequences of this element (see Section 4). To find necessary and sufficient conditions for the stability of the sequence A n it is necessary to segue to certain C * -subalgebras and quotient algebras (cf. Proposition 2). In Table 1 we give an overview on the efforts done in the literature during the last 25 years to equations of type (1) or (2), where we ignore the Fredholm integral operator K .
The aim of the present paper is to extend the possible choices of ρ 0 and τ 0 in comparison to Table 1, where here we restrict to collocation-quadrature methods and the case a 1 and b 0 , i.e., the Cauchy singular integral operator is absent in Equation (1). To reach this aim we use the zeros of classical Jacobi polynomials associated with weights v γ , δ ( x ) as collocation and quadrature nodes, not only the zeros of Chebyshev polynomials. Unfortunately, for this general choice of collocation and quadrature nodes here we are only able to prove the necessity of the stability conditions. Their sufficiency will be the topic of a forthcoming paper.
There exists a series of papers (see, for example, refs. [17,18,19,20,21]) devoted to the application of the Nyström method to Fredholm integral equations of the form
( I + K 1 + K 2 ) u = f
with non-compact integral operators K 1 , for which the operators M H ± are respective examples. Thereby, Equation (11) is studied in spaces of continuous or weighted continuous functions. However, since the idea of proving stability and convergence of the Nyström method is essentially based on the concept of collectively compact operator sequences, which works only for compact operators K 1 + K 2 , in the mentioned papers there is assumed that the norm of the operator K 1 is less than 1 and that this is true uniformly also for the approximating operators K 1 n . Then, for I + K 1 and I + K 1 n one can use the Neumann series argument. In the present paper, we are not constrained to apply such a condition on the norm of an operator.

2. Preliminaries

2.1. Properties of Integral Operators with Mellin Kernels

Let us start with collecting some statements on integral operators of interest here and already proved in the literature.
Lemma 1 
([22], Proposition 3.13). Let β ( 1 , 1 ) and H C ( R + ) . Moreover, we assume that there are real numbers p , q with p < q such that 1 + β 2 ( p , q ) and such that
lim t + 0 t p k ( t ) = 0 a n d lim t t q k ( t ) = 0 .
Then, for all α ( 1 , 1 ) , the integral operator M H : L α , β 2 L α , β 2 is bounded.
Note that M H + = R M H R , where R : L α , β 2 L α , β 2 is given by ( R f ) ( x ) = f ( x ) .
Lemma 2 
([22], Lemma 3.8). Let α , β R . If the condition
1 1 1 1 1 x 1 y α 1 + x 1 + y β | K ( x , y ) | 2 d y d x <
is fulfilled, then K : L α , β 2 L α , β 2 is a compact operator.
The following corollary is an immediate consequence of the previous lemma.
Corollary 1. 
Let α , β ( 1 , 1 ) and η , ζ , ψ , χ R such that
η < 1 + α 2 , ζ < 1 + β 2 a n d ψ < 1 α 2 , χ < 1 β 2 .
If the function
( 1 , 1 ) × ( 1 , 1 ) C , ( x , y ) v η , ζ ( x ) K ( x , y ) v ψ , χ ( y )
is continuous and bounded, then the operator K : L α , β 2 L α , β 2 is compact.
For ρ R , we introduce the weighted L 2 -space L ρ 2 : = L ρ 2 ( R + ) defined by the norm
f ρ , 2 : = 0 | f ( t ) | 2 t ρ d t 1 2 .
Lemma 3. 
Let H L 2 s 1 2 and f L 0 , 2 s 1 2 for some s R . Then v 0 , s ( M H f ) is a bounded function on ( 1 , 1 ) .
Proof. 
For x ( 1 , 1 ) , we have
( 1 + x ) s ( M H f ) ( x ) 1 1 H 1 + x 1 + y ( 1 + x ) s ( 1 + y ) s + 1 ( 1 + y ) s | f ( y ) | d y 1 1 H 1 + x 1 + y 2 ( 1 + x ) 2 s ( 1 + y ) 2 s + 1 d y 1 2 1 1 ( 1 + y ) 2 s 1 | f ( y ) | 2 d y 1 2 = 1 + x 2 | H ( t ) | 2 t 2 s 1 d t 1 2 f 0 , 2 s 1 H 2 s 1 , 2 f 0 , 2 s 1 ,
from which the assertion follows. □
For n N 0 and H C n ( R + ) , we define the operators n M H by
( n M H f ) ( x ) = 1 1 H ( n ) 1 + x 1 + y f ( y ) ( 1 + y ) n + 1 d y .
Let R = R ( 1 , 1 ) and C = C ( 1 , 1 ) denote the sets of all functions f : ( 1 , 1 ) C being bounded and Riemann integrable as well as continuous on each closed subinterval of ( 1 , 1 ) , respectively. For S R , C and ψ , χ R with ψ , χ 0 , by S ˜ ψ , χ b we refer to the set of all functions f S , for which the function v ψ , χ f is bounded on ( 1 , 1 ) . If we introduce the norm
f ψ , χ , = sup v ψ , χ ( x ) | f ( x ) | : 1 < x < 1 ,
then S ˜ ψ , χ b , . ψ , χ , becomes a Banach space. Moreover, by S ˜ ψ , χ we denote the set of all functions f S ( 1 , 1 ) , for which the finite limits
lim x 1 0 ( 1 x ) ψ f ( x ) and lim x 1 + 0 ( 1 + x ) χ f ( x )
exist, and by S ψ , χ the subspace of S ˜ ψ , χ of those functions f S ( 1 , 1 ) , for which the limits in (12) are equal to zero if ψ > 0 or χ > 0 , respectively. The spaces S ˜ ψ , χ and S ψ , χ are closed subspaces of S ˜ ψ , χ b and, consequently, also Banach spaces. Finally, for ψ 0 , χ 0 > 0 , set
S ψ 0 , χ 0 0 = 0 ψ < ψ 0 , 0 χ < χ 0 S ψ , χ .
Note that
S ψ 0 , χ 0 0 = 0 ψ < ψ 0 , 0 χ < χ 0 S ψ , χ = 0 ψ < ψ 0 , 0 χ < χ 0 S ˜ ψ , χ = 0 ψ < ψ 0 , 0 χ < χ 0 S ˜ ψ , χ b = f S ( 1 , 1 ) : C > 0 , ε > 0 with | f ( x ) | < C v ε ψ 0 , ε χ 0 ( x ) x ( 1 , 1 ) .
Lemma 4. 
Let α , β < 1 and η , ζ , ψ , χ R such that
0 η , ζ a n d 0 ψ < 1 α 2 , 0 χ < 1 β 2 .
Moreover, assume the map
[ 1 , 1 ] × [ 1 , 1 ] C , ( x , y ) v η , ζ ( x ) K ( x , y ) v ψ , χ ( y )
to be continuous. Then K : L α , β 2 C ˜ η , ζ is a compact operator.
Proof. 
Let u L α , β 2 and ε > 0 . Then, there exists a δ = δ ( ε ) > 0 such that
( v η , ζ K u ) ( x ) ( v η , ζ K u ) ( x ) 1 1 v η , ζ ( x ) K ( x , y ) v η , ζ ( x ) K ( x , y ) | u ( y ) | d y ε 1 1 v 2 ψ α , 2 χ β ( y ) d y 1 2 u α , β const ε u α , β
for all x , x [ 1 , 1 ] with | x x | < δ . Thus, the set
v η , ζ K u : u L α , β 2 , u α , β 1
is equicontinuous. In the same way one can show
sup K u η , ζ , : u L α , β 2 , u α , β 1 const .
Applying the Arzela-Ascoli theorem delivers the assertion. □
For z C and a measurable function f : ( 0 , ) C , for which t z 1 f ( t ) is integrable on each compact subinterval of ( 0 , ) , the Mellin transform f ^ ( z ) is defined as
f ^ ( z ) = lim R R 1 R t z 1 f ( t ) d t ,
if this limit exists. Moreover, for ξ R and p < q , let Γ ξ = { z C : Re z = ξ } and
Γ p , q = z C : p < Re z < q , C 0 ( Γ ξ ) = f C ( Γ ξ ) : lim | η | f ( ξ + i η ) = 0 .
Lemma 5. 
Let p , q R with p < q and f L 2 ξ 1 2 C ( R + ) for every ξ ( p , q ) . Then
(a)
the Mellin transform f ^ belongs to the space C 0 ( Γ ξ ) for every ξ ( p , q ) and is holomorphic in the strip Γ p , q .
Moreover, if f ^ ( z ) satisfies
sup ( 1 + | z | ) 1 + τ f ^ ( z ) : z Γ p 0 , q 0 <
for all closed intervals [ p 0 , q 0 ] ( p , q ) and some τ = τ ( p 0 , q 0 ) > 0 , then
(b)
lim t + 0 t p + ε f ( t ) = 0 a n d lim t t q ε f ( t ) = 0
for every ε > 0 .
Proof. 
From f L 2 ξ 1 2 , p < ξ < q , we can conclude that f L ξ 1 1 , p < ξ < q (see [23] (Lemma 3.4)). This implies f ^ C 0 ( Γ ξ ) , p < ξ < q (cf. [23] (page 7)). By [22] (Lemma 2.14) we get that f ^ ( z ) is holomorphic in the strip Γ p , q , and (a) is proved. For assertion (b) we have only to refer to [23] (Corollary 3.3). □
For a function H C ( R + ) and a real number ξ , we formulate the following conditions:
( A 0 )
There exist real numbers p and q with p < q such that H L 2 p 1 2 ( R + ) L 2 q 1 2 ( R + ) ,   ξ ( p , q ) , and
sup ( 1 + | z | ) 1 + τ H ^ ( z ) : z Γ p 0 , q 0 <
for all intervals [ p 0 , q 0 ] ( p , q ) and some τ = τ ( p 0 , q 0 ) > 0 .
( A 1 )
There exist real numbers p and q with p < q such that H L 2 p 1 2 ( R + ) L 2 q 1 2 ( R + ) ,   ξ ( p , q ) , and
sup ( 1 + | z | ) s + τ H ^ ( s ) ( z ) : z Γ p 0 , q 0 < , s 0 , 1 ,
for all intervals [ p 0 , q 0 ] ( p , q ) and some τ = τ ( p 0 , q 0 ) > 0 .
We set
A : = I + c M H + c + M H + + ,
where c ± L ( 1 , 1 ) and H ± C ( R + ) . The following lemma is an application of [10] (Theorem 4.12) to the operator in (16).
Lemma 6. 
Let α , β ( 1 , 1 ) , c ± PC , and H ± C ( R + ) . If the function H = H ± satisfies condition ( A 0 ) for ξ = ξ ± , where ξ + = 1 + α 2 and ξ = 1 + β 2 , then the integral operator A : L α , β 2 L α , β 2 defined in (16) is Fredholm if and only if the closed curve
Γ A : = Γ A Γ A +
does not contain the point 0 , where
Γ A ± : = 1 + c ± ( ± 1 ) H ^ ± ( ξ ± i t ) : t R ¯ .
In this case, the Fredholm index of A is equal to the negative winding number of the curve Γ A , where the orientation of Γ A is due to the above given parametrizations of Γ A ± .
Lemma 7 
([10], Proposition 6.1). Let α , β ( 1 , 1 ) , a C \ 0 and k C ( R + ) . If the function H satisfies condition ( A 0 ) for ξ = 1 + β 2 , then the homogeneous equations ( a I + M H ) u = 0 in the space L α , β 2 or ( a I + M H ) v = 0 in the space L α , β 2 have only the trivial solution.

2.2. Marcinkiewicz Inequalities

For n N 0 and α , β > 1 , by P ^ n α , β ( x ) we denote the monic orthogonal polynomial of degree n with respect to the weight v α , β ( x ) . Furthermore, let x k n α , β , k = 1 , , n , be the zeros of P ^ n α , β ( x ) . It is well known that these zeros are real, simple, and contained in ( 1 , 1 ) . Therefore, we can write
x k n α , β = cos θ k n α , β , θ k n α , β ( 0 , π ) ,
and order them as follows
1 < x n n α , β < x n 1 , n α , β < < x 1 n α , β < 1 , i . e . , 0 < θ 1 n α , β < < θ n n α , β < π .
We set θ 0 n α , β = 0 and θ n + 1 , n α , β = π as well as x n + 1 , n α , β = 1 and x 0 n α , β = 1 .
The monic Jacobi polynomials P ^ n α , β ( x ) satisfy the three-term recurrence relation (cf. [24] (Chapter V, (2.7),(2.29))
P ^ n + 1 α , β ( x ) = ( x α n ) P ^ n α , β ( x ) β n 2 P ^ n 1 α , β ( x ) , n N 0 ,
where P ^ 1 α , β ( x ) 0 , P ^ 0 α , β ( x ) 1 , and
α n = α n α , β = β 2 α 2 ( 2 n + α + β ) ( 2 n + α + β + 2 ) , n N , β n = β n α , β = 4 n ( n + α ) ( n + β ) ( n + α + β ) ( 2 n + α + β 1 ) ( 2 n + α + β ) 2 ( 2 n + α + β + 1 ) , n N \ 1 ,
and
α 0 = α 0 α , β = α + β 2 + α + β , β 1 = β 1 α , β = 4 ( 1 + α ) ( 1 + β ) ( 2 + α + β ) 2 ( 3 + α + β ) .
By p n α , β ( x ) we refer to the normalized (with respect to the inner product . , . α , β ) polynomials with positive leading coefficient. If we set
β 0 = β 0 α , β : = 1 1 v α , β ( x ) d x = 2 α + β + 1 Γ ( α + 1 ) Γ ( β + 1 ) Γ ( α + β + 2 ) ,
then (see [24] (Chapter I, (4.10))
p n α , β ( x ) = γ n α , β P ^ n α , β ( x ) with γ n α , β = β 0 β 1 β 2 β n 1
and from (17) we get
β n + 1 p n + 1 α , β ( x ) = ( x α n ) p n α , β ( x ) β n p n 1 α , β ( x ) , n N 0 ,
where p 1 α , β ( x ) 0 and p 0 α , β ( x ) β 0 1 . In view of
1 1 ( 1 x 2 ) p n 1 α + 1 , β + 1 ( x ) x k v α , β ( x ) d x = 0 , k = 0 , 1 , , n 2 ,
we have
( 1 x 2 ) p n 1 α + 1 , β + 1 ( x ) = A n p n 1 α , β ( x ) + B n p n α , β ( x ) + C n p n + 1 α , β ( x )
with certain real numbers A n , B n , and C n . Together with (20) this yields
1 x n k α , β 2 p n 1 α + 1 , β + 1 ( x k n α , β ) = A n C n β n β n + 1 p n 1 α , β ( x k n α , β ) .
Let β n = β n α , β and δ n = β n α + 1 , β + 1 . With the help of the relations (18) and Γ ( z + 1 ) = z Γ ( z ) we get
δ 0 β 0 = 2 ( α + 1 ) ( β + 1 ) ( α + β + 3 ) ( α + β + 2 ) .
Furthermore,
δ k β k + 1 = k ( k + α + β + 2 ) ( k + 1 ) ( k + α + β + 1 ) , k = 1 , 2 ,
Using the orthogonality properties of p n α , β ( x ) and (19) we obtain, for n > 1 ,
A n = A n ( α , β ) = 1 1 ( 1 x 2 ) p n 1 α + 1 , β + 1 ( x ) p n 1 α , β ( x ) v α , β ( x ) d x = γ n 1 α , β 1 1 p n 1 α + 1 , β + 1 ( x ) x n 1 v α + 1 , β + 1 ( x ) d x = γ n 1 α , β γ n 1 α + 1 , β + 1 = δ 0 δ 1 δ n 1 β 0 β 1 β n 1 = δ 0 δ 1 β 0 β 1 = 2 4 + α + β ( 2 + α ) ( 2 + β ) ( 2 + α + β + 1 ) ( 4 + α + β 1 ) ( 4 + α + β + 1 ) : n = 2 δ 0 β 0 β 1 · δ 1 δ n 2 β 2 β n 1 = 2 + α + β n + α + β ( 2 + α + β ) ( n 1 ) · · 4 ( n 1 ) ( n + α ) ( n + β ) ( n + α + β + 1 ) ( 2 n + α + β 1 ) ( 2 n + α + β ) 2 ( 2 n + α + β + 1 ) : n > 2 = 2 2 n + α + β ( n + α ) ( n + β ) ( n + α + β ) ( n + α + β + 1 ) ( 2 n + α + β 1 ) ( 2 n + α + β + 1 ) ,
and, for n = 1 ,
A 1 = A 1 ( α , β ) = δ 0 β 0 = 2 ( α + 1 ) ( β + 1 ) ( α + β + 2 ) ( α + β + 3 ) ,
as well as
C n = 1 1 ( 1 x 2 ) p n 1 α + 1 , β + 1 ( x ) p n + 1 α , β ( x ) v α , β ( x ) d x = γ n 1 α + 1 , β + 1 1 1 x n + 1 p n + 1 α , β ( x ) v α , β ( x ) d x = γ n 1 α + 1 , β + 1 γ n + 1 α , β = β 0 β 1 β n + 1 δ 0 δ 1 δ n 1 = β n β n + 1 A n , n N .
Hence, relation (22) can be written in the form
p n 1 α , β ( x k n α , β ) = c n 1 x k n α , β 2 p n 1 α + 1 , β + 1 ( x k n α , β ) ,
where
c n = 1 A n + β n 2 A n = A n A n 2 + β n 2 = 2 n + α + β 2 ( n + α + β + 1 ) ( 2 n + α + β 1 ) ( n + α ) ( n + β ) ( n + α + β ) ( 2 n + α + β + 1 ) , n > 1 ,
and
c 1 = α + β + 2 2 ( α + β + 2 ) ( α + 1 ) ( β + 1 ) ( α + β + 3 ) .
In what follows, by C we will denote a positive constant, which can assume different values at different places, and we will write C C ( x , n , ) to indicate that C does not depend on the parameters x , n , … If A = A ( x , n , ) and B = B ( x , n , ) are two positive functions depending on certain variables x , n , , then we will write A x , n , B if there is a constant C C ( x , n , ) > 0 such that
C 1 B ( x , n , ) A ( x , n , ) C B ( x , n , ) x , n ,
Let λ n α , β ( x ) stand for the n-th Christoffel function with respect to the weight v α , β ( x ) , i.e.,
λ n α , β ( x ) = j = 0 n 1 p j α , β ( x ) 2 1 , x [ 1 , 1 ] .
It is well-known that
λ n α , β ( x k n α , β ) = 1 1 k n α , β ( x ) v α , β ( x ) d x ,
where k n α , β ( x ) are the nth fundamental Lagrange interpolation polynomials
k n α , β ( x ) = p n α , β ( x ) ( x x k n α , β ) ( p n α , β ) ( x k n α , β ) = j = 1 , j k n x x j n α , β x k n α , β x j n α , β
with respect to the nodes x k n α , β , k = 1 , , n . For x [ 1 , 1 ] and n N , let φ ( x ) : = 1 x 2 , x [ 1 , 1 ] ,   Δ n ( x ) = φ ( x ) n + 1 n 2 , and v n α , β ( x ) : = 1 x + 1 n 2 α 1 + x + 1 n 2 β . Then (cf. [25] (Theorem 5))
λ n α , β ( x ) Δ n ( x ) v n α , β ( x ) x , n 1 .
Since Jacobi weights are so-called doubling weights (see, for example, [26] (Section 3.2.1, Exercise 3.2.4), we also have (see [27] (Theorem 1))
x k n α , β x k + 1 , n α , β Δ n ( t ) t , k , n 1 , k = 0 , 1 , , n , t x k n α , β , x k + 1 , n α , β , n N ,
which can equivalently be written as (see [28] (Theorem 3.2) and cf. [26] (Exercise 3.2.25))
θ k + 1 , n α , β θ k n α , β k , n 1 n , k = 0 , , n , n N .
Note that, due to (28),
θ k n α , β = j = 1 k θ j n α , β θ j 1 , n α , β k , n k n , k = 1 , , n + 1 ,
and
π θ k n α , β = j = k n θ j + 1 , n α , β θ j , n α , β k , n n k + 1 n , k = 0 , , n .
Hence, for k = 1 , , n ,
1 x k n α , β = 2 sin θ k n α , β 2 k , n θ k , n α , β k , n k n
and
1 + x k n α , β = 2 cos θ k n α , β 2 = 2 sin π θ k n α , β 2 k , n π θ k , n α , β k , n n k + 1 n
as well as
1 n k ( n k + 1 ) n 2 k , n 1 x k n α , β 2 .
which implies, due to (28),
x k n α , β x k + 1 , n α , β t , k , n φ ( t ) n , k = 1 , , n 1 , t x k + 1 , n α , β , x k n α , β .
Lemma 8. 
For n N , k = 1 , , n , and α , β , γ , δ > 1 , we have
λ n α , β ( x k n γ , δ ) k , n v α + 1 2 , β + 1 2 ( x k n γ , δ ) n .
Proof. 
In view of (26) we have
λ n α , β ( x k n γ , δ ) k , n 1 x k n γ , δ 2 n + 1 n 2 1 x k n γ , δ + 1 n 2 α 1 + x k n γ , δ + 1 n 2 β ,
and it remains to take into account
1 x k n γ , δ < 1 x k n γ , δ + 1 n ( 31 ) , ( 32 ) C 1 x k n γ , δ
and
1 x k n γ , δ 2 < 1 x k n γ , δ 2 + 1 n ( 33 ) C 1 x k n γ , δ 2 .
Lemma 9 
([29], Theorem 2.6). Let 2 γ 0 + α > 1 , 2 δ 0 + β > 1 , and consider a system
1 = x n n < x n 1 , n < < x 1 n = 1
of nodes x k n = cos ϑ k n , ϑ k n [ 0 , π ] , satisfying ϑ k n ϑ k 1 , n k , n 1 n for k = 2 , , n and n N . Moreover, let m be a fixed positive integer. Then there exists a positive constant C C ( n , Q ) , such that
k = 1 n λ n α , β ( x k n ) v n γ 0 , δ 0 ( x k n ) Q ( x k n ) 2 C v γ 0 , δ 0 Q α , β 2
holds true for all Q P m n , n N , where P n denotes the set of all algebraic polynomials of degree less than n .
Corollary 2. 
Assume 2 γ 0 + α > 1 , 2 δ 0 + β > 1 , and consider a system
1 < x n n < x n 1 , n < < x 1 n < 1
of nodes x k n = cos ϑ k n , θ k n [ 0 , π ] , satisfying ϑ k n ϑ k 1 , n k , n 1 n for k = 1 , , n + 1 and n N , where ϑ n + 1 , n = π and ϑ 0 n = 0 . Moreover, let m be a fixed positive integer. Then there exists a positive constant C C ( n , Q ) , such that (36) is satisfied for all Q P m n , n N .
Proof. 
Obviously
v n + 2 α , β ( x ) x , n v n α , β ( x ) a n d Δ n + 2 ( x ) x , n Δ n ( x ) , x [ 1 , 1 ] , n N .
Using (26) we get
λ n + 2 α , β ( x ) x , n λ n α , β ( x ) , x [ 1 , 1 ] .
Now, we consider the following system of nodes
x ˜ n + 2 , n + 2 : = 1 , x ˜ 1 , n + 2 : = 1 , x ˜ k + 1 , n + 2 = x k n , k = 1 , , n .
If we apply Lemma 9 to this system, we immediately arrive at our assertion. □
In the particular case x k n = x k n γ , δ and γ 0 = δ 0 = 0 , from (35) and Corollary 2 we get
1 n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) Q ( x k n γ , δ ) 2 C Q α , β 2 , n N , α , β , γ , δ > 1 ,
for all Q P m n and with C C ( n , Q ) .
Lemma 10 
([29], Theorem 2.7). If α , β , γ , δ > 1 and 1 2 < α γ , β δ < 3 2 , then there exists a positive constant C C ( n , Q ) such that
C Q α , β 2 k = 1 n λ n α , β x k n γ , δ Q x k n γ , δ 2 , Q P n , n N .
Corollary 3. 
If α , β , γ , δ > 1 and
(a)
1 2 < α γ , β δ < 3 2 ,
then there exists a positive constant C C ( n , Q ) , such that
C 1 Q α , β 2 1 n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) Q ( x k n γ , δ ) 2 C Q α , β 2
for all Q P n , n N , where the second inequality holds true without condition ( a ) .
Proof. 
Relation (37) and Lemma 10 together with (35) deliver the assertion. □

2.3. The Algebra alg T ( PC )

By alg T ( PC ) we denote the smallest C * -subalgebra of the algebra L ( 2 ) of all linear and bounded operators on the Hilbert space 2 generated by the Toeplitz matrices
T ( g ) = g ^ j k j , k = 0
with piecewise continuous generating functions
g ( t ) : = Z g ^ t
defined on the unit circle T : = { t C : | t | = 1 } and continuous on T \ ± 1 .
Of course, alg T ( PC ) is a C * -subalgebra of the C * -algebra L T ( 2 ) L ( 2 ) generated by all Toeplitz matrices T ( f ) with piecewise continuous generating function f : T C . It is well known (see Chapter 16 in [30]) that there exists an isometrical isomorphism smb from the quotient algebra L T ( 2 ) / K ( 2 ) ( K ( 2 ) —the ideal in L ( 2 ) of compact operators) onto the algebra C ( M ) , . of all complex valued and continuous functions on the compact space M = T × [ 0 , 1 ] , where the topology on M is defined by the neighborhoods (cf. Theorem 16.1 in [30])
U ε , δ ( e i η 0 , 0 ) : = ( e i η , λ ) : η 0 δ < η < η 0 , 0 λ 1 ( e i η 0 , λ ) : 0 λ < ε , U ε , δ ( e i η 0 , 1 ) : = ( e i η , λ ) : η 0 < η < η 0 + δ , 0 λ 1 ( e i η 0 , λ ) : ε λ 1 , U δ 1 , δ 2 ( e i η 0 , λ 0 ) : = ( e i η 0 , λ ) : λ 0 δ 1 < λ < λ 0 + δ 2
with 0 < δ 1 < λ 0 < 1 δ 2 < 1 , 0 < δ < 2 π , and 0 < ε < 1 .
Proposition 1 
([30], Theorem 16.2, [31], Theorem 4.97). The mapping smb has the following properties:
(a)
R alg T ( PC ) is Fredholm if and only if smb R ( t , λ ) 0 for every ( t , λ ) T × [ 0 , 1 ] , where smb R : = smb ( R ) .
(b)
If R alg T ( PC ) is Fredholm, then the index is equal to the negative winding number of the closed curve
Γ R : = { smb R ( e i s , 0 ) : 0 < s < π } { smb R ( 1 , s ) : 0 s 1 } { smb R ( e i s , 0 ) : 0 < s < π } { smb R ( 1 , s ) : 0 s 1 } ,
where the orientation of Γ R is due to the above given parametrization.
Lemma 11 
([9], Lemma 2.10). Let H C ( R + ) . If H fulfils condition ( A 1 ) for ξ = 1 2 , then, for every s > 0 , the matrix
M : = H j + s k + s 1 k + s j , k = 0
defines an operator M L ( 2 ) , which belongs to the algebra alg T ( PC ) , and its symbol is given by
smb M ( t , λ ) = H ^ 1 2 + i 2 π log λ 1 λ : t = 1 , 0 : t T \ 1 .

3. The Collocation-Quadrature Method

We consider the integral Equation (cf. (16))
A u = ( I + c M H + c + M H + + + K ) u = f ,
where c ± PC , f L α , β 2 , H ± C ( R + ) , and the kernel function K ( x , y ) of the integral operator K (cf. (4)) is supposed to be continuous on ( 1 , 1 ) × ( 1 , 1 ) . In order to get approximate solutions, we use a polynomial collocation-quadrature method. To introduce that method, we need some further notations. Let n N and γ , δ , ρ , τ > 1 be real numbers. For u : ( 1 , 1 ) C , the Lagrange interpolation operator L n γ , δ is defined by
L n γ , δ u = j = 1 n u ( x j n γ , δ ) j n γ , δ ,
j n γ , δ ( x ) = p n γ , δ ( x ) ( x x j n γ , δ ) p n γ , δ ( x j n γ , δ ) = k = 1 , k j n x x k n γ , δ x j n γ , δ x k n γ , δ .
To the integral operators M H ± and K , we associate the quadrature operators
M n , H ± : C ( 1 , 1 ) C ( 1 , 1 ) , u k = 1 n λ k n γ , δ H 1 · 1 x k n γ , δ v γ , δ u ( x k n γ , δ ) 1 x k n γ , δ
and
K n : C ( 1 , 1 ) C ( 1 , 1 ) , u k = 1 n λ k n γ , δ K ( · , x k n γ , δ ) v γ , δ u ( x k n γ , δ ) ,
respectively. For certain ρ , τ > 1 , the collocation-quadrature method seeks for approximations u n L α , β 2 of the form
u n = ϑ p n : = v ρ α 2 , τ β 2 p n , p n P n
to the exact solution of (40) by solving
I + c M n , H + c + M n , H + + + K n u n ( x k n γ , δ ) = f n ( x k n γ , δ ) , k = 1 , , n ,
where P n stands for the set of all algebraic polynomials of degree less than n and the functions f n : ( 1 , 1 ) C are continuous and satisfy ϑ 1 f n P n as well as
lim n f f n α , β = 0 .
We set
p ˜ n : = ϑ p n ρ , τ , n = 0 , 1 , 2 ,
Note that p ˜ n n = 0 forms a complete orthonormal system in L α , β 2 . Using the weighted fundamental Lagrange interpolation polynomials
˜ k n γ , δ ( x ) : = ϑ ( x ) k n γ , δ ( x ) ϑ ( x k n γ , δ ) , k = 1 , 2 , , n ,
we can write u n as
u n = j = 0 n 1 α j n p ˜ j = k = 1 n ξ k n ˜ k n γ , δ .
If we introduce the Fourier projections
L n : L α , β 2 L α , β 2 , u j = 0 n 1 u , p ˜ j α , β p ˜ j
and the weighted Lagrange interpolation operator
L ˜ n γ , δ : = ϑ L n γ , δ ϑ 1 I ,
then the collocation system (42) can be written as
A n u n = f n
with
A n = L ˜ n γ , δ I + c M n , H + c + M n , H + + + K n L n .
Note also that, with the introduced notations, the assertion of Corollary 3 remains true for Q im L n , i.e.,
C 1 Q α , β 2 1 n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) Q ( x k n γ , δ ) 2 C Q α , β 2 , Q im L n , n N ,
if 1 2 < ρ γ , τ δ < 3 2 . Indeed, for Q = ϑ P im L n , we have
C 1 Q α , β 2 = C 1 P ρ , τ 2 π n k = 1 n v ρ + 1 2 , τ + 1 2 ( x k n γ , δ ) P ( x k n γ , δ ) 2 C P ρ , τ 2 = C Q α , β 2
and
v ρ + 1 2 , τ + 1 2 ( x k n γ , δ ) P ( x k n γ , δ ) 2 = v α + 1 2 , β + 1 2 ( x k n γ , δ ) Q ( x k n γ , δ ) 2 .
It is well known that, in the investigation of numerical methods for operator equations, the stability of the respective operator sequences plays an essential role.
Definition 1. 
We call the sequence A n in (44) stable (in L α , β 2 ) if, for all sufficiently large n, the operators A n : im L n im L n are invertible and if the norms A n 1 L n L ( L α , β 2 ) are uniformly bounded.
If the method is stable and if A n L n converges strongly to A L ( L α , β 2 ) , then the operator A is injective. If additionally the image of A equals L α , β 2 , then (43) implies the L α , β 2 -convergence of the solution u n of (44) to the (unique) solution u L α , β 2 of (40). This can be seen from the estimate
L n u u n α , β = A n 1 L n ( A n L n u A n u n ) α , β A n 1 L n L ( L α , β 2 ) A n L n u A u α , β + f f n α , β .

4. C * -Algebra Framework

In order to investigate the stability of the collocation-quadrature method, we use specific C * -algebra techniques. With the help of those tools, we are able to transform our stability problem into an invertibility problem in an appropriate C * -algebra. The sequence ( A n ) is considered as an element of such a C * -algebra. To define that algebra, we need some operators and spaces.
Let n N and
X 1 = L α , β 2 , X 2 = X 3 = 2 , L n ( 1 ) = L n , L n ( 2 ) = L n ( 3 ) = P n ,
where
P n : 2 2 , ( ξ j ) j = 0 ( ξ 0 , , ξ n 1 , 0 , ) .
Let the operators F n L ( im P n ) be given by
F n ( ξ 0 , ξ 1 , , ξ n 1 , 0 , ) = ( F n ) 1 ( ξ 0 , ξ 1 , , ξ n 1 , 0 , ) = ( ξ n 1 , ξ n 2 , , ξ 0 , 0 , ) .
Moreover, for t T : = { 1 , 2 , 3 } , we define E n ( t ) : im L n im L n ( t ) by
E n ( 1 ) : = L n , E n ( 2 ) : = V n , E n ( 3 ) : = F n V n ,
where
V n u : = n 1 2 v α 2 + 1 4 , β 2 + 1 4 ( x 1 n γ , δ ) u ( x 1 n γ , δ ) , , n 1 2 v α 2 + 1 4 , β 2 + 1 4 ( x n n γ , δ ) u ( x n n γ , δ ) , 0 , .
All operators E n ( t ) are invertible with
( E n ( 1 ) ) 1 = L n | im L n ,
( E n ( 2 ) ) 1 = n k = 1 n ξ k 1 v α 2 1 4 , β 2 1 4 ( x k n γ , δ ) ˜ k n γ , δ ,
( E n ( 3 ) ) 1 = n k = 1 n ξ n k v α 2 1 4 , β 2 1 4 ( x k n γ , δ ) ˜ k n γ , δ .
Moreover,
E n ( t ) u 2 2 = 1 n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) u ( x k n γ , δ ) 2 , u im L n , t = 2 , 3 .
Whence, in view of (46), the following lemma is proved.
Lemma 12. 
The operators E n ( j ) : im L n im L n ( t ) , n N , t = 1 , 2 , 3 , are uniformly bounded together with their inverses if 1 2 < ρ γ , τ δ < 3 2 .
Lemma 13. 
For r , t T with r t , the operators E n r , t : = E n ( r ) E n ( t ) 1 L n ( t ) converge together with their adjoints weakly to the zero operator if 1 2 < ρ γ , τ δ < 3 2 .
Proof. 
Since the operators E n r , t are uniformly bounded, it suffices to verify the convergence on a dense subset. At first, we consider the operator E n 2 , 1 = V n L n . Let k , m N 0 be arbitrary but fixed and n > max m , k as well as e k = ( δ j , k ) j = 0 2 , k N 0 . Using (31), we get
V n L n p ˜ m , e k 2 = n 1 2 v α 2 + 1 4 , β 2 + 1 4 ( x k + 1 , n γ , δ ) p ˜ m ( x k + 1 , n γ , δ ) C n ρ 1 .
Taking into account ρ > 1 , we conclude the weak convergence of E n 2 , 1 to the zero operator. Similarly, using (32) instead of (31), we get
F n V n L n p ˜ m , e k 2 = n 1 2 v α 2 + 1 4 , β 2 + 1 4 ( x n k , n γ , δ ) p ˜ m ( x n k , n γ , δ ) C n τ 1 .
Hence, E n 3 , 1 converges weakly to the zero operator. Fix k N 0 and p ˜ ( x ) = ϑ ( x ) ( 1 x 2 ) p ( x ) with p ( x ) being a polynomial. Then, for all sufficiently large n N ,
V n 1 P n e k , p ˜ α , β = n 1 2 v α 2 1 4 , β 2 1 4 ( x k + 1 , n γ , δ ) ˜ k + 1 , n γ , δ , ϑ v 1 , 1 p α , β = n 1 2 v ρ 2 1 4 , τ 2 1 4 ( x k + 1 , n γ , δ ) k + 1 , n γ , δ , v ρ γ , τ δ v 1 , 1 p γ , δ = n 1 2 v ρ 2 1 4 , τ 2 1 4 ( x k + 1 , n γ , δ ) k + 1 , n γ , δ , S n γ , δ v ρ γ , τ δ v 1 , 1 p γ , δ = n 1 2 v ρ 2 1 4 , τ 2 1 4 ( x k + 1 , n γ , δ ) λ k + 1 , n γ , δ S n γ , δ v ρ γ , τ δ v 1 , 1 p ( x k + 1 , n γ , δ )
(for the definition of S n γ , δ , see (51) below). Choose ψ , χ R such that
γ 2 1 4 < ψ < min 1 + γ , γ 2 + 3 4 , γ ρ 2 + 1 2
and
δ 2 1 4 < χ < min 1 + δ , δ 2 + 3 4 , δ τ 2 + 1 2 ,
which is possible since γ 2 1 4 < γ ρ 2 + 1 2 is equivalent to ρ γ < 3 2 . In virtue of [32] (2.2) there is a constant C C ( n , k , f ) such that
S n γ , δ f ( x k n γ , δ ) v ψ , χ ( x k n γ , δ ) C sup v ψ , χ ( x ) | f ( x ) | : 1 < x < 1 ln ( n + 1 ) , n N .
Note that, due to (35),
λ k n γ , δ v γ + 1 2 , δ + 1 2 ( x k n γ , δ ) = λ n γ , δ ( x k n γ , δ ) v γ + 1 2 , δ + 1 2 ( x k n γ , δ ) n , k 1 n , k = 1 , , n ,
Consequently,
V n 1 P n e k , p ˜ α , β C n 1 2 v γ ρ 2 + 1 4 ψ , δ τ 2 + 1 4 χ ( x k + 1 , n γ , δ ) v ψ , χ ( x k + 1 , n γ , δ ) S n γ , δ v ρ γ , τ δ v 1 , 1 p ( x k + 1 , n γ , δ ) ( 31 ) C n ρ 2 γ 1 + 2 ψ ln ( n + 1 ) 0
because of 2 ψ < 2 γ ρ + 1 . Thus, E n 1 , 2 converges weakly to zero. Analogously, we get the same for E n 1 , 3 .
Let s , t { 2 , 3 } with s t . The weak convergence of E n s , t follows by the relations E n 2 , 3 = V n ( F n V n ) 1 P n = F n = F n V n 1 V n P n = E n 3 , 2 and F n e k , e m = e n 1 k , e m = 0 if n > m + k + 1 . Finally, note that with E n r , t also E n r , t * converges weakly to the zero operator, which follows immediately from the definition of the weak convergence. The lemma is proved. □
For all what follows we assume that ρ γ , τ δ 1 2 , 3 2 . By F we denote the set of all sequences ( A n ) : = ( A n ) n = 1 of linear operators A n : im L n im L n for which the strong limits
W t ( A n ) : = lim n E n ( t ) A n E n ( t ) 1 L n ( t ) and W t ( A n ) * = lim n E n ( t ) A n E n ( t ) 1 L n ( t ) *
exist for all t T . If we provide F with the algebraic operations
λ 1 ( A n ) + λ 2 ( B n ) : = ( λ 1 A n + λ 2 B n ) , λ 1 , λ 2 C ,
( A n ) ( B n ) : = ( A n B n ) , ( A n ) * = ( A n * ) ,
and the supremum norm
( A n ) F : = sup n 1 A n L n L ( L α , β 2 ) ,
one can easily check, that F becomes a C * -algebra with the identity element ( L n ) .
Corollary 4. 
Let T t L X t , t T be compact operators, i.e., T t K X t . Then the sequences A n ( t ) with A n ( t ) = E n ( t ) 1 L n ( t ) T t E n ( t ) belong to F , where
E n ( r ) A n ( t ) E n ( r ) 1 L n ( r ) 0 a n d E n ( r ) A n ( t ) E n ( r ) 1 L n ( r ) * 0 , n ,
strongly in X r for r T , r t .
Proof. 
This is due to
E n ( r ) A n ( t ) E n ( r ) 1 L n ( r ) = E n r , t T t E n t , r a n d E n ( r ) A n ( t ) E n ( r ) 1 L n ( r ) * = E n t , r * T t * E n r , t * ,
the compactness of T t : X t X t and the weak convergence of E n r , t as well as E n r , t * to the zero operator if r t (see Corollary 4). □
Corollary 5. 
The mappings W t : F L X t ,   t T , are unital * -homomorphisms with norm 1.
Proof. 
Of course W t : F L X ( t ) , t T are*- homomorphisms. The relation W t L X ( t ) = 1 follows from the fact that * -homomorphism are bounded by 1 and that W t ( L n ) L X ( t ) = 1 .
The convergence
E n ( t ) A n ( t ) E n ( t ) 1 L n ( t ) = L n ( t ) T t L n ( t ) T t i n X t ( even in L ( X t ) )
together with Corollary 4 deliver that
J : = t T E n ( t ) 1 L n ( t ) T t E n ( t ) + C n : C n N , T t K X t , t T
is a subset of F , where N F is the two-sided closed ideal of F of all sequences ( A n ) F with A n L n L ( L α , β 2 ) 0 .
Proposition 2 
([33,34], Theorem 10.33). The set J forms a two-sided closed ideal in the C * -algebra F . Moreover, a sequence ( A n ) F is stable if and only if the operators W t ( A n ) : X t X t , t T , and the coset ( A n ) + J F / J are invertible.

5. The Limit Operators of the Collocation-Quadrature Methods

In this section, under certain conditions, we prove that the sequence A n of the collocation quadrature method, defined in (45), belongs to the algebra F from the previous section. We do this by determine the limit operators W t A n , t = 1 , 2 , 3 , and proving that also the sequences of the respective adjoint operators converge strongly.
The following Lemma is due to [35] (Satz III.2.1). Recall the definition of R ψ , χ 0 in (13) and the equalities (14).
Lemma 14. 
Let γ , δ > 1 .
(a)
If f R 1 + γ , 1 + δ 0 , then lim n k = 1 n λ k n γ , δ f ( x k n γ , δ ) = 1 1 f ( x ) v γ , δ ( x ) d x .
(b)
If f R 1 + γ 2 , 1 + δ 2 0 , then lim n L n γ , δ f f γ , δ = 0 .
Lemma 15. 
Let ψ , χ R , α , β , γ , δ > 1 , and f R 1 + α 2 , 1 + β 2 0 , i.e.,
| f ( x ) | C ( 1 x ) ε 1 + α 2 ( 1 + x ) ε 1 + β 2 , x ( 1 , 1 )
for some ε > 0 (cf. (14)). Then we have
v ψ , χ L n γ , δ v ψ , χ f f i n L α , β 2 ,
for n , if
1 2 < α + 2 ψ γ < 3 2 a n d α + 2 ψ > 1 ,
as well as
1 2 < β + 2 χ δ < 3 2 a n d β + 2 χ > 1 .
Proof. 
Let δ 1 > 0 be arbitrary chosen and n N . Since α + 2 ψ , β + 2 χ > 1 and v ψ , χ ( x ) f ( x ) 2 v α + 2 ψ , β + 2 χ ( x ) C ( 1 x 2 ) 2 ε 1 , we can choose a polynomial p ( x ) , such that f v ψ , χ p α , β 2 = v ψ , χ f p α + 2 ψ , β + 2 χ 2 < δ 1 . For n > deg p , we have
v ψ , χ L n γ , δ v ψ , χ f f α , β 2 2 v ψ , χ L n γ , δ ( v ψ , χ f p ) α , β 2 + v ψ , χ p f α , β 2 .
By using Corollary 3 for α + 2 ψ and β + 2 χ instead of α and β , relation (35) for γ , δ instead of α , β , and Lemma 14, (a) we get
v ψ , χ L n γ , δ ( v ψ , χ f p ) α , β 2 = L n γ , δ v ψ , χ f p α + 2 ψ , β + 2 χ 2 C n k = 1 n v 2 ψ + α + 1 2 , 2 χ + β + 1 2 ( x k n γ , δ ) ( v ψ , χ f p ) ( x k n γ , δ ) 2 C k = 1 n λ k n γ , δ v 2 ψ + α γ , 2 χ + β δ ( x k n γ , δ ) ( v ψ , χ f p ) ( x k n γ , δ ) 2 C 1 1 ( v ψ , χ f p ) ( x ) 2 v 2 ψ + α , 2 χ + β ( x ) d x = C f v ψ , χ p α , β 2
for n and some constant C C ( n , f , p ) , where we took into account
v 1 + γ , 1 + δ ( x ) v 2 ψ + α γ , 2 χ + β δ ( x ) v ψ , χ ( x ) f ( x ) 2 C ( 1 x 2 ) 2 ε
and
v 1 + γ , 1 + δ ( x ) v 2 ψ + α γ , 2 χ + β δ ( x ) | p ( x ) | 2 C v 1 + 2 ψ + α , 1 + 2 χ + β ( x ) C ( 1 x 2 ) ε 1
with ε 1 = min 1 + 2 ψ + α , 1 + 2 χ + β > 0 , which shows the applicability of Lemma 14, (a). Thus,
lim sup n v ψ , χ L n γ , δ v ψ , χ f f α , β 2 2 ( C + 1 ) δ 1 .
Since this is true for all δ 1 > 0 , we get the assertion. □
The following corollary is an immediate consequence of the previous lemma and concerned with the case ψ = ρ α 2 , χ = τ β 2 .
Corollary 6. 
Let α , β , γ , δ , ρ , τ > 1 . If f R 1 + α 2 , 1 + β 2 0 and ρ γ , τ δ 1 2 , 3 2 , then L ˜ n γ , δ f f in L α , β 2 .
Lemma 16. 
Let 0 ψ < 1 + α 2 , 0 χ < 1 + β 2 , and α γ , β δ 1 2 , 3 2 . Moreover, let
[ 1 , 1 ] × [ 1 , 1 ] C , ( x , y ) f ( x , y ) v ψ , χ ( y )
be a continuous function. Then
lim n sup L n γ , δ f ( x , · ) f ( x , · ) α , β : 1 x 1 = 0 .
Proof. 
Fix ψ 0 , χ 0 such that ψ < ψ 0 < 1 + α 2 and χ < χ 0 < 1 + β 2 . By assumption f x C ψ 0 , χ 0 for all x [ 1 , 1 ] , where f x ( y ) : = f ( x , y ) . Moreover,
lim x x 0 f x f x 0 ψ 0 , χ 0 , = 0 f o r a l l x 0 [ 1 , 1 ] .
Suppose the assertion of the lemma is not true. Then, there are an ε > 0 and a sequence n 1 < n 2 < of natural numbers, satisfying
sup L n k γ , δ f x f x α , β : 1 x 1 2 ε for   all k N .
Hence, for every k N , there is an x k [ 1 , 1 ] such that L n k γ , δ f x k f x k α , β ε , and we can assume that x k x * [ 1 , 1 ] for k . Due to our assumptions we have
M 0 : = 1 1 ( 1 x ) α 2 ψ 0 ( 1 + x ) β 2 χ 0 d x <
and, by Lemma 15 (choose ψ = χ = 0 and use α γ , β δ 1 2 , 3 2 ) and the Banach-Steinhaus theorem,
M 1 : = sup L n γ , δ C ψ 0 , χ 0 L α , β 2 : n N < .
Moreover, there is an k 0 N such that
L n k γ , δ f x * f x * α , β < ε 3 and f x k f x * ψ 0 , χ 0 , < ε 3 max M 0 , M 1 for   all k k 0 .
For k k 0 , we get the contradiction
ε L n k γ , δ f x k f x k α , β L n k γ , δ ( f x k f x * ) α , β + L n k γ , δ f x * f x * α , β + f x * f x k α , β M 1 f x k f x * ψ 0 , χ 0 , + ε 3 + M 0 f x * f x k ψ 0 , χ 0 , < ε ,
and the lemma is proved. □
Lemma 17. 
Let α , β ( 1 , 1 ) and η , ζ , ψ , χ R such that
0 η < 1 + α 2 , 0 ζ < 1 + β 2 a n d 0 ψ < 1 α 2 , 0 χ < 1 β 2 .
Moreover, let the map
[ 1 , 1 ] × [ 1 , 1 ] C , ( x , y ) v η , ζ ( x ) K ( x , y ) v ψ , χ ( y )
be continuous. If ρ γ , τ δ 1 2 , 3 2 holds true, then
lim n L ˜ n γ , δ K L n L n K L n L L α , β 2 = 0 .
That means L ˜ n γ , δ K L n J with
W 1 L ˜ n γ , δ K L n = K a n d W t L ˜ n γ , δ K L n = Θ , t { 2 , 3 } .
Proof. 
By Corollary 6 we have L ˜ n γ , δ f f in L α , β 2 for all f C ˜ η , ζ . Together with the compactness of the operator K : L α , β 2 C ˜ η , ζ (see Lemma 4) and the strong convergence L n * = L n I in L α , β 2 , this leads to the limit relation (48). Hence, L ˜ n γ , δ K L n = L n K L n + C n with C n N and consequently, by definition, L ˜ n γ , δ K L n J . Corollary 4 yields (49). □
Define
Ω = ( ω 1 , ω 2 ) R 2 : ω 1 , ω 2 > 1 , 1 2 < ω 1 ω 2 < 1 2 , ω 1 < 2 ω 2 + 1 ( ω 1 , ω 2 ) R 2 : ω 1 , ω 2 > 1 , 1 2 < ω 1 ω 2 < 3 2 .
Furthermore, for γ , δ > 1 , by S n γ , δ we refer to the Fourier operator given by
S n γ , δ f = j = 0 n 1 f , p j γ , δ γ , δ p j γ , δ .
From [36] (Theorem 1) we infer the following lemma.
Lemma 18. 
Let α , β , γ , δ , > 1 . Then, there is a constant C C ( n , f ) such that, for all n N and f L α , β 2 , the inequality
S n γ , δ f α , β C f α , β
holds true if and only if
| α γ | < min 1 2 , 1 + γ a n d | β δ | < min 1 2 , 1 + δ ,
which is equivalent to
1 2 < α γ , β δ < 1 2 , α < 2 γ + 1 , a n d β < 2 δ + 1 .
Corollary 7. 
Let a PC and α , β , γ , δ , γ 0 , δ 0 > 1 . If
1 2 < α γ , β δ < 3 2
and
1 2 < α γ 0 , β δ 0 < 1 2 , α < 2 γ 0 + 1 , β < 2 δ 0 + 1
are satisfied, then L n γ , δ a S n γ 0 , δ 0 a I strongly in L α , β 2 .
Proof. 
The set P of all algebraic polynomials is dense in L α , β 2 . Moreover, for every p P , we have L n γ , δ a S n γ 0 , δ 0 p = L n γ , δ a p for all sufficiently large n and, by Lemma 15, L n γ , δ a p a p in L α , β 2 . Additionally, in view of Corollary 3 and Lemma 18, for f L α , β 2 we can
L n γ , δ a S n γ 0 , δ 0 u α , β 2 C n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) a S n γ 0 , δ 0 f ( x k n γ , δ ) 2 a 2 C n k = 1 n v α + 1 2 , β + 1 2 ( x k n γ , δ ) S n γ 0 , δ 0 f ( x k n γ , δ ) 2 C a 2 S n γ 0 , δ 0 f α , β 2 C a 2 u α , β 2 .
Hence, the operators L n γ , δ a S n γ 0 , δ 0 : L α , β 2 L α , β 2 , n N , are uniformly bounded. Now, the Banach-Steinhaus theorem gives the assertion. □
In the following Corollary we apply the previous Corollary to operators of the form
ϑ v γ 0 ρ , δ 0 τ L n γ , δ a S n γ 0 , δ 0 ϑ 1 v ρ γ 0 , τ δ 0 I = v γ 0 ρ + α 2 , δ 0 τ + β 2 L n γ , δ a S n γ 0 , δ 0 v ρ + α 2 γ 0 , τ + β 2 δ 0 I .
Corollary 8. 
Let a PC and α , β , γ , δ , γ 0 , δ 0 , ρ , τ > 1 . Then,
ϑ v γ 0 ρ , δ 0 τ L n γ , δ a S n γ 0 , δ 0 ϑ 1 v ρ γ 0 , τ δ 0 I a I
strongly in L α , β 2 if
1 2 < 2 γ 0 ρ γ , 2 δ 0 2 τ δ < 3 2 ,
and
1 2 < ρ γ 0 , τ δ 0 < 1 2 , ρ < 2 γ 0 + 1 , τ < 2 δ 0 + 1 .
Proof. 
The claimed strong convergence is equivalent to the strong convergence of the operators L n γ , δ a S n γ 0 , δ 0 to a I in L 2 γ 0 ρ , 2 δ 0 τ 2 . For α = 2 γ 0 ρ and β = 2 δ 0 τ , the conditions α , β > 1 , (53) and (54) are equivalent to (55) and (56). Hence, Corollary 7 is applicable. □
Lemma 19. 
Let α , β ( 1 , 1 ) and η , ζ , ψ , χ R such that
0 η < 1 + α 2 , 0 ζ < 1 + β 2 a n d 0 ψ < 1 α 2 , 0 χ < 1 β 2 ,
and such that the function
[ 1 , 1 ] × [ 1 , 1 ] C , ( x , y ) v η , ζ ( x ) K ( x , y ) v ψ , χ ( y )
is continuous. Let ( ρ , γ ) , ( τ , δ ) Ω . Then lim n L ˜ n γ , δ K n L n L n K L n L L α , β 2 = 0 . That means (cf. Lemma 17) L ˜ n γ , δ K n L n J with
W 1 L ˜ n γ , δ K n L n = K a n d W t L ˜ n γ , δ K n L n = Θ , t { 2 , 3 } .
Proof. 
First of all, we notice that K : L α , β 2 L α , β 2 and K : L α , β 2 C ˜ η , ζ are well-defined and compact (cf. Corollary 1 and Lemma 4). Choose η 0 η , 1 + α 2 and ζ 0 ζ , 1 + β 2 . By using Corollary 6 we get
L ˜ n γ , δ K n L n L n K L n L L α , β 2 L ˜ n γ , δ K n L n L ˜ n γ , δ K L n L L α , β 2 + L ˜ n γ , δ K L n L n K L n L L α , β 2 C K n L n K L n L L α , β 2 , C η 0 , ζ 0 + L ˜ n γ , δ K L n L n K L n L L α , β 2 .
Since, due to Lemma 17, L ˜ n γ , δ K L n L n K L n L L α , β 2 0 holds true, it suffices to verify the convergence
K n L n K L n L L α , β 2 , C η 0 , ζ 0 0 .
We define r , s 0 , 1 by
r = 0 : 1 2 < ρ γ < 1 2 , 1 : 1 2 < ρ γ < 3 2 , a n d s = 0 : 1 2 < τ δ < 1 2 , 1 : 1 2 < τ δ < 3 2 .
Recall that the application of the operator K n to a function u C ( 1 , 1 ) can be written as (see (41))
K n u ( x ) = 1 1 L n γ , δ v γ , δ K ( x , · ) u ( y ) v γ , δ ( y ) d y .
We define
K ˜ n u ( x ) = 1 1 L n γ , δ v γ , δ K ( x , · ) ϑ v r , s ( y ) ϑ 1 u ( y ) v γ + r , δ + s ( y ) d y = 1 1 L n γ , δ K ( x , · ) v ρ α 2 γ r , τ β 2 δ s ( y ) ϑ 1 u ( y ) v γ + r , δ + s ( y ) d y .
Let u n = ϑ p n im L n , i.e., p n P n . Due to the algebraic accuracy of the Gaussian rule, in case of r + s 1 as well as in case of deg p n < n 1 and r = s = 1 , we have
K ˜ n u n ( x ) = 1 1 L n γ , δ v γ , δ K ( x , · ) ϑ v r , s ( y ) v r , s p n ( y ) v γ , δ ( y ) d y = K n u n ( x ) .
In case of deg p n = n 1 and r = s = 1 , we write p n = ε n p n 1 γ + 1 , δ + 1 + p n 1 with deg p n 1 < n 1 , i.e., p n 1 = S n 1 γ + 1 , δ + 1 p n and ε n p n 1 γ + 1 , δ + 1 = S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 p n . We get, due to the previous considerations, K ˜ n u n = ε n K ˜ n ϑ p n 1 γ + 1 , δ + 1 + K n ϑ p n 1 . Moreover,
K ˜ n ϑ p n 1 γ + 1 , δ + 1 ( x ) = 1 1 L n γ , δ v γ , δ K ( x , · ) ϑ v 1 , 1 ( y ) ( 1 y 2 ) p n 1 γ + 1 , δ + 1 ( y ) v γ , δ ( y ) d y = A n ( γ , δ ) 1 1 L n γ , δ v γ , δ K ( x , · ) ϑ v 1 , 1 ( y ) p n 1 γ , δ ( y ) v γ , δ ( y ) d y = A n ( γ , δ ) k = 1 n λ k n γ , δ v γ , δ ( x k n γ , δ ) K ( x , x k n γ , δ ) ϑ ( x k n γ , δ ) p n 1 γ , δ ( x k n γ , δ ) 1 x k n γ , δ 2 = A n ( γ , δ ) c n ( γ , δ ) k = 1 n λ k n γ , δ v γ , δ ( x k n γ , δ ) K ( x , x k n γ , δ ) ϑ ( x k n γ , δ ) p n 1 γ + 1 , δ + 1 ( x k n γ , δ ) = A n ( γ , δ ) c n ( γ , δ ) K n ϑ p n 1 γ + 1 , δ + 1 ( x ) = κ n K n ϑ p n 1 γ + 1 , δ + 1 ( x )
with κ n = n + γ + δ + 1 2 n + γ + δ + 1 , where we took into account relations (21) and (25) as well as the orthogonality properties of p n γ , δ ( x ) and p n + 1 γ , δ ( x ) . Consequently, for u L α , β 2 , we have K ˜ n L ˜ n u = K n L n u , where
L ˜ n = L n : r + s 1 , ϑ 1 κ n S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 + S n 1 γ + 1 , δ + 1 ϑ 1 L n : r = s = 1 .
We show, that L ˜ n as well as L ˜ n * converge strongly in L α , β 2 to the identity operator. For the convergence of L ˜ n , it suffices to show that, in case r = s = 1 , ϑ S n γ + 1 , δ + 1 ϑ 1 I I strongly in L α , β 2 , which is equivalent to S n γ + 1 , δ + 1 I strongly in L ρ , τ 2 , being again equivalent to
1 2 < ρ γ < 3 2 a n d 1 2 < τ δ < 3 2 ,
by Lemma 18. In the present situation, the last conditions are equivalent to the conditions which are satisfied in case r = s = 1 . Since, for f , g L α , β 2 ,
ϑ S n γ , δ ϑ 1 f , g α , β = S n γ , δ ϑ 1 f , ϑ v α γ , β δ g γ , δ = f , ϑ 1 v γ α , δ β S n γ , δ ϑ v α γ , β δ g α , β = f , v γ ρ + α 2 , δ τ + β 2 S n γ , δ v ρ + α 2 γ , τ + β 2 δ g α , β ,
we have to check if
ϑ S n γ + 1 , δ + 1 ϑ 1 I * = v γ + 1 ρ + α 2 , δ + 1 τ + β 2 S n γ + 1 , δ + 1 v ρ + α 2 γ 1 , τ + β 2 δ 1 I I i n L α , β 2
strongly, which is equivalent to S n γ + 1 , δ + 1 I strongly in L 2 γ + 2 ρ , 2 δ + 2 τ 2 . Due to Lemma 18, this is again equivalent to (60). As a consequence of these considerations we have that K L ˜ n K L n L L α , β 2 , C η 0 , ζ 0 0 if n . Hence, in order to prove (58) it suffices to show that
lim n K ˜ n K L ˜ n L L α , β 2 , C η 0 , ζ 0 = 0 .
We have
| ( K ˜ n K ) L ˜ n u ( x ) | 1 1 v γ + r ρ 2 , δ + s τ 2 ( y ) L n γ , δ K ( x , · ) v ρ α 2 γ r , τ β 2 δ s ( y ) K ( x , y ) v ρ α 2 γ r , τ β 2 δ s ( y ) · · v α 2 , β 2 ( y ) L ˜ n u ( y ) d y C L n γ , δ I K ( x , · ) v ρ α 2 γ r , τ β 2 δ s 2 γ + 2 r ρ , 2 δ + 2 s τ u α , β
Hence, in order to show relation (61) we can try to apply Lemma 16 for α 0 = 2 γ + 2 r ρ and β 0 = 2 δ + 2 s τ instead of α and β , respectively, as well as for
f ( x , y ) = v η 0 , ζ 0 ( x ) K ( x , y ) v ρ α 2 γ r , τ β 2 δ s ( y ) .
Since ( ρ , γ ) Ω and
α 0 = ρ + 2 [ r ( ρ γ ) ] 2 γ ρ : r = 0 , ρ : r = 1 , 1 2 < ρ γ < 1 , ρ 1 : r = 1 , 1 ρ γ < 3 2 ,
we have α 0 > 1 and, analogously, β 0 > 1 . Moreover,
α 0 γ = γ ρ 1 2 , 1 2 : r = 0 , 2 ( ρ γ ) 1 2 , 3 2 : r = 1 ,
such that α 0 γ and, analogously, β 0 δ belong to the interval 1 2 , 3 2 . The conditions
ψ 0 : = ψ + γ + r ρ α 2 < 1 + α 0 2 and χ 0 : = χ + δ + s τ β 2 < 1 + β 0 2
are equivalent to ψ < 1 + α 2 and χ < 1 β 2 . This all together implies that the function
f ( x , y ) v ψ 0 , χ 0 ( y ) = v η 0 , ζ 0 ( x ) K ( x , y ) v ψ 0 + ρ α 2 γ r , χ 0 + τ β 2 δ s ( y ) = v η 0 , ζ 0 ( x ) K ( x , y ) v ψ , χ ( y )
is continuous on [ 1 , 1 ] 2 and that Lemma 16 is applicable to f ( x , y ) with ψ 0 and χ 0 instead of ψ and χ . Hence,
K ˜ n K L ˜ n L L α , β 2 , C η 0 , ζ 0 sup L n γ , δ I v η 0 , ζ 0 ( x ) K ( x , · ) v ρ α 2 γ r , τ β 2 δ s 2 γ + 2 r ρ , 2 δ + 2 s τ : 1 x 1 0 i f n .
For the proof of (57) it remains to refer to (49). □
Corollary 9. 
Let α , β ( 1 , 1 ) and p , q R with p < 1 + β 2 < q , as well as H C ( R + ) such that
lim t 0 t p H ( t ) = 0 a n d lim t t q H ( t ) = 0 .
Moreover, let ξ R with 1 < α + 2 ξ < 1 and χ 1 , χ 2 : [ 1 , 1 ] [ 0 , 1 ] be continuous functions which vanish in a neighbourhood of the point 1 and are identically 1 in a neighbourhood of the point 1 . Then, for ( ρ , γ ) , ( τ , δ ) Ω , we have
L ˜ n γ , δ v ξ , 0 M n , H χ 1 M n , H χ 2 v ξ , 0 L n J ,
where the limit operators are given by
W 1 L ˜ n γ , δ v ξ , 0 M n , H χ 1 M n , H χ 2 v ξ , 0 L n = v ξ , 0 M H χ 1 M H χ 2 v ξ , 0 I
and
W t L ˜ n γ , δ v ξ , 0 M n , H χ 1 M n , H χ 2 v ξ , 0 L n = 0 , t 2 , 3 .
Proof. 
Due to our assumptions, the following functions are continuous on [ 1 , 1 ] 2 ,
( x , y ) v ζ + ξ , 0 ( x ) v χ ξ , 0 ( y ) H 1 + x 1 + y 1 1 + y 1 χ 1 ( x ) 1 χ 2 ( y ) ,
( x , y ) v ζ + ξ , η ( x ) v χ ξ , 0 ( y ) H 1 + x 1 + y 1 1 + y χ 1 ( x ) 1 χ 2 ( y ) ,
( x , y ) v ζ + ξ , 0 ( x ) v χ ξ , ψ ( y ) H 1 + x 1 + y 1 1 + y 1 χ 1 ( x ) χ 2 ( y ) ,
where ζ = max ξ , 0 , χ = max ξ , 0 , and η : = max p , 0 , ψ : = max 1 q , 0 . Since
L ˜ n γ , δ v ξ , 0 M n , H χ 1 M n , H χ 2 v ξ , 0 L n = L ˜ n γ , δ v ξ , 0 ( 1 χ 1 ) M n , H ( 1 χ 2 ) v ξ , 0 L n + L ˜ n γ , δ v ξ , 0 χ 1 M n , H ( 1 χ 2 ) v ξ , 0 L n + L ˜ n γ , δ v ξ , 0 ( 1 χ 1 ) M n , H χ 2 v ξ , 0 L n
and
0 ζ < 1 + α 2 , 0 η < 1 + β 2 , a n d 0 χ < 1 α 2 , 0 ψ < 1 β 2 ,
we are able to apply Lemma 19. □
Analogously we can prove the following corollary.
Corollary 10. 
Let α , β ( 1 , 1 ) and p , q R with p < 1 + α 2 < q , as well as H C ( R + ) such that
lim t 0 t p H ( t ) = 0 a n d lim t t q H ( t ) = 0 .
Moreover, let ξ R with 1 < β + 2 ξ < 1 and χ 1 , χ 2 : [ 1 , 1 ] [ 0 , 1 ] be continuous functions which vanish in a neighbourhood of the point 1 and are identically 1 in a neighbourhood of the point 1 . Then, for ( ρ , γ ) , ( τ , δ ) Ω , we have
L ˜ n γ , δ v 0 , ξ M n , H + χ 1 M n , H + χ 2 v 0 , ξ L n J ,
where the limit operators are given by
W 1 L ˜ n γ , δ v 0 , ξ M n , H + χ 1 M n , H + χ 2 v 0 , ξ L n = v 0 , ξ M H + χ 1 M H + χ 2 v 0 , ξ I
and
W t L ˜ n γ , δ v 0 , ξ M n , H + χ 1 M n , H + χ 2 v 0 , ξ L n = 0 , t 2 , 3 .
In what follows we identify an element ( ξ 0 , , ξ n 1 , 0 , ) im P n with the respective element ξ k k = 0 n 1 C n and the linear operator A n : im P n im P n with its matrix representation A n = a j k j , k = 1 n C n × n , i.e.,
A n ξ k k = 0 n 1 = k = 0 n 1 a j + 1 , k + 1 ξ k j = 0 n 1 .
For example, we have the representation
V n L ˜ n γ , δ K n L n V n 1 = λ k n γ , δ v γ , δ ( x k n γ , δ ) v α 2 + 1 4 , β 2 + 1 4 ( x j n γ , δ ) K x j n γ , δ , x k n γ , δ v α 2 1 4 , β 2 1 4 ( x k n γ , δ ) j , k = 1 n .
Let us formulate the following condition for a positive kernel function.
(B)
For the function H : R + R + , there are a positive constant c M and a real number κ such that
(a)
0 < t s implies H ( t ) t κ c M H ( s ) s κ
or
(b)
0 < t s implies H ( t ) t κ c M H ( s ) s κ .
Corollary 11. 
Let α , β ( 1 , 1 ) and H C ( R + ) be a positive function, which satisfies condition ( A 1 ) for ξ = 1 + β 2 and condition (B). Then, for A n = L ˜ n γ , δ M n , H L n the sequences of operators E n ( t ) A n E n ( t ) 1 L n ( t ) : X t X t ,   t 1 , 2 , 3 , n N , are uniformly bounded in case of ( ρ , γ ) , ( τ , δ ) Ω .
Proof. 
Due to Lemma 12 it suffices to prove the assertion for t = 2 . Furthermore, in view of Corollary 9 in combination with Lemma 5,(b), we have only to prove the uniform boundedness of the operators H n χ : im P n im P n with (cf. also (62))
H n χ : = V n L ˜ n γ , δ χ M n , H χ L n V n 1 = χ ( x j n γ , δ ) χ ( x k n γ , δ ) λ k n γ , δ v α 2 + 1 4 , β 2 + 1 4 ( x j n γ , δ ) v γ , δ ( x k n γ , δ ) H 1 + x j n γ , δ 1 + x k n γ , δ v α 2 1 4 , β 2 1 4 ( x k n γ , δ ) 1 + x k n γ , δ j , k = 1 n ,
where χ : [ 1 , 1 ] [ 0 , 1 ] is a continuous function, which vanishes in a neighbourhood of the point 1 and is identically 1 in a neighbourhood of the point 1 . We use the notation H n χ = h j k n , χ j , k = 1 n . Let us note that the entries h j k n , χ of the matrix H n χ are positive and, in view of the choice of the function χ : [ 1 , 1 ] [ 0 , 1 ] and (47),
h j k n , χ C χ ( 1 + x j n γ , δ ) β 2 + 1 4 n ( 1 + x k n γ , δ ) β 2 + 3 4 H 1 + x j n γ , δ 1 + x k n γ , δ
with a constant C χ C χ ( j , k , n ) . Due to (31), there is a constant C 1 C 1 ( k , n ) > 0 , such that
C 1 1 n + 1 k n 2 1 + x k n γ , δ C 1 n + 1 k n 2 .
We conclude, by additionally using condition (B),
h j k n , χ C n n n + 1 k n + 1 j n + 1 k β + 1 2 1 + x j n γ , δ 1 + x k n γ , δ κ H 1 + x j n γ , δ 1 + x k n γ , δ 1 + x j n γ , δ 1 + x k n γ , δ κ C a n + 1 k n + 1 j n + 1 k β + 1 2 H C 1 2 ( n + 1 j ) 2 ( n + 1 k ) 2 = : h n + 1 j , n + 1 k a
if condition (B),(a) is in force, and
h j k n , χ C b n + 1 k n + 1 j n + 1 k β + 1 2 H 1 C 1 2 ( n + 1 j ) 2 ( n + 1 k ) 2 = : h n + 1 j , n + 1 k b
if condition (B),(b) is fulfilled, where C d C d ( j , k , n ) ,   d a , b . As a consequence of these estimates we have
H n χ L ( 2 ) J n H d J n L ( 2 ) ,
where H d = h j k d j , k = 1 n , d a , b , and J n = δ n 1 j , k j , k = 0 n 1 . For d a , b , consider the function g d : R + R with
g d ( x ) = c M a x β + 1 2 H C 1 2 x 2 : d = a , c M b x β + 1 2 H C 1 2 x 2 : d = b .
Due to our assumptions, there exist p , q R with p < q such that H L 2 p 1 2 ( R + ) L 2 q 1 2 ( R + ) and ξ = 1 + β 2 ( p , q ) . From that, we derive
g d L 2 ( 2 p β 1 2 ) 1 2 ( R + ) L 2 ( 2 q β 1 2 ) 1 2 ( R + )
and 1 2 2 p β 1 2 , 2 q β 1 2 . Moreover, the Mellin transforms of g a ( x ) and g b ( x ) are up to a multiplicative constant equal to C 1 z H ^ z + β 2 + 1 4 and C 1 z H ^ z + β 2 + 1 4 , respectively. Furthermore, p < 1 2 + β 2 + 1 4 < q . Thus, g d fulfils condition ( A 1 ) for ξ = 1 2 . Because of H d = 1 k + 1 g d j + 1 k + 1 j , k = 0 the uniform boundedness of the operators H n χ : 2 2 is a consequence of Lemma 11, inequality (66) and the fact that the norm of the operators J n : 2 2 is equal to 1. □
Analogously to the previous one, we can prove the following corollary.
Corollary 12. 
Let α , β ( 1 , 1 ) and H C ( R + ) be a positive function, which satisfies condition ( A 1 ) for ξ = 1 + α 2 and condition (B). Then, for A n = L ˜ n γ , δ M n , H + L n , the sequences of operators E n ( t ) A n E n ( t ) 1 L n ( t ) : X t X t ,   t 1 , 2 , 3 , n N , are uniformly bounded in case of ( ρ , γ ) , ( τ , δ ) Ω .
Remark 1. 
Let | γ | = | δ | = 1 2 . Then the assertions of Lemmas 11 and 12 remain true, if we only assume that H C ( R + ) satisfies condition ( A 1 ) for σ = 1 + β 2 resp. ξ = 1 + α 2 . Thus, we do not need both the positivity of H and condition (B).
Proof. 
First of all, we notice that only the verification of (66) is necessary. Moreover, it is well known that
θ j n 1 2 , 1 2 = j 1 2 n π , θ j n 1 2 , 1 2 = j n + 1 π , θ j n 1 2 , 1 2 = j 1 2 n + 1 2 π , θ j n 1 2 , 1 2 = j n + 1 2 π .
Hence, we can give up the usage of a cutting-off function χ , and the matrices H n = H n 1 have the form H n = P n H P n or J n H J n with
H = 2 k j k σ + 1 2 H j 2 k 2 j , k = 1
or
H = 2 k 1 2 j 1 2 k 1 2 σ + 1 2 H j 1 2 2 k 1 2 2 j , k = 1 ,
where σ { α , β } . Finally, we have only to recall Lemma 11. □
Let s > 1 and J s be the Bessel function of the first kind and of order s . We have
J s ( x ) = r = 0 ( 1 ) r r ! x 2 2 r + s Γ ( s + r + 1 ) , x > 0 .
It is well-known, that J s has countable infinitely many positive simple zeros, which accumulate only at infinity. We denote these zeros in increasing order by ψ s , k , k = 1 , 2 , By using Legendre’s duplication formula we get
J 1 2 ( x ) = 2 π x sin x and J 1 2 ( x ) = 2 π x cos x ,
such that
ψ 1 2 , k = k π and ψ 1 2 , k = ( k 1 2 ) π .
Lemma 20 
([37], Theorem 4.1). For k N fixed, the nodes x k n α , β = cos θ k n α , β ( 0 < θ k n α , β < π ) of the Jacobi polynomial p n α , β ( x ) admit the representation
θ k n α , β = ψ α , k c n α , β + O n 5 ,
where
c n α , β = n + α + β + 1 2 2 + 1 α 2 3 β 2 12 .
The relation p n α , β ( cos θ ) = ( 1 ) n p n β , α ( cos ( π θ ) ) yields θ k n α , β = π θ n k + 1 , n β , α . Hence,
x k n α , β = x n k + 1 , n β , α .
Corollary 13. 
For k N fixed, we have
θ n k + 1 , n α , β = π ψ β , k c n β , α + O n 5 .
Lemma 21. 
Let j , k N , and ζ 1 , ζ 2 be real numbers which satisfy ζ 1 + ζ 2 = 1 . Moreover, let H C ( R + ) . Then, for n tending to infinity, we have
v ζ 1 2 , 0 ( x j n α , β ) v ζ 2 2 , 0 ( x k n α , β ) n ( 1 x k n α , β ) H 1 x j n α , β 1 x k n α , β 2 ( ψ α , j ) ζ 1 ( ψ α , k ) 2 ζ 2 H ψ α , j ψ α , k 2
and
v 0 , ζ 1 2 ( x n j + 1 , n α , β ) v 0 , ζ 2 2 ( x n k + 1 , n τ + , τ ) n ( 1 + x n k + 1 , n α , β ) H 1 + x n j + 1 , n α , β 1 + x n k + 1 , n α , β 2 ( ψ β , j ) ζ 1 ( ψ β , k ) 2 ζ 2 H ψ β , j ψ β , k 2 .
Proof. 
The first relation is a consequence of
v ζ 1 2 , 0 ( x j n α , β ) v ζ 2 2 , 0 ( x k n α , β ) n ( 1 x k n α , β ) H 1 x j n α , β 1 x k n α , β = 2 sin ζ 1 ( θ j n α , β / 2 ) ( θ j n α , β / 2 ) ζ 1 sin ζ 2 ( θ k n α , β / 2 ) ( θ k n α , β / 2 ) ζ 2 ( θ k n α , β / 2 ) 2 sin 2 ( θ k n α , β / 2 ) ( n θ j n α , β ) ζ 1 ( n θ k n α , β ) 2 ζ 2 · · H sin 2 ( θ j n α , β / 2 ) ( θ j n α , β / 2 ) 2 ( θ k n α , β / 2 ) 2 sin 2 ( θ k n α , β / 2 ) ( n θ j n α , β ) 2 ( n θ k n α , β ) 2
and Lemma 20. By applying (67) we get the second one. □
For an arbitrary a PC , let us compute the sequence of the adjoints of the operators L ˜ n γ , δ a L n : L α , β 2 L α , β 2 . We define integers r , s 0 , 1 as in (59) and obtain, for functions f , g L α , β 2 ,
L ˜ n γ , δ a L n f , g α , β = L n γ , δ ϑ 1 a L n f , ϑ 1 g ρ , τ = L n γ , δ ϑ 1 a L n f , v ρ γ r , τ δ s ϑ 1 g γ + r , δ + s = L n γ , δ ϑ 1 a L n f , S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g γ + r , δ + s = L n γ , δ ϑ 1 a L n f , v r , s S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g γ , δ .
In case of r + s 1 , we can proceed as follows using the algebraic accuracy of the Gaussian rule
L ˜ n γ , δ a L n f , g α , β = k = 1 n λ k n γ , δ ϑ 1 a L n f ( x k n γ , δ ) v r , s ( x k n γ , δ ) S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g ( x k n γ , δ ) ¯ = k = 1 n λ k n γ , δ ϑ 1 L n f ( x k n γ , δ ) v r , s ( x k n γ , δ ) a ¯ S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g ( x k n γ , δ ) ¯ = ϑ 1 L n f , v r , s L n γ , δ a ¯ S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g γ , δ = L n f , ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g α , β = f , L n ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 g α , β .
Hence,
L ˜ n γ , δ a L n * = L n ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ S n γ + r , δ + s v ρ γ r , τ δ s ϑ 1 I .
If r = s = 1 we use relations (21) and (25) as well as the orthogonality properties of the polynomials p n γ , δ ( x ) and p n + 1 γ , δ ( x ) . With g 0 = v ρ γ 1 , τ δ 1 ϑ 1 g we get
L ˜ n γ , δ a L n f , g α , β = L n γ , δ ϑ 1 a L n f , v 1 , 1 S n γ + 1 , δ + 1 g 0 γ , δ = L n γ , δ ϑ 1 a L n f , v 1 , 1 S n 1 γ + 1 , δ + 1 g 0 γ , δ + L n γ , δ ϑ 1 a L n f , v 1 , 1 g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ + 1 , δ + 1 γ , δ = f , L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 g α , β + A n ( γ , δ ) L n γ , δ ϑ 1 a L n f , g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ , δ γ , δ
and
A n ( γ , δ ) L n γ , δ ϑ 1 a L n f , g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ , δ γ , δ = A n ( γ , δ ) k = 1 n λ k n γ , δ ϑ 1 a L n f ( x k n γ , δ ) g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ , δ ( x k n γ , δ ) ¯ = A n ( γ , δ ) c n ( γ , δ ) k = 1 n λ k n γ , δ ϑ 1 L n f ( x k n γ , δ ) · · a ¯ g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 v 1 , 1 ( x k n γ , δ ) p n 1 γ + 1 , δ + 1 ( x k n γ , δ ) ¯ = A n ( γ , δ ) c n ( γ , δ ) ϑ 1 L n f , L n γ , δ a ¯ v 1 , 1 g 0 , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ + 1 , δ + 1 γ , δ = κ n f , L n ϑ v γ ρ , δ τ L n γ , δ a ¯ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 g 0 α , β = κ n f , L n ϑ v γ ρ , δ τ L n γ , δ a ¯ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 g α , β ,
where κ n = n + γ + δ + 1 2 n + γ + δ + 1 . Consequently, in case r = s = 1 we have
L ˜ n γ , δ a L n * = L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I + κ n L n ϑ v γ ρ , δ τ L n γ , δ a ¯ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I .
It is easy to see that L ˜ n γ , δ a L n converges strongly to a I in L α , β 2 (see the beginning of the proof of Lemma 22 below). Let us discuss the convergence of the adjoint operators. At first, we consider the operators on the right hand side of (68). For γ 0 = γ + r and δ 0 = δ + s , the conditions (55) and (56) are equivalent to 1 2 < γ ρ + 2 r δ τ + 2 s < 3 2 as well as 1 2 < ρ γ r , τ δ s < 1 2 ,   ρ < 2 γ + 1 ,   τ < 2 δ + 1 , respectively. Hence, by definition (59) of the numbers r , s 0 , 1 , we can apply Corollary 8 together with the strong convergence of L n I in L α , β 2 and get the strong convergence
L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I a ¯ I i n L α , β 2
for r , s 0 , 1 . In particular, we have the strong convergence of L ˜ n γ , δ a L n * for r + s 1 . From (70) we also get the strong convergence of
L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I = L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I · · ϑ v γ + 1 ρ , δ + 1 τ L n 1 γ , δ S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I a ¯ I .
Thus, due to formula (69), to prove the strong convergence of the operators L ˜ n γ , δ a L n * in case r = s = 1 , it remains to show that the operators
C n : = ϑ v γ ρ , δ τ L n γ , δ a ¯ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I
converge in L α , β 2 strongly to the zero operator. Up to now we know that the operators L ˜ n γ , δ a L n * and L n ϑ v γ + 1 ρ , δ + 1 τ L n γ , δ a ¯ S n 1 γ + 1 , δ + 1 v ρ γ 1 , τ δ 1 ϑ 1 I are uniformly bounded. Thus, in view of relation (69), also the operators C n : L α , β 2 L α , β 2 are uniformly bounded, and it suffices to show their convergence on a dense subset of L α , β 2 . Such a subset is the space v γ + 1 ρ , δ + 1 τ ϑ P : P P because of the relations
f v γ + 1 ρ , δ + 1 τ ϑ P α , β = v ρ γ 1 , τ δ 1 ϑ 1 f P 2 γ + 2 ρ , 2 δ + 2 τ f L α , β 2 ,
2 γ + 2 ρ > 1 , 2 δ + 2 τ > 1 (note that ρ γ , τ δ < 3 2 ), and the density of P in L α , β 2 . These results can be used for the proof of the following lemma. Nevertheless, we will give a shorter proof of the strong convergence of the adjoint operators.
Lemma 22. 
Let a PC and ( ρ , γ ) , ( τ , δ ) Ω . Then we have the strong convergences
L ˜ n γ , δ a L n a I a n d L ˜ n γ , δ a L n * a ¯ I
in the space L α , β 2 .
Proof. 
Since L n = ϑ S n ρ , τ ϑ 1 I and consequently L ˜ n γ , δ a L n = ϑ L n γ , δ a S n ρ , τ ϑ 1 I , The strong convergence of L ˜ n γ , δ a L n to a I in L α , β 2 is equivalent to the strong convergence of L n γ , δ a S n ρ , τ to a I in L ρ , τ 2 . To prove this, we can apply Corollary 7 for α = γ 0 = ρ and β = δ 0 = τ .
Let us turn to the convergence of the adjoint operators. We define integers r , s 0 , 1 as in (59). Since we already know that these operators L ˜ n γ , δ a L n * : L α , β 2 L α , β 2 are uniformly bounded it suffices to prove their convergence on a dense subset of L α , β 2 . As such a subset we can take the set
X r s = ϑ v γ + r ρ , δ + s τ P : P P ,
since 2 γ + 2 r ρ , 2 δ + 2 s τ > 1 for ( ρ , γ ) , ( τ , δ ) Ω and P is dense in L 2 γ + 2 r ρ , 2 δ + 2 s τ 2 as well as
f ϑ v γ + r ρ , δ + s τ P α , β = ϑ 1 v ρ γ r , τ δ s f P 2 γ + 2 r ρ , 2 δ + 2 s τ f L α , β 2 .
For f L α , β 2 , g = ϑ v γ + r ρ , δ + s τ P with P P , and all sufficiently large n , we obtain
L ˜ n γ , δ a L n f , g α , β = L n γ , δ ϑ 1 a L n f , ϑ 1 g ρ , τ = L n γ , δ ϑ 1 a L n f , v γ + r ρ , δ + s τ P ρ , τ = L n γ , δ ϑ 1 a L n f , v r , s P γ , δ = k = 1 n λ k n γ , δ ϑ 1 a L n f ( x k n γ , δ ) v r , s ( x k n γ , δ ) P ( x k n γ , δ ) ¯ = k = 1 n λ k n γ , δ ϑ 1 L n f ( x k n γ , δ ) v r , s ( x k n γ , δ ) a ¯ P ( x k n γ , δ ) ¯ = v r , s ϑ 1 L n f , L n γ , δ a ¯ P γ , δ = L n f , ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ P α , β = f , L n ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ P α , β .
It remains to show that ϑ v γ + r ρ , δ + s τ L n γ , δ a ¯ P converges to a ¯ g in L α , β 2 , which is equivalent to L n γ , δ a ¯ P a ¯ P in L 2 γ + 2 r ρ , 2 δ + 2 s τ 2 . This follows from Lemma 15 by choosing the parameters α : = 2 γ + 2 r ρ , β : = 2 δ + 2 s τ , as well as ψ = χ = 0 , and taking into account that the conditions
1 2 < γ + 2 r ρ , δ + 2 s τ < 3 2 a n d 2 γ + 2 r ρ , 2 δ + 2 s > 1
are satisfied for ( ρ , γ ) , ( τ , δ ) Ω by the definition of r , s 0 , 1 .
For H C ( R + ) , we define the integral operator
( N H f ) ( x ) = 1 1 H 1 + y 1 + x f ( y ) d y 1 + x , x ( 1 , 1 ) .
We have the following relation
N H = M H 1 w i t h H 1 ( t ) = H ( t 1 ) t 1 , t R + .
Furthermore, one can easily show that the adjoint operator of M H : L α , β 2 L α , β 2 is equal to M H * = v α , β N H ¯ v α , β I .
Lemma 23. 
Let α , β ( 1 , 1 ) , p < 1 + β 2 < q and H C 1 ( R + ) satisfy H L 2 p 1 2 L 2 q 1 2 as well as H L 2 p + 1 2 L 2 q + 1 2 . If the operators L ˜ n γ , δ M n , H L n : L α , β 2 L α , β 2 are uniformly bounded for some γ , δ with ( ρ , γ ) , ( τ , δ ) Ω , then we have the strong convergences
L ˜ n γ , δ M n , H L n M H a n d L ˜ n γ , δ M n , H L n * M H * i n L α , β 2 .
Proof. 
Due to the Banach-Steinhaus theorem, it suffices to show the convergence on a dense subset of L α , β 2 . We consider functions g ( x ) of the form
g ( x ) = ϑ ( x ) ( 1 x 2 ) 2 p ( x ) , x ( 1 , 1 ) ,
where p ( x ) is an arbitrary polynomial. By choosing χ p , 1 + β 2 0 , 1 + β 2 and taking Corollary 6 into account, we get sup L ˜ n γ , δ L ( C 0 , χ , L α , β 2 ) : n N < . Consequently, since g im L n ,
L ˜ n γ , δ M n , H L n g M H g α , β = L ˜ n γ , δ M n , H g M H g + L ˜ n γ , δ M H g M H g α , β C M n , H g M H g 0 , χ , + L ˜ n γ , δ M H g M H g α , β .
At first, we consider the term M n , H g M H g 0 , χ , . For f : ( 1 , 1 ) C , we denote by R n γ , δ ( f ) the error
R n γ , δ ( f ) = 1 1 f ( y ) v γ , δ ( y ) d y k = 1 n λ k n γ , δ f ( x k n γ , δ ) .
From [38] ((5.1.35)) follows
R n γ , δ ( f ) C n 1 1 | f ( y ) | v γ , δ ( y ) φ ( y ) d y ,
where φ ( x ) = 1 x 2 and the constant does not depend on f and n . We have
( M n , H g ) ( x ) ( M H g ) ( x ) = R n γ , δ ( f x ) , x ( 1 , 1 )
with
f x ( y ) = H 1 + x 1 + y g 0 ( y ) 1 + y a n d g 0 ( y ) = v γ , δ ( y ) g ( y ) .
With the help of Lemma 3, we get
1 1 | f x ( y ) | v γ , δ ( y ) φ ( y ) d y 1 1 H 1 + x 1 + y φ ( y ) | g ( y ) | ( 1 + y ) 2 d y + 1 1 H 1 + x 1 + y ( 1 + x ) φ ( y ) | g ( y ) | ( 1 + y ) 3 d y + 1 1 H 1 + x 1 + y v γ , δ ( y ) φ ( y ) | g 0 ( y ) | 1 + y d y C 1 1 H 1 + x 1 + y d y 1 + y + 1 1 H 1 + x 1 + y 1 + x 1 + y d y 1 + y C ( 1 + x ) χ , x ( 1 , 1 ) .
Thus,
1 n sup v 0 , χ ( x ) 1 1 | f x ( y ) | v γ , δ ( y ) 1 y 2 d y : 1 < x < 1 0 , n ,
and (74) delivers lim n M n , H g M H g 0 , χ , = 0 .
Secondly, we deal with L ˜ n γ , δ M H g M H g α , β . Again by Lemma 3 we can estimate
M H g ( x ) C 1 1 H 1 + x 1 + y d y 1 + y C ( 1 + x ) χ
and get M H g C ˜ 0 , χ b , which implies, due to Corollary 6,
lim n L ˜ n γ , δ M H g M H g α , β = 0 .
Let us turn to the strong convergence of the adjoint operators. Let f L α , β 2 , define r , s 0 , 1 as in (59), and take g from the dense subset X r s defined in (71), what means g = v γ + r ρ , δ + s τ ϑ P with P P . We set
h ( x , y ) : = H 1 + x 1 + y 1 1 + y , x , y ( 1 , 1 ) ,
and get, for all sufficiently large n ,
L ˜ n γ , δ M n , H L n f , g α , β = L n γ , δ ϑ 1 M n , H L n f , v r , s P γ , δ = j = 1 n λ j n γ , δ ϑ 1 ( x j n γ , δ ) M n , H L n f ( x j n γ , δ ) v r , s P ( x j n γ , δ ) ¯ = j = 1 n λ j n γ , δ ϑ 1 ( x j n γ , δ ) k = 1 n λ k n γ , δ h ( x j n γ , δ , x k n γ , δ ) v γ , δ L n f ( x k n γ , δ ) v r , s P ( x j n γ , δ ) ¯ = k = 1 n λ k n γ , δ v γ , δ L n f ( x k n γ , δ ) j = 1 n λ j n γ , δ h ( x j n γ , δ , x k n γ , δ ) ¯ v r , s ϑ 1 P ( x j n γ , δ ) ¯ = k = 1 n λ k n γ , δ v γ , δ L n f ( x k n γ , δ ) j = 1 n λ j n γ , δ h ( x j n γ , δ , x k n γ , δ ) ¯ v γ , δ v α , β g ( x j n γ , δ ) ¯ = k = 1 n λ k n γ , δ v γ , δ L n f ( x k n γ , δ ) N n , H ¯ v α , β g ( x j n γ , δ ) ¯ = k = 1 n λ k n γ , δ v r , s ϑ 1 L n f ( x k n γ , δ ) v γ r , δ s ϑ N n , H ¯ v α , β g ( x k n γ , δ ) ¯
where
N n , H : C ( 1 , 1 ) C ( 1 , 1 ) , u j = 1 n λ j n γ , δ H 1 + x j n γ , δ 1 + · ( v γ , δ u ) ( x j n γ , δ ) 1 + · .
In case of r + s 1 , we use the algebraic accuracy of the Gaussian rule and obtain
L ˜ n γ , δ M n , H L n f , g α , β = v r , s ϑ 1 L n f , L n γ , δ v γ r , δ s ϑ N n , H ¯ v α , β g γ , δ = ϑ 2 L n f , v r , s ϑ L n γ , δ ϑ 1 v ρ γ r , τ δ s v α , β N n , H ¯ v α , β g γ , δ = f , L n v γ + r ρ , δ + s τ L ˜ n γ , δ v ρ γ r , τ δ s v α , β N n , H ¯ v α , β g α , β .
Thus,
L ˜ n γ , δ M n , H L n * g = L n v γ + r ρ , δ + s τ L ˜ n γ , δ v ρ γ r , τ δ s v α , β N n , H ¯ v α , β g = L n v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N n , H ¯ v α , β g
if r + s 1 . In case of r = s = 1 , we write
L ˜ n γ , δ M n , H L n f , g α , β = k = 1 n λ k n γ , δ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f ( x k n γ , δ ) N ( x k n γ , δ ) ¯ + k = 1 n λ k n γ , δ v 1 , 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f ( x k n γ , δ ) N ( x k n γ , δ ) ¯ ,
where we used the abbreviation N ( x ) = v γ 1 , δ 1 ϑ N n , H ¯ v α , β g ( x ) and where, for the second sum in (77), we get
k = 1 n λ k n γ , δ v 1 , 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f ( x k n γ , δ ) v γ 1 , δ 1 ϑ N n , H ¯ v α , β g ( x k n γ , δ ) ¯ = S n 1 γ + 1 , δ + 1 ϑ 1 L n f , L n γ , δ v γ 1 , δ 1 ϑ N n , H ¯ v α , β g γ + 1 , δ + 1 = ϑ 1 L n f , v 1 , 1 S n 1 γ + 1 , δ + 1 L n γ , δ v γ 1 , δ 1 ϑ N n , H ¯ v α , β g γ , δ = ϑ 2 L n f , v 1 , 1 ϑ S n 1 γ + 1 , δ + 1 L n γ , δ ϑ 1 v ρ γ 1 , τ δ 1 v α , β N n , H ¯ v α , β g γ , δ = f , L n v γ + 1 ρ , δ + 1 τ ϑ S n 1 γ + 1 , δ + 1 L n γ , δ ϑ 1 v ρ γ 1 , τ δ 1 v α , β N n , H ¯ v α , β g α , β .
For the first sum in (77) we use the relations
v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f ( x k n γ , δ ) = ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 v 1 , 1 ( x k n γ , δ ) p n 1 γ + 1 , δ + 1 ( x k n γ , δ ) = ( 25 ) c n 1 ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ , δ ( x k n γ , δ )
and conclude with the help of (21)
k = 1 n λ k n γ , δ v 1 , 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f ( x k n γ , δ ) N ( x k n γ , δ ) ¯ = c n 1 ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 k = 1 n λ k n γ , δ p n 1 γ , δ ( x k n γ , δ ) N ( x k n γ , δ ) ¯ = c n 1 ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ , δ , L n γ , δ N γ , δ = c n 1 A n 1 ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 v 1 , 1 p n 1 γ + 1 , δ + 1 , L n γ , δ N γ , δ = κ n 1 ϑ 1 L n f , p n 1 γ + 1 , δ + 1 γ + 1 , δ + 1 p n 1 γ + 1 , δ + 1 , L n γ , δ N γ + 1 , δ + 1 = κ n 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 ϑ 1 L n f , L n γ , δ N γ + 1 , δ + 1 = κ n 1 ϑ 1 L n f , S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 L n γ , δ N γ + 1 , δ + 1 ,
which is equal to
κ n 1 f , L n v γ + 1 ρ , δ + 1 τ ϑ S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 L n γ , δ ϑ 1 v ρ γ 1 , τ δ 1 h α , β
with h = v α , β N n , H ¯ v α , β g . That means, together with (76)–(78),
L ˜ n γ , δ M n , H L n * g = L n v γ + r ρ + α 2 , δ + s τ + β 2 S ˜ n L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N n , H ¯ v α , β g ,
where
S ˜ n = I : r + s 1 , κ n 1 S n γ + 1 , δ + 1 S n 1 γ + 1 , δ + 1 + S n 1 γ + 1 , δ + 1 : r = s = 1 .
If we apply Corollary 8 with a 1 ,   γ 0 = γ + r , and δ 0 = δ + s , then we see that v γ + r ρ + α 2 , δ + s τ + β 2 S ˜ n v ρ + α 2 γ r , τ + β 2 δ s I converges in L α , β 2 strongly to the identity operator. Hence, it remains to show the convergence
v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N n , H ¯ v α , β g v α , β N H ¯ v α , β g
in L α , β 2 . For this, we can take g from the subset X r s 0 = ϑ v γ + r ρ , δ + s τ v 2 , 2 P : P P of X r s , which is also dense in L α , β 2 . At first we remark that, since the conditions of Lemma 15 are satisfied for ψ = γ + r ρ + α 2 and χ = δ + s τ + β 2 , we have, for all f R 1 + α 2 , 1 + β 2 0 ,
v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s f f i n L α , β 2 .
Choose ψ 0 1 + β q , 1 + β 2 0 , 1 + β 2 β , 1 + β 2 . With the help of (80) we can estimate
v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N n , H ¯ v α , β g v α , β N H ¯ v α , β g α , β C v α , β N n , H ¯ v α , β g v α , β N H ¯ v α , β g max 0 , α , ψ 0 , + v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N H ¯ v α , β g v α , β N H ¯ v α , β g α , β .
We have
N n , H ¯ v α , β g ( x ) N H ¯ v α , β g ( x ) = R n γ , δ ( f x ) , x ( 1 , 1 ) ,
with
f x ( y ) = H 1 + y 1 + x ¯ g 0 ( y ) d y 1 + x and g 0 ( y ) = v α γ , β δ ( y ) g ( y ) = v α ρ 2 + r + 2 , β τ 2 + s + 2 P .
Set H 2 ( t ) = H ( t 1 ) t 2 , t R + . Since 1 + β q < ψ 0 < 1 + β p holds true, we get 2 p 1 < 1 + 2 β 2 ψ 0 < 2 q 1 . Hence, due to our assumptions H L 2 p 1 2 L 2 q 1 2 and H L 2 p + 1 2 L 2 q + 1 2 , we have H L 2 ( β ψ 0 ) + 1 2 and H L 2 ( β ψ 0 ) + 3 . Consequently, H 1 , H 2 L 2 ( ψ 0 β ) 1 2 , where H 1 is defined in (72). By
α ρ 2 + γ + r + 2 + 1 2 = γ ρ 2 + α + γ 2 + r + 5 2 > 1 2 r 1 + r + 5 2 = 1 ,
β τ 2 + δ + s + 2 + 1 2 = δ τ 2 + β + δ 2 + s + 5 2 > 1 2 s 1 + s + 5 2 = 1 ,
and again using Lemma 3 we obtain
1 1 | f x ( y ) | v γ , δ ( y ) φ ( y ) d y 1 1 H 1 + y 1 + x v α , β ( y ) φ ( y ) | g ( y ) | ( 1 + x ) 2 d y + 1 1 H 1 + y 1 + x v γ , δ ( y ) φ ( y ) | g 0 ( y ) | 1 + x d y C 1 1 H 1 + y 1 + x 1 + y 1 + x d y 1 + x + 1 1 H 1 + y 1 + x d y 1 + x C 1 1 H 2 1 + x 1 + y d y 1 + y + 1 1 H 1 1 + x 1 + y d y 1 + y C ( 1 + x ) β ψ 0 , x ( 1 , 1 ) .
Thus,
1 n sup υ 0 , ψ 0 β ( x ) · 1 1 | f x ( y ) | d y : 1 < x < 1 0 , n .
Relation (74) yields
lim n v α , β N n , H ¯ v α , β g N H ¯ v α , β g max 0 , α , ψ , = 0 .
As above we see that v α , β N H v α , β g belongs to C ˜ max 0 , α , ψ 0 b , and (80) delivers
lim n v γ + r ρ + α 2 , δ + s τ + β 2 L n γ , δ v ρ + α 2 γ r , τ + β 2 δ s v α , β N H v α , β g v α , β N H v α , β g α , β = 0 .
The lemma is proved. □
The proof of the following lemma is analogous.
Lemma 24. 
Let α , β ( 1 , 1 ) , p < 1 + α 2 < q and H C 1 ( R + ) satisfy H L 2 p 1 2 L 2 q 1 2 as well as H L 2 p + 1 2 L 2 q + 1 2 . If the operators L ˜ n γ , δ M n , H + L n : L α , β 2 L α , β 2 are uniformly bounded for some γ , δ with ( ρ , γ ) , ( τ , δ ) Ω , then we have the strong convergences
L ˜ n γ , δ M n , H + L n M H + a n d L ˜ n γ , δ M n , H + L n * M H + * i n L α , β 2 .
The following lemma is a version of the dominated convergence theorem and will be useful in proofs of the strong convergence of operator sequences in the Hilbert space 2 .
Lemma 25. 
Let ξ , η 2 , ξ n = ( ξ j n ) j = 0 , | ξ j n | | η j | j = 0 , 1 , 2 , , n n 0 and
lim n ξ j n = ξ j j = 0 , 1 , 2 ,
Then lim n ξ n ξ 2 = 0 .
For what follows we set
b j k n , ± : = φ ( x k n γ , δ ) n v α 2 + 1 4 , β 2 + 1 4 ( x j n γ , δ ) v α 2 + 1 4 , β 2 + 1 4 ( x k n γ , δ ) H 1 x j n γ , δ 1 x k n γ , δ 1 1 x k n γ , δ .
Moreover, we need the limit relations (see Lemma 20)
lim n n 2 ( 1 x k n γ , δ ) = lim n sin 2 θ k n γ , δ 2 2 θ k n γ , δ 2 2 n 2 θ k n γ , δ 2 = ψ γ , k 2 2
and ([39] (15.3.11))
lim n n γ + 1 λ k n γ , δ = 2 δ γ + 1 2 ψ γ , k γ J γ ( ψ γ , k ) ,
which are true for all γ , δ > 1 and fixed k N .
Lemma 26. 
Let the conditions of Corollary 11 be satisfied. Then the strong limits of the operators
V n L ˜ n γ , δ M n , H L n V n 1 P n : 2 2 a n d V n L ˜ n γ , δ M n , H L n V n 1 P n * : 2 2
as well as
F n V n L ˜ n γ , δ M n , H L n V n 1 F n P n : 2 2 a n d F n V n L ˜ n γ , δ M n , H L n V n 1 F n P n * : 2 2
exist, where W 2 L ˜ n γ , δ M n , H L n = 0 and
W 3 L ˜ n γ , δ M n , H L n = ψ δ , j ψ δ , k β + 1 2 4 ψ δ , k J δ ( ψ δ , k ) 2 H ψ δ , j ψ δ , k 2 j , k = 1 .
Proof. 
Due to Corollary 11 and the uniform boundedness of the operators E n ( t ) ± 1 (see Lemma 12), all sequences of operators under consideration here are uniformly bounded. Thus, in view of the Banach-Steinhaus theorem, it suffices to verify the convergence on the set { e m = ( δ j , m ) j = 0 : m = 0 , 1 , } 2 . Moreover, in view of Corollary 9 and Lemma 5 we can replace the operators M n , H by χ M n , H χ I , where χ : [ 1 , 1 ] [ 0 , 1 ] is a continuous function, which vanishes in a neighbourhood of the point 1 and is identically 1 in a neighbourhood of the point 1 .
Fix k N . Regarding the proof of Corollary 11, for n > k , we have
F n V n L ˜ n γ , δ χ M n , H χ L n V n 1 F n P n e k 1 = h n + 1 j , n + 1 k n , χ j = 1 n = χ ( x n + 1 j , n γ , δ ) n λ n + 1 k , n γ , δ v γ + 1 2 , δ + 1 2 ( x n + 1 k , n γ , δ ) b n + 1 j , n + 1 k n , χ ( x n + 1 k , n γ , δ ) j = 1 n
with b j k n , defined in (81) and the entries h j k n , χ of the matrix H n χ defined in the proof of Corollary 11. For fixed j , Lemma 21 yields ( ζ 1 = β + 1 2 , ζ 2 = β + 1 2 )
χ ( x n + 1 j , n γ , δ ) b n + 1 j , n + 1 k n , χ ( x n + 1 k , n γ , δ ) 2 ( ψ δ , j ) β + 1 2 ( ψ δ , k ) β + 3 2 H ψ δ , j ψ δ , k 2 i f n ,
since lim n x n + 1 j , n γ , δ = 1 . Moreover, taking into account (82) and (83) we get
n λ n + 1 k , n γ , δ v γ + 1 2 , δ + 1 2 ( x n + 1 k , n γ , δ ) = n δ + 1 λ k n δ , γ 2 1 + x k n δ , γ γ + 1 2 n 2 ( 1 x k n δ , γ ) δ + 1 2 2 γ δ + 1 2 ψ δ , k δ J δ ( ψ δ , k ) 2 1 2 γ + 1 2 ψ δ , k 2 2 δ + 1 2 = 2 J δ ( ψ δ , k ) 2 ψ δ , k
if n . Due to (64) and (65), respectively, we can estimate
h n + 1 j , n + 1 k n , χ h j k d ,
where h j k d j = 1 = H d e k 1 2 , since H d L ( 2 ) (cf. the end of the proof of Corollary 11). Hence, it remains to apply Lemma 25 with ξ j n = h n j , n + 1 k n , χ and η j = h j + 1 , k d for k = 1 , 2 , to get formula (84).
On the other hand, since lim n x j n γ , δ = 1 and χ ( x ) = 0 in a neighbourhood of the point 1 , we have χ ( x j n γ , δ ) = 0 for all sufficiently large n . Hence,
V n L n γ , δ χ M n , H χ V n 1 P n e k 1 j = h j k n , χ 0 for all j N
if n tends to infinity. Moreover, again due to the choice of χ ( x ) , we have h j k n , χ = 0 for all n > n 0 = n 0 ( k ) and j = 1 , , n . Consequently, if we set
M : = max h j k n , χ : n = 1 , , n 0 , j = 1 , , n ,
then h j k n , χ f j 1 , j N , where f j = M , j = 0 , , n 0 1 and f j = 0 , j = n 0 , n 0 + 1 , , such that f = f j j = 0 2 . Thus, the application of Lemma 25 with ξ j n = h j + 1 , k n , χ and η j = f j yields W 2 L ˜ n γ , δ M n , H L n = 0 . The proof of the strong convergence of the adjoint oprators follows the same ideas by using that H d * belongs to L ( 2 ) (see (85)). □
In the same way we can prove the following.
Lemma 27. 
Let the conditions of Corollary 12 be in force. Then the strong limits of the operators
V n L ˜ n γ , δ M n , H + L n V n 1 P n : 2 2 a n d V n L ˜ n γ , δ M n , H + L n V n 1 P n * : 2 2
as well as
F n V n L ˜ n γ , δ M n , H + L n V n 1 F n P n : 2 2 a n d F n V n L ˜ n γ , δ M n , H + L n V n 1 F n P n * : 2 2
exist, where W 3 L ˜ n γ , δ M n , H + L n = 0 as well as
W 2 L ˜ n γ , δ M n , H + L n = ψ γ , j ψ γ , k α + 1 2 4 ψ γ , k J γ ( ψ γ , k ) 2 H ψ γ , j ψ γ , k 2 j , k = 1 .
Lemma 28. 
Let a PC and ( ρ , γ ) , ( τ , δ ) Ω . Then the strong limits of the operators
V n L ˜ n γ , δ a L n V n 1 P n : 2 2 a n d V n L ˜ n γ , δ a L n V n 1 P n * : 2 2
as well as
F n V n L ˜ n γ , δ a L n V n 1 F n P n : 2 2 a n d F n V n L ˜ n γ , δ a L n V n 1 F n P n * : 2 2
exist, where
W 2 ( L ˜ n γ , δ a L n ) = a ( 1 ) I a n d W 3 ( L ˜ n γ , δ a L n ) = a ( 1 ) I
as well as I : 2 2 is the identity operator.
Proof. 
We are only going to show the first two convergences. The proof of the other convergences can be done in the same way. In view of Lemmas 12 and 22 it suffices to show the convergence on a dense subset of the space 2 . Let e m = ( δ j , m ) j = 0 , m N 0 , For n > m + 1 , we have
V n L ˜ n γ , δ a L n V n 1 P n e m = a ( x m + 1 , n γ , δ ) δ j , m + 1 j = 1 n = a ( x m + 1 , n γ , δ ) e m .
Since a x m + 1 , n γ , δ a ( 1 ) as n tends to infinity, we obtain W 2 ( L ˜ n γ , δ a L n ) = a ( 1 ) I . Moreover, from (86) we infer V n L ˜ n γ , δ a L n V n 1 P n = diag a ( x k n γ , δ ) k = 1 n . Hence,
V n L ˜ n γ , δ a L n V n 1 P n * = diag a ( x k n γ , δ ) ¯ k = 1 n ,
and the strong convergence of these adjoint operators follows as before. □

6. The Stability Theorem

We recall that A n = L ˜ n γ , δ I + c M n , H + c + M n , H + + + K n L n .
Theorem 1. 
Let α , β ( 1 , 1 ) , c ± PC , and H = H ± C 1 ( R + ) be positive functions, satisfying condition ( A 1 ) for ξ = ξ ± and real numbers p = p ± , q = q ± , where ξ + = 1 + α 2 and ξ = 1 + β 2 , as well as H ± L 2 p ± + 1 2 L 2 q ± + 1 2 and condition (B), respectively. Moreover, let K : [ 1 , 1 ] × [ 1 , 1 ] C be a function, which fulfils the requirements of Lemma 19. Then A n belongs to the algebra F for all γ , δ with ( ρ , γ ) , ( τ , δ ) Ω . If this is the case, then for the sequence A n to be stable it is necessary that the operators
W 1 ( A n ) = I + c M H + c + M H + + + K : L α , β 2 L α , β 2
and
W 2 ( A n ) = I + c + ( 1 ) ψ γ , j ψ γ , k α + 1 2 4 ψ γ , k J γ ( ψ γ , k ) 2 H + ψ γ , j ψ γ , k 2 j , k = 1 : 2 2 ,
W 3 ( A n ) = I + c ( 1 ) ψ δ , j ψ δ , k β + 1 2 4 ψ δ , k J δ ( ψ δ , k ) 2 H ψ δ , j ψ δ , k 2 j , k = 1 : 2 2 .
are invertible.
Proof. 
At first let K be the zero operator. We notice that
L ˜ n γ , δ I + c M n , H + c + M n , H + + L n = L n + L ˜ n γ , δ c L n L ˜ n γ , δ M n , H L n + L ˜ n γ , δ c + L n L ˜ n γ , δ M n , H + + L n .
From Corollaries 11 and 12, Lemmas 22–24 follows
A n I + c M H + c + M H + + = W 1 ( A n ) and ( A n ) * W 1 ( A n ) * .
Lemmas 26–28 deliver the existence of
W 2 ( A n ) = I + c + ( 1 ) ψ γ , j ψ γ , k α + 1 2 4 ψ γ , k J γ ( ψ γ , k ) 2 H + ψ γ , j ψ γ , k 2 j , k = 1 ,
and
W 3 ( A n ) = I + c ( 1 ) ψ δ , j ψ δ , k β + 1 2 4 ψ δ , k J γ ( ψ δ , k ) 2 H ψ δ , j ψ δ , k 2 j , k = 1 ,
as well as the strong convergence of the sequences of the respective adjoint operators. Hence ( A n ) F . This allows us to apply Proposition 2, which immediately delivers the assertion. If the integral operator K does not vanish, then the assertion follows in combination with Lemma 19. □
In case of Chebyshev nodes, we can formulate the following theorem.
Theorem 2. 
Let α , β ( 1 , 1 ) , c ± PC , and H ± C 1 ( R + ) be positive functions, satisfying condition ( A 1 ) for ξ = ξ ± and real numbers p ± , q ± , where ξ + = 1 + α 2 , ξ = 1 + β 2 , and let H ± L 2 p ± + 1 2 L 2 q ± + 1 2 . Moreover, assume that K : [ 1 , 1 ] × [ 1 , 1 ] C fulfils the requirements of Lemma 19. Then, A n belongs to the algebra F for all γ , δ with | γ | = | δ | = 1 2 and ( ρ , γ ) , ( τ , δ ) Ω . In that case, for the sequence A n to be stable it is necessary that
(a)
the kernel of the operator W 1 ( A n ) = I + c M H + c + M H + + + K : L α , β 2 L α , β 2 is trivial,
(b)
the curves
1 + c ( 1 ) H ^ ( ξ i t ) : t R ¯ a n d 1 + c + ( 1 ) H ^ + ( ξ + i t ) : t R ¯
do not contain the point 0 and their winding numbers are equal to zero,
(c)
the kernels of the operators
W 2 ( A n ) = I + c + ( 1 ) ψ γ , j ψ γ , k α + 1 2 2 π ψ γ , k H + ψ γ , j ψ γ , k 2 j , k = 1 : 2 2 ,
W 3 ( A n ) = I + c ( 1 ) ψ δ , j ψ δ , k β + 1 2 2 π ψ δ , k H ψ δ , j ψ δ , k 2 j , k = 1 : 2 2
are trivial.
If K is the zero operator and one of the functions c ± vanishes, then condition (a) is automatically a consequence of condition (b), which due to Lemma 7.
Proof. 
In view of Remark 1, in comparison with Theorem 1 we can omit the positivity of the Mellin kernel functions H ± and condition (B), and we only have to show that the operators W 1 ( A n ) , W 2 ( A n ) , and W 3 ( A n ) are invertible. Since K : L α , β 2 L α , β 2 is compact (cf. Lemma 4), we can make use of Lemma 6. Thus, conditions (a) and (b) deliver the invertibility of the operator W 1 ( A n ) : L α , β 2 L α , β 2 . It remains to check the invertibility of W 2 ( A n ) , W 3 ( A n ) : 2 2 . Without loss of generality, we assume that γ = δ = 1 2 . Then,
W 2 ( A n ) = I + c + ( 1 ) 2 k j k α + 1 2 H + j k 2 j , k = 1 ,
W 3 ( A n ) = I + c ( 1 ) 2 k j k β + 1 2 H j k 2 j , k = 1 .
We consider the function g ± : R + R with
g d ( x ) = x α + 1 2 H + x 2 : d = + , x β + 1 2 H x 2 : d = .
Due to our assumption we have H ± L 2 p ± 1 2 ( R + ) L 2 q ± 1 2 ( R + ) and ξ ± ( p ± , q ± ) . From that, we derive
g + L 2 ( 2 p + α 1 2 ) 1 2 ( R + ) L 2 ( 2 q + α 1 2 ) 1 2 a n d 1 2 2 p + α 1 2 , 2 q + α 1 2
as well as
g L 2 ( 2 p β 1 2 ) 1 2 ( R + ) L 2 ( 2 q β 1 2 ) 1 2 a n d 1 2 2 p β 1 2 , 2 q β 1 2 .
Moreover, the Mellin transforms of g + ( x ) and g ( x ) are equal to
H ^ + z + α 2 + 1 4 a n d H ^ z + β 2 + 1 4 ,
respectively, and p + < 1 2 + α 2 + 1 4 < q + and p < 1 2 + β 2 + 1 4 < q . Thus, g ± fulfils condition ( A 1 ) for ξ = 1 2 . This allows us to apply Lemma 11. Consequently, W 2 ( A n ) , W 3 ( A n ) alg T ( PC ) with
Γ W 2 ( A n ) = 1 + c + ( 1 ) H ^ + ( ξ + + i t ) : t R ¯
and
Γ W 3 ( A n ) = 1 + c ( 1 ) H ^ ( ξ + i t ) : t R ¯ .
In view of condition (b), the curves Γ W 2 ( A n ) , Γ W 3 ( A n ) do not contain the point 0 and their winding number is zero. From Proposition 1 we derive that W 2 ( A n ) and W 3 ( A n ) are Fredholm operators with vanishing index. Thus, condition (c) delivers the invertibility of those operators. □

7. Final Remarks

Finally, let us discuss the progress we have made in the present paper for possible representations of endpoint singularities of the approximate solution (cf. (5))
u n ( x ) = ( 1 x ) ρ 0 ( 1 + x ) τ 0 p n ( x )
in comparison with the paper [7]. Recall that for our method we can choose the parameters γ , δ > 1 for the nodes x j n γ , δ , the parameters α , β ( 1 , 1 ) for the space L α , β 2 , and the parameters ρ , τ > 1 for the orthonormal system v ρ 0 , τ 0 p n ρ , τ n = 0 in L α , β 2 , where
ρ 0 = ρ α 2 and τ 0 = τ β 2 .
In [7], the case γ = δ = ρ = τ = 1 2 together with α , β ( 1 , 1 ) is considered. That means that the range for ρ 0 and τ 0 is given by the interval 1 4 , 3 4 . If, in the present paper, we choose γ = δ = 1 2 , then for the choice of ρ and τ we have to fulfil the condition
ρ , 1 2 , τ , 1 2 Ω ,
where (cf. (50))
Ω = ( ω 1 , ω 2 ) R 2 : ω 1 , ω 2 > 1 , 1 2 < ω 1 ω 2 < 1 2 , ω 1 < 2 ω 2 + 1 ( ω 1 , ω 2 ) R 2 : ω 1 , ω 2 > 1 , 1 2 < ω 1 ω 2 < 3 2 .
Condition (88) is equivalent to ρ , τ ( 0 , 1 ) ( 1 , 2 ) , such that, for ρ 0 and τ 0 we have the possible ranges
ρ 0 ρ α 2 : ρ ( 0 , 1 ) ( 1 , 2 ) , α ( 1 , 1 ) = α ( 1 , 1 ) α 2 , 1 α 2 1 α 2 , 2 α 2 = 1 2 , 3 2
and
τ 0 τ β 2 : τ ( 0 , 1 ) ( 1 , 2 ) , β ( 1 , 1 ) = β ( 1 , 1 ) β 2 , 1 β 2 1 β 2 , 2 β 2 = 1 2 , 3 2 .
These possible ranges for ρ 0 and τ 0 can be extended, if we do not fix γ and δ . We see that, for every ρ > 1 and τ > 1 , there exist γ > 1 and δ > 1 such that ( ρ , γ ) Ω and ( τ , δ ) Ω . Consequently, for every ρ 0 > 1 and τ 0 > 1 , we can choose parameters ρ , τ > 1 and γ , δ > 1 such that the respective collocation-quadrature method (42) looks for approximate solutions of the form (87) with a polynomial p n ( x ) .
Another distinction between [7] and the present investigations is that in [7] the collocation method is studied, the implementation of which is much more expansive (cf. [12,16]) than the collocation-quadrature method considered here.
Of course, the advantage of the results in [7] is that there also the sufficiency of the stability conditions is proved and that in (1) also the case b ( x ) 0 is considered. These problems will be studied for the collocation-quadrature methods considered here in forthcoming papers.
Finally, we can conclude that we were able to prove necessary conditions for the stability of, in comparison with the existing literature, a wider class of collocation-quadrature methods based on the zeros of classical Jacobi polynomials. In this way we can enlarge the range of endpoint singularities of the solutions of singular integral equations of Mellin type, which we can represent in the respective approximate solutions. The questions on the sufficiency of the formulated stability conditions and on the extension of the results presented here to Cauchy singular integral equations remain open for further studies.

Author Contributions

Both authors have contributed to all parts of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Junghanns, P.; Weber, U. Banach algebra techniques for Cauchy singular equations on an interval. In Boundary Element Technology, XII; Southampton: Boston, MA, USA, 1997; Volume 1, pp. 217–256. [Google Scholar]
  2. Junghanns, P.; Kaiser, R. Collocation for Cauchy singular integral equations. Linear Algebra Appl. 2013, 439, 729–770. [Google Scholar] [CrossRef]
  3. Junghanns, P.; Roch, S.; Silbermann, B. Collocation methods for systems of Cauchy singular integral equations on an interval. Vychisl. Tekhnol. 2001, 6, 88–124. [Google Scholar]
  4. Junghanns, P.; Weber, U. Local theory of projection methods for Cauchy singular equations on an interval. In Boundary Integral Methods: Numerical and Mathematical Aspects; WIT Press/Computational Mechanics Publication: Boston, MA, USA, 1999; Volume 1, pp. 217–256. [Google Scholar]
  5. Junghanns, P.; Rathsfeld, A. A polynomial collocation method for Cauchy singular integral equations over the interval. Electron. Trans. Numer. Anal. 2002, 14, 79–126. [Google Scholar]
  6. Junghanns, P.; Rogozhin, A. Collocation methods for Cauchy singular integral equations on the interval. Electron. Trans. Numer. Anal. 2004, 17, 11–75. [Google Scholar]
  7. Junghanns, P.; Rathsfeld, A. On polynomial collocation for Cauchy singular integral equations with fixed singularities. Integral Equ. Oper. Theory 2002, 43, 155–176. [Google Scholar] [CrossRef]
  8. Junghanns, P.; Kaiser, R. A numerical approach for a special crack problem. Dolomites Res. Notes Approx. 2017, 10, 56–67. [Google Scholar]
  9. Junghanns, P.; Kaiser, R. On a collocation-quadrature method for the singular integral equation of the notched half-plane problem. In Large Truncated Toeplitz Matrices, Toeplitz Operators, and Related Topics; Birkhäuser/Springer: Cham, Switzerland, 2017; Volume 259, pp. 413–462. [Google Scholar]
  10. Junghanns, P.; Kaiser, R. A note on Kalandiya’s method for a crack problem. Appl. Numer. Math. 2020, 149, 52–64. [Google Scholar] [CrossRef]
  11. Junghanns, P.; Kaiser, R.; Mastroianni, G. Collocation for singular integral equations with fixed singularities of particular Mellin type. Electron. Trans. Numer. Anal. 2014, 41, 190–248. [Google Scholar]
  12. Capobianco, M.R.; Criscuolo, G.; Junghanns, P. On the numerical solution of a hypersingular integral equation with fixed singularities. In Recent Advances in Operator Theory and Applications; Birkhäuser: Basel, Switzerland, 2009; Volume 187, pp. 95–116. [Google Scholar]
  13. Junghanns, P.; Kaiser, R.; Potts, D. Collocation–quadrature methods and fast summation for Cauchy singular integral equations with fixed singularities. Linear Algebra Appl. 2016, 491, 187–238. [Google Scholar] [CrossRef]
  14. Junghanns, P.; Rost, K. Krylov subspace methods for Cauchy singular integral equations. Facta Univ. Ser. Math. Inform. 2004, 19, 93–108. [Google Scholar]
  15. Junghanns, P.; Rost, K. Matrix representations associated with collocation methods for Cauchy singular integral equations. Math. Methods Appl. Sci. 2007, 30, 1811–1821. [Google Scholar] [CrossRef]
  16. Mastroianni, G.; Frammartino, C.; Rathsfeld, A. On polynomial collocation for second kind integral equations with fixed singularities of Mellin type. Numer. Math. 2003, 94, 333–365. [Google Scholar] [CrossRef]
  17. Elschner, J. On spline approximation for a class of noncompact integral equations. Math. Nachr. 1990, 146, 271–321. [Google Scholar] [CrossRef]
  18. Laurita, C. A new stable numerical method for Mellin integral equations in weighted spaces with uniform norm. Calcolo 2020, 57, 25. [Google Scholar] [CrossRef]
  19. Mastroianni, G.; Monegato, G. Nyström interpolants based on the zeros of Legendre polynomials for a noncompact integral operator equation. IMA J. Numer. Anal. 1994, 14, 81–95. [Google Scholar] [CrossRef]
  20. Monegato, G. A stable Nyström interpolant for some Mellin convolution equations. Orthogonal polynomials and numerical analysis (Luminy, 1994). Numer. Algorithms 1996, 11, 271–283. [Google Scholar] [CrossRef]
  21. Monegato, G.; Prössdorf, S. On the numerical treatment of an integral equation arising from a cruciform crack problem. Math. Methods Appl. Sci. 1990, 12, 489–502. [Google Scholar] [CrossRef]
  22. Junghanns, P.; Kaiser, R. A note on the Fredholm theory of singular integral operators with Cauchy and Mellin kernels. In Operator Theory, Analysis and the State Space Approach; Birkhäuser/Springer: Cham, Switzerland, 2018; Volume 271, pp. 291–325. [Google Scholar]
  23. Junghanns, P.; Kaiser, R. A note on the Fredholm theory of singular integral operators with Cauchy and Mellin kernels, II. In Operator Theory, Functional Analysis and Applications; Birkhäuser/Springer: Cham, Switzerland, 2021; Volume 282, pp. 353–392. [Google Scholar]
  24. Chihara, T.S. An introduction to orthogonal polynomials. In Mathematics and Its Applications; Gordon and Breach, Science Publishers: New York, NY, USA, 1978; Volume 13. [Google Scholar]
  25. Mastroianni, G. Some weighted polynomial inequalities. In Proceedings of the International Conference on Orthogonality, Moment Problems and Continued Fractions, Delft, The Netherlands, 31 October–4 November 1994; Volume 65, pp. 279–292. [Google Scholar]
  26. Junghanns, P.; Mastroianni, G.; Notarangelo, I. Weighted polynomial approximation and numerical methods for integral equations. In Pathways in Mathematics; Birkhäuser/Springer: Cham, Switzerland, 2021. [Google Scholar]
  27. Mastroianni, G.; Totik, V. Uniform spacing of zeros of orthogonal polynomials. Constr. Approx. 2010, 32, 181–192. [Google Scholar] [CrossRef]
  28. Mastroianni, G.; Vértesi, P. Some applications of generalized Jacobi weights. Acta Math. Hungar. 1997, 77, 323–357. [Google Scholar] [CrossRef]
  29. Mastroianni, G.; Russo, M.G. Lagrange interpolation in weighted Besov spaces. Constr. Approx. 1999, 15, 257–289. [Google Scholar] [CrossRef]
  30. Krupnik, N.Y. Banach algebras with symbol and singular integral operators. In Operator Theory: Advances and Applications; Iacob, A., Translator; Birkhäuser Verlag: Basel, Switzerland, 1987; Volume 26. [Google Scholar]
  31. Böttcher, A.; Silbermann, B. Analysis of Toeplitz Operators, 2nd ed.; Prepared jointly with Alexei Karlovich; Springer Monographs in Mathematics; Springer: Berlin, Germany, 2006. [Google Scholar]
  32. Luther, U.; Mastroianni, G. Fourier projections in weighted L-spaces. In Problems and Methods in Mathematical Physics (Chemnitz, 1999); Birkhäuser: Basel, Switzerland, 2001; Volume 121, pp. 327–351. [Google Scholar]
  33. Prössdorf, S.; Silbermann, B. Numerical analysis for integral and related operator equations. In Mathematische Lehrbücher und Monographien, II. Abteilung: Mathematische Monographien [Mathematical Textbooks and Monographs, Part II: Mathematical Monographs]; Akademie-Verlag: Berlin, Germany, 1991; Volume 84. [Google Scholar]
  34. Prössdorf, S.; Silbermann, B. Numerical analysis for integral and related operator equations. In Operator Theory: Advances and Applications; Birkhäuser Verlag: Basel, Switzerland, 1991; Volume 52. [Google Scholar]
  35. Freud, G. Orthogonale Polynome; Lehrbücher und Monographien aus dem Gebiete der Exakten Wissenschaften, Mathematische Reihe, Band 33; Birkhäuser Verlag: Basel, Switzerland, 1969. [Google Scholar]
  36. Muckenhoupt, B. Mean convergence of Jacobi series. Proc. Am. Math. Soc. 1969, 23, 306–310. [Google Scholar] [CrossRef]
  37. Gatteschi, L. On the zeros of Jacobi polynomials and Bessel functions. In Proceedings of the International Conference on Special Functions: Theory and Computation, Turin, Italy, 10–12 October 1984; pp. 149–177. [Google Scholar]
  38. Mastroianni, G.; Milovanović, G.V. Interpolation processes—Basic theory and applications. In Springer Monographs in Mathematics; Springer: Berlin, Germany, 2008. [Google Scholar]
  39. Szego, G. Orthogonal Polynomials, 4th ed.; American Mathematical Society, Colloquium Publications: Providence, RI, USA, 1975; Volume XXIII. [Google Scholar]
Table 1. Cases already considered in the literature.
Table 1. Cases already considered in the literature.
Equation (1)Method
[1] a , b C [ 1 , 1 ] , ρ 0 , τ 0 1 4 , 3 4 , coll. with Jacobi nodes
b ( ± 1 ) = 0
[2,3,4] c ± 0 ρ 0 = τ 0 = 1 2 resp. ρ 0 = τ 0 = 1 2 ,
coll. with Chebyshev nodes
[5,6] c ± 0 ρ 0 , τ 0 1 4 , 3 4 , coll. with Chebyshev nodes
[7] c ± 0 ρ 0 , τ 0 1 4 , 3 4 , coll. with Chebyshev nodes
[8,9,10,11] c ± 0 ρ 0 = τ 0 ± 1 2 ,
coll.-quadr. with Chebyshev nodes
[12,13,14,15,16] c ± 0 numerical aspects, fast algorithms
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Junghanns, P.; Kaiser, R. Quadrature Methods for Singular Integral Equations of Mellin Type Based on the Zeros of Classical Jacobi Polynomials. Axioms 2023, 12, 55. https://doi.org/10.3390/axioms12010055

AMA Style

Junghanns P, Kaiser R. Quadrature Methods for Singular Integral Equations of Mellin Type Based on the Zeros of Classical Jacobi Polynomials. Axioms. 2023; 12(1):55. https://doi.org/10.3390/axioms12010055

Chicago/Turabian Style

Junghanns, Peter, and Robert Kaiser. 2023. "Quadrature Methods for Singular Integral Equations of Mellin Type Based on the Zeros of Classical Jacobi Polynomials" Axioms 12, no. 1: 55. https://doi.org/10.3390/axioms12010055

APA Style

Junghanns, P., & Kaiser, R. (2023). Quadrature Methods for Singular Integral Equations of Mellin Type Based on the Zeros of Classical Jacobi Polynomials. Axioms, 12(1), 55. https://doi.org/10.3390/axioms12010055

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