Varun Thesis
Varun Thesis
net/publication/370340849
CITATIONS READS
0 789
1 author:
Varun ..
University of Lucknow
20 PUBLICATIONS 12 CITATIONS
SEE PROFILE
All content following this page was uploaded by Varun .. on 28 April 2023.
Ph.D. THESIS
SUBMIT T ED T O T HE
UN IVERSIT Y OF LUCKN OW
FOR T HE AWARD OF T HE DEGREE OF
Doctor of Philosophy
IN
MATHEMATICS
BY
VARUN
M.Sc., CSIR(UGC)-NET(JRF)
UN DER
T HE SUPERVISION OF
Dr. SWARNIMA BAHADUR
Self Declaration
I, Varun certify that the work embodied in this Ph.D. thesis entitled:
“A study of different types of interpolating polynomials and their
convergence behaviour on non-uniformly distributed set of nodes
on the unit circle” is my own bona-fide work, which was carried out by
me under the supervision of Dr. Swarnima Bahadur for a period of
January 2017 to October 2022 at the University of Lucknow. The content
embodied in this Ph.D. thesis has not been submitted for the award of any
other degree/diploma. I declare that I have faithfully acknowledged, given
credit to and referred to the researchers wherever their works have been cited
in the text and the body of the thesis. I further certify that I have not wilfully
used any other person’s work, para, text, data, results, etc. reported in the
journals, books, magazines, reports, dissertations, theses, etc., or available
at web-sites/ Internet and included them in this Ph.D. thesis and cited that
as my own work.
Varun
Department of Mathematics
and Astronomy, University of
Lucknow, Lucknow–226 007,
Date:
Uttar Pradesh, India.
Place: Lucknow, India.
Certificate
Dedicated
to my
Father Shri Prakash Kanaujia
Mother Mrs. Prabha
and
Sister Dr. Rimjhim
Preface
The present thesis, which is being submitted for the award of degree of Doc-
tor of Philosophy of the University of Lucknow, Lucknow, is the outcome
of the research work carried by the author since 2017 as a Junior Research
Fellow (CSIR UGC-JRF) and since 2019 as a Senior Research Fellow
(CSIR UGC-SRF) in University of Lucknow.
Organization of Thesis
This thesis is organized in eight chapters as follows:
Chapter 1: Introduction
This chapter presents an overview of some basic definitions, theorems and
historical survey of various types of interpolatory polynomials and their con-
vergence behaviour on some set of nodes, which we will use in this thesis.
Acknowledgements
First and foremost, I am incredibly grateful to my supervisor, Dr. Swarn-
ima Bahadur, for their invaluable advice, continuous support, and patience
during my research work. Her knowledge and experience have encouraged
me in my academic research and daily life. Her insightful feedback pushed
me to sharpen my thinking and brought my work to a higher level. Without
her persistent help, the goal of this work would not have been realized.
I want to express my deep and sincere gratitude towards Prof. Pankaj
Mathur and Prof. Meena Sahai for their valuable suggestions and support
on my research works.
I am grateful to Prof. Bhaskar Srivastava, Prof. Poonam Sharma
and Prof. Vivek Sahai, Head, Department of Mathematics and As-
tronomy, University of Lucknow, Lucknow, for providing the necessary
research facilities during my study. I am thankful to University and Depart-
mental library for providing me the required facilities.
I want to offer my special thanks to my friends Dr. Harendra Verma,
Dhanesh Kumar, Shivam Agarwal, Satish Kumar and Vivek Kumar Yadav
for their unwavering support during all times. Also, I would like to thank my
friends outside of the Mathematics world Er. Mayank Agarwal, Er. Sachin
Nair, Er. Shashank Gaurav and C.A. Apurva Patel for their valuable support
in my good and bad time and always encouraging me to achieve my goal and
move forward in the life.
This acknowledgement would be incomplete without my sincere gratitude
to my parents, Shri Prakash Kanaujia, Mrs. Prabha and my sister
Dr. Rimjhim for their love and understanding, and continuous support to
complete my research work.
University of Lucknow
Lucknow, U.P. Mr. Varun
2021
Page No.
1 Introduction 1
1.1 Best approximation . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Orthogonal Polynomials . . . . . . . . . . . . . . . . . . . . . 4
1.3 Classical Orthogonal Polynomials . . . . . . . . . . . . . . . . 4
1.3.1 Hermite Polynomials . . . . . . . . . . . . . . . . . . . 4
1.3.2 Laguerre Polynomials . . . . . . . . . . . . . . . . . . . 5
1.3.3 Jacobi Polynomials . . . . . . . . . . . . . . . . . . . . 5
1.4 Approximations . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4.1 Types of Approximations . . . . . . . . . . . . . . . . . 6
1.5 Polynomial Interpolation . . . . . . . . . . . . . . . . . . . . . 8
1.6 Nodal systems . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.6.1 Uniform Nodal System . . . . . . . . . . . . . . . . . . 10
1.6.2 Non-Uniform Nodal System . . . . . . . . . . . . . . . 10
1.7 Interpolation on the unit circle . . . . . . . . . . . . . . . . . . 11
1.8 Different Types of Polynomial Interpolation . . . . . . . . . . 12
1.8.1 Lagrange Interpolation . . . . . . . . . . . . . . . . . . 12
1.8.2 Hermite Interpolation . . . . . . . . . . . . . . . . . . . 16
1.8.3 Hermite-Fejér Interpolation . . . . . . . . . . . . . . . 20
1.8.4 Lacunary Interpolation . . . . . . . . . . . . . . . . . . 21
1.8.5 Pál-type Interpolation . . . . . . . . . . . . . . . . . . 27
1.9 What does notation (0,1) and (0,1;2) type Interpolation
means? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.10 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
ix
Interpolation on non-uniform nodes of the unit circle
Introduction
1
Interpolation on non-uniform nodes of the unit circle
where
| x |, if | x |< 1
τ (x) =
1,
if | x |> 1.
then
lim ω(f ; δ) = 0.
δ→∞
Theorem 1.1.1. For every bounded function f(x) defined on the finite inter-
val [a,b] and for any integer n, there exists an ordinary polynomial
Pn (f, x) = a0 + a1 x + a2 x2 + · · · + an xn ,
of degree not greater than n, which deviates least from it on this interval.
The least value of its deviation from f(x) is given by the quantity.
En (f ) =| f (x) − Pn (x) |,
Approximation Theorems
b − a
En (f ) ≤ c ω f ; ,
n
b − a r
(r) b − a
En (f ) ≤ Mr ω f ; ,
n n
If f (x) and g(x) are two real functions with x ∈ [a, b], such that their
Lebesgue-Stieltjes integral is zero, i.e.
Z b
(f, g) = f (x)g(x) dα(x) = 0 (1.1)
a
where f (x), g(x) ∈ Lα (a, b) and dα(x) is the continous or discontinous dis-
tribution of mass in [a, b], then we say that “f (x) is orthogonal to g(x)”.
If we orthogonalize the set of non-negative powers of 1, x, . . . , xn . . . . We ob-
tain a set of polynomials p0 (x), p1 (x), . . . , pn (x), . . . satisfying the conditions:
(i) pn (x) is a polynomial of degree n, in which the coefficients of xn is posi-
tive.
(ii) The sequence {pn (x)} is orthonormal i.e.
Z b
pm (x)pn (x) dα(x) = δmn ; n, m = 0, 1, 2, . . .
a
xy 00 + (α + 1 − x)y 0 + ny = 0.
−1 −α x dn −x n+α
L(α)
n (x) = (n!) x e [e x ].
dxn
dn
Pn(α,β) (x) = (−2)−n (n!)−1 (1 − x)−α (1 + x)−β [(1 − x)n+α (1 + x)n+β ]. (1.3)
dxn
−1
2n
2n (− 1 ,− 1 )
Tn (x) = 2 Pn 2 2 (x).
n
−1
2n + 1
2n (1,1)
Un (x) = 2 Pn 2 2 (x).
n+1
1.4 Approximations
• Linear Approximation
Let f (x) be the function, which we wish to approximate using the class
of functions {gn (x)}, n = 0, 1, 2, . . . , m. Suppose we approximate the
function f (x) by linear combination
• Rational Approximation
Let f (x) be the function, which we wish to approximate using the class
of functions {gn (x)}, n = 0, 1, 2, . . . . We approximate the function
• Interpolation
In the mathematical field of approximation theory and numerical anal-
ysis, one can define interpolation as the technique of obtaining the most
likely estimate of a certain assumptions. In simple words, if a set of
values of a function f (x) at a set of values of x are given, then the
process to determine either the values of x, for which f (x) has some
specified value or the value of f (x) at an intermediate value of set of
values of x is known as interpolation.
For e.g.:- Consider the data in the following table 1.
x 1 2 3 4
The type of interpolation in which data set is interpolated using the polyno-
mial of lowest possible degree that passes through the points of the data set
is known as the polynomial interpolation. The theory of interpolation may
be traced from the well - known theorem of Weierstrass.
| f (x) − Pn (x) |→ 0 as n → ∞.
If the arc length between the consecutive points of the nodal system on the
unit circle remains same, then the nodal system is called uniform nodal sys-
tem.
2kπ
For e.g. zk = exp(iθ) where θ = n
and k = 0, 1, . . . , n − 1.
The nodal system, in which arc length between the consecutive points of the
nodal system on the unit circle can be same or not is a non-uniform nodal
system.
For e.g. zk = exp(iθk ), where k = 1(1)n. It is interesting to note that in
above example if θ1 = θ2 = ... = θn , then the set of non-uniform nodes
transforms to a uniform nodal system.
In the interpolation theory, there are many interesting problems that require
domains other than the real interval. Interpolation problems, in which nodes
of the interpolation lie on the unit circle is known as interpolation on the
unit circle.
Kiŝ [70] was the initiator of interpolation processes on the unit circle.
He considered the (0, 2) and (0, 1, ..., r − 2, r) interpolation for an integer
r ≥ 2 on the nth roots of unity. Darius [43] took the roots of complex
numbers with modulus one as the nodal points and using a suitable modulus
of continuity obtained a result about convergence of the interpolants for
continuous functions.
Recently, Berriochoa et.al [36] studied the Lagrange interpolation on the
unit circle (LIUC) for piecewise continuous functions with suitable properties
and zeros of the para-orthogonal polynomials were considered as the nodal
tion
z11
z21 z22
Q=
z31 z32 z33
.. .. ..
. . .
zn1 zn2 . . . znn
n
X
Ln (f, z) = f (zk )lk (z),
k=1
q(z)
lk (z) = , k = 1, 2, ..., n ,
(z − zkn )q 0 (zkn )
where
q(z) = C(z − z1n )(z − z2n ) . . . (z − znn ).
n
X
lim sup lkn (z) = +∞,
n→∞
k=1
holds densely in a closed set of points that lie on the unit circle |z| = 1.
Erdös and Vertesi [48] in 1980 showed that for a proper function f ∈
[−1, 1], Ln (f, z) will diverge almost everywhere for any set of nodes. In the
same year, A. H. German [56] considered the Lagrange interpolation problem
in the region D̄ = D ∪ T.
Lagrange interpolation on an infinite interval was considered by J. Sz-
abadós and A. K. Varma [114] in 1985 and they also obtained the almost
! !
1
||f − Ln (f )||∞ ≤ cω f, log(n + 1)
n+1
Also, Tian Zi Jiang [65] obtained two infinite asymptotic expansions for
the Lebesgue constant associated with Lagrange polynomial interpolation on
unit circle.In 1990, M. R. Akhlaghi, A. Jakjmovski and A. Sharma [3] con-
sidered Lagrange interpolating polynomial in the unit circle and extended
Walsh’s over convergence theorem to several complex interpolating polyno-
mial. In 1993, G. Mastroianni and P. Vertesi [78] established the mean
convergence of Lagrange interpolation on the general set of nodes.
Also in 1996, R. Brück [38] considered the case of Lagrange interpolation
on non-uniformly distributed nodes on the unit circle. He also studied its
convergence on the nodes, which are the images of the roots of unity, under
a Möbius transform of unit disc onto itself. In 2000, Darius [4] took the
roots of complex numbers with modulus one as the nodal points and using
a suitable modulus of continuity obtained a result about convergence of the
interpolants for continuous functions.
Trefethen [124] breaking the myths about the polynomial interpolation
quoted that polynomial interpolants always converge, if function is smooth
a little bit in Chebyshev points.Thus, function is required to be imposed by
some restrictions such as altering conditions on its modulus of continuity
in order that sequence of Lagrange interpolation polynomials attains finer
properties of the convergence.
In 2012, S. Bahadur [13] studied Lagrange interpolation on the unit circle
by projecting vertically the zeros of (1 − x2 )Pn (x), where Pn (x) denotes the
nth Legendre polynomial and studied the convergence for the same. Recently,
Berriochoa et.al [36] studied the Lagrange interpolation on the unit circle for
piecewise continuous functions with suitable properties and zeros of the para-
orthogonal polynomials were considered as the nodal points in accordance
with analytic weights.
∞
X
f (z) = ak z k ,
k=0
then
lim (Pn−1 (z) − P̆n−1 (z)) = 0, (1.4)
n→∞
Further in 1991, Zi Yu Wang and Shen Ji Tian [133] considered the zeros
0 0
of (1 − x2 )Pn−1 (x), where Pn−1 (x) is the derivative of (n − 1)th Legendre
polynomial and obtained the estimates for the same.
In 1991, R. Sakai [99] considered f ∈ C[−1, 1] and assumed
(α,β)
be the set of zeros of the Jacobi polynomial Pn (x); α, β > −1; with
n+α
Pn(α,β) (1) = and Sn∗ = {Sn } ∪ {±1}
n
Also in 1993, J. Prasad [97] studied the Hermite and Hermite - Fejér
interpolation of higher order on the n distinct zeros of (1 − x2 )Pn−2 (x) where
Pn−2 (x) is the (n − 2)th Legendre polynomial.
In 1993, B. Della Vecchia and G. Mastroianni [131] considered the case
of Hermite interpolation based on the zeros of Jacobi and generalized Jacobi
polynomials in C - space and obtained the estimates in uniform and weighted
Lp −norm.
In 1996, T. N. T. Goodman, K. G. Ivanov and A. Sharma [58] considered
the behaviour of Hermite interpolant on the roots of unity.
In 1997, G. Y. Shi [108] obtained the convergence in Lp - norm of the trun-
cated Hermite interpolants on an arbitrary system of nodes, where truncated
Hermite polynomial is
r X
X n
Lnmr (f, x) = f (j) (xk )Ajk (x),
j=0 k=1
(p)
Ajk (x) = δjp δkq , j, p = 0(1)m − 1 and k, q = 1(1)n.
Q
projecting vertically the zeros of n (x) on the unit circle.
Also, E. Berriochoa, A. Cachafeiro and E. Martinez Breyb [32] studied
the convergence of the Hermite - Fejér and the Hermite interpolation on
polynomials, which are constructed by taking equally spaced nodes on the
unit circle. Consequently, they achieved some improvements on Hermite in-
terpolation problems on the real line. In 2014, E. Berriochoa, A. Cachaferiro
and J. Diaz [35] deals with the order of convergence of the Laurent polyno-
mials of Hermite - Fejér interpolation on the unit circle with nodal system.
The supremum norm of the error of interpolation is obtained for analytic
functions as well as the corresponding asymptotic constants.
In the same year, S. Bahadur and M. Shukla [17] considered Hermite
(α,β)
interpolation on the vertically projected zeros of (1 − x2 )Pn (x), where
(α,β)
Pn (x) stands for Jacobi’s polynomial and proved the following Theorem
“Let f (z) be continuous in |z| ≤ 1 and analytic for |z| < 1 . Let the arbitrary
numbers βk0 s be such that:
2n+1
X 2n+1
X
Hn (z) = f (zk )Ak (z) + βk Bk (z),
k=0 k=0
the Hermite interpolation on the above said nodes.They [18] also considered
the Hermite – interpolation on the vertically projected zeros of the derivative
of the Jacobi polynomial in place of Jacobi’s polynomial. S. Bahadur and
Sariya Bano [25] in 2017 read modified Hermite interpolation problem on
the nodes, which are obtained by projecting vertically the zeros of Legendre
polynomial onto the unit circle.
Recently, E. Berriochoa et al. [37] studied the Hermite interpolation on
the nodal arrays that meet certain separation properties. The novelty of
this article was to develop a interpolation theory that make a practical use
of these nodal systems linked to certain mechanical models that fit these
distribution.
Fejér [49] in 1916 took an important step and proved a theorem, where the
values of the derivatives in the Hermite scheme were equal to zero.
Fejér’s theorem : If f ∈ C[−1, 1], then Hn (f, x) converges to f (x) uni-
formly on [-1,1] as n tends to infinty. Interpolation polynomials Hn (f, x) is
defined by
n
!2
X Tn (x)
Hn (f, x) = f (xkn )(1 − xkn x) ,
k=1
n(x − xkn )
where xkn are the zeros of the Chebyshev polynomial of the first kind.
Hn (f, x) satisfies the below given conditions, where k = 1, 2, ..., n.
Mills [89] in his paper highlights Hermite and Hemite Fejér interpolation
y 00 + f y = 0.
0
πn (x) = (1 − x2 )Pn−1 (x),
where Pn−1 (x) is the (n − 1)th Legendre polynomial. At the outset the exis-
tence and uniqueness was discovered on the abscissas as the zeros of Ultras-
pherical polynomial
1
Pn(λ) (x), λ≥− , n ≥ 4,
2
whereas the explicit representation and convergence theorem were proved for
πn (x) abscissas only.
Also, K. K. Mathur and A. Sharma extended the study of P. Turán on
the zeros of Hermite polynomials, although they confined themselves to the
Theorem 1.8.1. Let f (z) be continuous in |z| ≤ 1 and analytic for |z| < 1.
Further, let ω(δ) be the modulus of continuity of f and
1
lim ω(δ) log =0
δ→0 δ
n
X n
X
Rn (z) = f (zk )Ak (z) + βk Bk (z),
k=1 k=1
Abbasov [1] extended the above result to the two dimensional case. Sharma
[103] has considered the case of (0, 2, 3) - interpolation on the abscissas con-
sidered by Kiŝ. Later on he considered (0, M )- interpolation at equidistant
nodes on unit circle and proved that if a given function satisfies Dini - Lip-
1
lim ω(δ) log =0
δ→0 δ
where ! !
n! n!
βkn =O and vkn = O ,
n3 log n n2 log n
Rn (f, zk,n ) = αk,n ,
0
Rn (f, zk,n ) = βk,n ,
000
Rn (f, zk,n ) = γk,n ,
0
Wn (x) = (1 − x2 )Pn−1 (x)
to be the vertical projections on the unit circle of the zeros of (1 − x2 )Pn (x),
where Pn (x) stands for nth Legendre polynomial having the zeros xk =
cos θk , k = 1(1)n such that
n
Y
Wn (x) = (x − xi,n )
1=1
The roots of yi,n (i = 1(1)n) of Wn0 (x) are interscaled between the roots of
Wn (x) i.e.
−∞ < xn,n < yn−1,n < · · · < x2,n < y1,n < +∞.
He proved that for given arbitrary numbers {a∗i,n }ni=1 and {b∗i,n }n−1
i=1 , there
with initial condition Tn (a) = 0, where a is a given point, different from the
nodal points xi,n (i = 1(1)n, n = 1, 2, . . . ).
In 1985, L. Szilli [120], was the first to apply this method on an infinite in-
terval for even values of n and obtained a uniquely determined polynomial
of degree ≤ 2n − 1. She also established a convergence theorem and showed
that the uniqueness fails to hold for the odd values of n. In 1990, M. R.
Akhlaghi and A. Sharma [2] studied the (0;2) - interpolation problem and
proved regularity and explicit representation of the interpolatory polynomial
when the function values were prescribed at the zeros of πn (x) and the second
derivative at the zeros of πn0 (x). They also considered the analogous prob-
lem by interchanging the nodes. Further in 1991, T. Kilgore and J. Prestin
[69] investigates the Hermite interpolant Hn (f ) of a function f ∈ C 2 [−1, 1],
where the derivatives at ±1 are interpolated upto the order r − 1. They
obtain estimates in the maximum norm over [-1,1] for (f − Hn (f ))(i) in all
possible cases, 0 ≤ i, r − 1 ≤ q. In 1995, K. K. Mathur and R. Srivastava
[83] studied the Pál - type Hermite interpolation on infinite interval. In the
same year, J. Szabados and A. K. Varma [117] studied the Pál - type (0,2)
Rx
interpolation on the mixed zeros of Pn−1 (x) and −1
Pn−1 (t)dt and proved
that the interpolatory polynomial approximates the continuous function f
on [-1,1] in Telyakovskii - Gopengauz order. Also, T. F. Xie [141] considered
(0; 1; 2) - interpolation polynomials for approximation.
Further, in 1996, I. Joó and L. G. Pál [68] obtained Lacunary (0;0,1) - inter-
polation on the roots of Jacobi polynomials and its derivatives respectively.
Also, I. Joó and V. E. S. Szabó [66] gave a generalization of Pál - type in-
terpolation process. Also, L. G. Pál [93] gave a general Lacunary (0;0,1)-
interpolation process. Later on, R. Srivastava and K. K. Mathur [84] consid-
ered the weighted (0;0,2) - interpolation on the roots of Hermite polynomials.
In 2003, H. P. Dikshit [46] also considered the Pál - type interpolation on non
- uniformly distributed nodes on the unit circle involving the Möbius trans-
forms of the zeros of z n +1 and z n −1 along with one additional point or two
additional points and 1. Marcel G.de Brüin and A. Sharma [39] considered
Lacunary Pál - type interpolation, where they considered (0, m1 , m2 , . . . , mq )
- interpolation problem on the zeros of
with
0 < α0 < · · · < αq
Bahadur [10] considered Pál – type (0,1’0) - interpolation on two disjoint set
of nodes obtained by projecting vertically the zeros of Legendre polynomial
(1 − x2 )Pn (x) and Pn0 (x) respectively on the unit circle. In another paper,
she [8] also considered (0;0,1) - interpolation by interchanging the nodes. In
2013, M. Lénard [76] considered two set of nodes (one consists of the zeros
of a polynomial pn of degree n, while the elements of the other one are the
zeros of p0n ) different interpolation conditions are prescribed simultaneously.
Weighted (0, 1, . . . , r − 2, r) - interpolation conditions (r ≥ 2) are given on
one of the sets of the nodes with Hermite conditions on the other set and she
studied the regularity of the problem and obtained the explicit representation
of the fundamental polynomials of interpolation. Also, V. Srivastava et. al.
[109] considered trigonometric polynomials
cos (2n+1)θ
2
sin (2n+1)θ
2
rn (x) = θ
and sn (x) =
cos 2 sin 2θ
such that their zeros are interscaled and neither of the polynomials is the
extrema of the other polynomial. On these interscaled set of nodes, they
studied an interpolation process, when Lagrange data is prescribed on one
set of nodes and Hermite data is prescribed on the other set of nodes. In an-
other similar paper they [110] considered the converse of the above problem
by interchanging the conditions on the function and its derivative. Later on
S. Bahadur {[14],[15]} considered (0;0,2) and (0,2;0) - interpolation on the
vertically projected zeros of (1 − x2 )Pn (x) and Pn0 (x) onto the unit circle. Af-
ter that, she (with M. Shukla) {[19],[20]} considered the Hermite - Lagrange
and Lagrange - Hermite interpolation on the vertically projected zeros of
(α,β) (α,β)0
(1 − x2 )Pn (x) and Pn (x) onto the unit circle, obtained the regularity,
explicit forms and established a convergence theorem for that interpolatory
polynomial.
Interpolation means?
Sometimes readers will come across the notation (0, 1) and (0, 1; 2) type in-
terpolation. The comma in the notation tells us that the function and their
derivatives have been prescribed on a particular set of nodes, whereas semi-
colon in the notation tells us that the different nodal systems have been
considered. For e.g.
(0,1) interpolation: The absence of semicolon tells that only one set of
nodes is considered, whereas comma tells us that function and its first deriva-
tive are prescribed on that same set of nodes.
(0,1;2) interpolation: Single semicolon tells us that two different nodal sys-
tems have been considered for the interpolation. Also, the first half before
the semicolon tells that conditions on the function and its first derivatives
are prescribed on the one set of nodal system, whereas the second half after
the semicolon tells that conditions on the second derivative of the function
are prescribed on the other set of nodal system.
1.10 Preliminaries
(α,β)
Throughout the thesis, Jacobi polynomial of degree n is denoted by Pn (x),
where α, β > −1.
Let us define two sets T = {z ∈ C : |z| = 1} and D = {z ∈ C : |z| < 1}.
(α,β)
The differential equation satisfied by Pn (x) is
00 0
(1 − x2 )Pn(α,β) (x) + [β − α − (α + β + 2)x]Pn(α,β) (x)
1 + z2
Using the Szegő transformation x = , we have
2z
00
(z 2 − 1)4 Pn(α,β) (x) + 4z(z 2 − 1) {(α + β + 2)z 2 + 1}(z 2 − 1)
0
−2z 3 (β − α) Pn(α,β) (x) − 16z 6 n(n + α + β + 1)Pn(α,β) (x) = 0.
(1.9)
Let Z2n and T2n−2 be two distinct sets of nodes, which are obtained by
(α,β) (α,β)0
projecting vertically the zeros of Pn (x) and Pn (x) respectively on the
unit circle.
k = 1, 2, ..., n ; xk , yk ∈ R} (1.10)
(α,β)
where, xk = cos θk , k = 1(1)n are the zeros of Pn (x) and x∗k = cos φk ,
(α,β)0
k = 1(1)n − 1 are the zeros of Pn (x) .
The nodal polynomials W (z) and W ∗ (z) defined on Z2n and T2n−2 are given
by (1.12) and (1.13) respectively.
2n
!
Y 1 + z2 n
W (z) = (z − zk ) = Kn Pn(α,β) z , (1.12)
k=1
2z
where
Γ(α + β + n + 1)
Kn = 22n n! .
Γ(α + β + 2n + 1)
and !
2n−2
Y 0 1 + z 2 n−1
W (z) =
∗
(z − tk ) = Kn∗ Pn(α,β) z , (1.13)
k=1
2z
where
Γ(α + β + n + 1)
Kn∗ = 22n−1 (n − 1)! .
Γ(α + β + 2n + 1)
Let
R(z) = (z 2 − 1)W (z). (1.14)
W (z)
Lk (z) = , k = 1, 2, ..., 2n. (1.15)
(z − zk )W 0 (zk )
W ∗ (z)
Lk∗ (z) = , k = 1, 2, ..., (2n − 2). (1.16)
(z − tk )W ∗0 (tk )
R(z)
Lk (z) = , k = 0, 1, 2, ..., 2n + 1, (1.17)
(z − zk )R 0 (zk )
√
| z 2 − 1 |= 2 1 − x2 . (1.18)
and r
√ √ q
| z − zk |= 2 1 − xxk − 1− x2 1 − x2k . (1.19)
0
(1 − x2 ) | Pn(α,β) (x) |= O(nα+1 ). (1.22)
α 1 1
(1 − x2 ) 2 + 4 | Pn(α,β) (x) |= O(n− 2 ). (1.23)
0
| Pn(α,β) (x) |= O(nα+2 ). (1.24)
00
| Pn(α,β) (x) |= O(nα+4 ). (1.25)
Considering set of nodes Z2n and T2n−2 such that for each k , xk , x∗k ∈ (−1, 1),
we have !−2
k
(1 − x2k )−1 ∼ . (1.26)
n
1
| Pn(α,β) (x∗k ) |∼ k −α− 2 nα . (1.27)
0 3
| Pn(α,β) (xk ) |∼ k −α− 2 nα+2 . (1.28)
00 5
| Pn(α,β) (xk ) |∼ k −α− 2 nα+4 . (1.29)
Convergence of
Lagrange-Hermite
Interpolation on unit circle
LHn (f, zk ) = f (zk ),
k = 0(1)2n + 1,
(2.1)
LH0n (f, ±1) = α±1 ,
36
Interpolation on non-uniform nodes of the unit circle
Let us consider
LHn (z) = R(z) q(z), (2.2)
Polynomials
We shall write
2n+1
X X
LHn (z) = f (zk )Ak (z) + α±1 Bk (z), (2.5)
k=0 k=0,2n+1
where Ak (z) and Bk (z) are fundamental polynomials, each of degree at-
most 2n + 3 satisfying the conditions (2.6) and (2.7) respectively.
Ak (zj ) = δkj
; j, k = 0(1)2n + 1,
(2.6)
A0k (zj ) = 0
; k = 0(1)2n + 1, j = 0, 2n + 1.
Bk (zj ) = 0
; j = 0(1)2n + 1, k = 0, 2n + 1
(2.7)
Bk0 (zj ) = δkj
; j, k = 0, 2n + 1.
R(z)(z + zk )
Bk (z) = , (2.8)
4Kn
where k = 0, 2n + 1.
where t(z) is a linear polynomial. From the first condition in (2.7), one can
verify that Bk (zj ) = 0 for j = 0(1)2n + 1.
Differentiating (2.9) with respect to z , we get
Bk0 (z) = 2zW (z)t(z) + (z 2 − 1)W 0 (z)t(z) + (z 2 − 1)W (z)t0 (z). (2.10)
For j = 0, 2n + 1, we have
Using the second set of condition of (2.7), (1.15) and (1.17), we get t(1) =
1 −1
and t(−1) = .
2Kn 2Kn
Therefore, we have
z + zk
t(z) = . (2.12)
4Kn
and for k = 0, 2n + 1
h 1 1 L0k (zk ) i
Ak (z) = (z + zk )Lk (z) − + (z − zk ) . (2.14)
2zk 4zk2 2zk
At z = zj , we get
To satisfy the second set of conditions in (2.6), A0k (zj ) = 0. Thus, we get
L0k (zj )
ck = − . (2.17)
(zj + zk )R 0 (zj )
0
Lk (z)(z − zk )R (zk ) = R(z). (2.18)
0 0
L0k (z)(z − zk )R (zk ) + Lk (z)R (zk ) = R 0 (z). (2.19)
0
L0k (zj )(zj − zk )R (zk ) = R 0 (zj ). (2.20)
1
ck = . (2.21)
R 0 (zk )(zk2 − zj2 )
1
ck = . (2.22)
R 0 (zk )(zk2 − 1)
1
t0 (z0 ) =
2 (2.24)
1
t2n+1 (z2n+1 ) = −
2
Similarly, satisfying the second set of conditions in (2.6), we get
11
t00 (z0 ) =− 0
+ L0 (z0 )
2 2 (2.25)
11
t02n+1 (z2n+1 ) = − L02n+1 (z2n+1 )
2 2
1
tk (zk ) =
2zk
(2.26)
0 tk (zk ) 0
tk (zk ) = − 1 + 2zk Lk (zk )
2zk
tk (z) = a1 + b1 (z − zk ), (2.27)
1
tk (zk ) = a1 =
2zk
(2.28)
1 L0k (zk )
t0k (zk ) = b1 = − +
4zk2
2zk
h 1 1 L0k (zk ) i
tk (z) = − + (z − zk .
)
2zk 4zk2 2zk
We get the desired result (2.14) by putting the value of tk (z) in (2.23). Hence
the theorem follows.
2n+1
X
max | Lk (z) | ≤ c log n, (2.29)
|z|=1
k=0
r
2
| Bk (z) | ≤ . (2.30)
πn
R(z)(z + zk )
| Bk (z) |= .
4Kn
(0,0)
|z 2 − 1||Pn (x)||z n ||z + zk |
| Bk (z) |= .
4
√
| Bk (z) |≤ 1 − x2 |Pn(0,0) (x)|.
2n+1
X
| Ak (z) |≤ c1 n1/2 log n,
k=0
(z + zk )R(z)
Ak (z) = Lk (z) + .
(zk2 − 1)R 0 (zk )
Then,
|z + zk ||R(z)|
|Ak (z)| ≤ |Lk (z)| + .
|zk2 − 1||R 0 (zk )|
(0,0)
|z + zk ||z 2 − 1||Pn (x)||z n |
|Ak (z)| ≤ |Lk (z)| + (0,0)0
.
|zk2 − 1|2 |Pn (xk )||zkn |
√ (0,0)
1 − x2 |Pn (x)|
|Ak (z)| ≤ |Lk (z)| + (0,0)0
.
(1 − x2k )|Pn (xk )|
n1/2
|Ak (z)| ≤ |Lk (z)| + .
k
2n 2n 2n
X X X n1/2
|Ak (z)| ≤ |Lk (z)| + ,
k=1 k=1 k=1
k
X2n X2n
|Ak (z)| ≤ |Lk (z)| + n1/2 log n. (2.31)
k=1 k=1
" ! #
1 1 L0k (zk )
Ak (z) = (z + zk )Lk (z) − + (z − zk ) .
2zk 4zk2 2zk
X X
|Ak (z)| ≤ |Lk (z)| + c2 + c3 n1/2 , (2.32)
k=0,2n+1 k=0,2n+1
2n+1
X 2n+1
X
|Ak (z)| ≤ |Lk (z)| + c2 + c3 n1/2 + n1/2 log n. (2.33)
k=0 k=0
2n+1
X
|Ak (z)| ≤ c log n + c2 + c3 n1/2 + n1/2 log n. (2.34)
k=0
2.4 Convergence
2n+1
X X
LHn (z) = f (zk )Ak (z) + α±1 Bk (z), (2.36)
k=0 k=0,2n+1
| LHn (z) − f (z) | = O(ω2 (f, n−1 )n1/2 log n), (2.37)
be expressed as
2n+1
0
X X
Fn (z) = Fn (zk )Ak (z) + Fn (zk )Bk (z). (2.41)
k=0 k=0,2n+1
2n+1
X
| LHn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Ak (z) |
k=0
X
+ | α±1 − Fn0 (zk ) || Bk (z) | + | Fn (z) − f (z) |,
k=0,2n+1
2n+1
X X
| LHn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Ak (z) | + | α±1 || Bk (z) |
k=0 k=0,2n+1
| {z } | {z }
N1 N2
X
+ | Fn0 (zk ) || Bk (z) | + | Fn (z) − f (z) | . (2.42)
| {z }
k=0,2n+1 N4
| {z }
N3
| LHn (z) − f (z) | = O(ω2 (f, n−1 )n1/2 log n). (2.47)
Convergence of Interpolatory
polynomial between Lagrange
and Hermite
HLn (f, zk ) = τk
; k = 0, 2n + 1,
(3.1)
HL0n (f, zk ) = βk
; k = 1(1)2n,
49
Interpolation on non-uniform nodes of the unit circle
Z z
HLn (z) = a1 W (t)dt. (3.4)
−1
Z 1
HLn (1) = a1 W (t)dt. (3.5)
−1
R1
Since, −1
W (t)dt 6= 0. This implies a1 = 0.
Hence, the theorem follows.
Polynomials
We shall write
X 2n
X
HLn (z) = τk Uk (z) + βk Vk (z), (3.6)
k=0,2n+1 k=1
where Uk (z) and Vk (z) are unique polynomials, each of degree atmost
Uk (zj ) = δkj
; j, k = 0, 2n + 1
(3.7)
Uk0 (zj ) = 0
; k = 0, 2n + 1, j = 1(1)2n,
and
Vk (zj ) = 0
; j = 0, 2n + 1, k = 1(1)2n
(3.8)
Vk0 (zj ) = δkj
; j, k = 1(1)2n.
Z z
2
Vk (z) = (z − 1)Lk (z) − (t2 − 1)Lk0 (t)dt
−1
Z z Z z
Sk (1)
+ (1 − 2zk ) Lk (t)dt − R 1 W (t)dt. (3.9)
−1 −1
W (t)dt −1
Proof. Let
Z z
2
Vk (z) = (z − 1)Lk (z) + Sk (z) + a W (t)dt, (3.10)
−1
where Vk (z) is atmost of degree 2n+1 satisfying the conditions given in (3.8),
then
Z 1
Vk (1) = Sk (1) + a W (t)dt,
−1
Z 1
0 = Sk (1) + a W (t)dt,
−1
Sk (1)
a = −R 1 . (3.11)
−1
W (t)dt
Vk0 (z) = 2zLk (z) + (z 2 − 1)Lk0 (z) + Sk0 (z) + aW (z). (3.12)
At z = zj , we have
Vk0 (zj ) = 2zj Lk (zj ) + (zj2 − 1)Lk0 (zj ) + Sk0 (zj ). (3.13)
Z z Z z
Sk (z) = − 2
(t − 1)Lk0 (t)dt + ak Lk (t)dt. (3.15)
−1 −1
Now, to determine the constant ak , we use the second set of conditions given
in (3.8) for j = k,
Vk0 (zk ) = 2zk Lk (zk ) + (zk2 − 1)Lk0 (zk ) + Sk0 (zk ). (3.16)
ak = 1 − 2zk . (3.17)
Hence, from (3.10), (3.11), (3.15) and (3.17), we have theorem 3.2.1 .
Proof. Let
Z z
Uk (z) = (−1) bkk
W (t)dt, (3.19)
−zk
Z zj
k
Uk (zj ) = (−1) bk W (t)dt. (3.20)
−zk
1
bk = R 1 .
−1
W (t)dt
2n
X
| Vk (z) | = O(log n).
k=1
Z z
2
Vk (z) = (z − 1)Lk (z) − (t2 − 1)Lk0 (t)dt
−1
Z z Z z
Sk (1)
+ (1 − 2zk ) Lk (t)dt − R 1 W (t)dt. (3.21)
−1 −1
W (t)dt −1
Z z Z z
2
|Vk (z)| ≤ |(z − 1)Lk (z)| + (1 − t 2
)Lk0 (t)dt + (1 − 2zk ) Lk (t)dt
| {z } −1 −1
I1 | {z }
I2
Z z
Sk (1)
+ R1 W (t)dt . (3.22)
W (t)dt −1
| −1 {z }
I3
2n
X 2n
X 2n
X 2n
X
|Vk (z)| ≤ I1 + I2 + I3 . (3.23)
k=1 k=1 k=1 k=1
(0,0)
|(z 2 − 1)||Pn (x)||z n |
I1 = (0,0)0
. (3.25)
|(z − zk )||Pn (xk )|
√ (0,0)
2 1 − x2 |Pn (x)|
I1 = √ q √ ,
(0,0)0
p
2 2
2 1 − xxk − 1 − x 1 − xk |Pn (xk )|
√ √ (0,0)
q √ p
2 1 − x |Pn (x)| 1 − xxk + 1 − x2 1 − x2k
2
= p (0,0)0
,
(1 − xxk )2 − (1 − x2 )(1 − x2k ) |Pn (xk )|
√ √ (0,0)
q p
2 1 − x |Pn (x)| 1 − xxk + (1 − xxk )2 − (x − xk )2
2
= (0,0)0
,
|x − xk | |Pn (xk )|
√ √ (0,0) √
2 1 − x2 |Pn (x)| 1 − xxk
≤ (0,0)0
.
|x − xk | |Pn (xk )|
1
For | x − xk | ≥ | 1 − x2k |, we have
2
√ √ (0,0) √
2 2 1 − x2 |Pn (x)| 1 − xxk
I1 ≤ (0,0)0
, (3.26)
(1 − x2k )|Pn (xk )|
√ (0,0)
4 1 − x2 |Pn (x)|
≤ (0,0)0
.
(1 − x2k )|Pn (xk )|
!
n−1
I1 = O .
(k 2 n−2 )(k −3/2 n2 )
1
I1 = O(n−1 k − 2 ). (3.27)
1
The estimate remains the same in the case, when | x − x∗k | < | 1 − x∗2
k |.
2
On taking summation , we get
2n
!
X log n
I1 = O 1 . (3.28)
k=1 n2
Now, to determine I2 , let us first determine Lk0 (z). From (1.15), we have
W (z)
Lk (z) = .
(z − zk )W 0 (zk )
0
Lk (z)(z − zk )W (zk ) = W (z).
0 0
Lk0 (z)(z − zk )W (zk ) + Lk (z)W (zk ) = W 0 (z).
0
W 0 (z) − Lk (z)W (zk )
Lk0 (z) = .
(z − zk )W 0 (zk )
W 0 (z) Lk (z)
|Lk0 (z)| ≤ 0 + . (3.29)
(z − zk )W (zk ) (z − zk )
Z z Z z
I2 = 2
(1 − t )Lk0 (t)dt + (1 − 2zk ) Lk (t)dt .
−1 −1
Z z Z z
I2 ≤ 2
| (t − 1) || Lk0 (t) | dt + (1 + 2 | zk |) | Lk (t) | dt.
−1 −1
Using (3.29), (1.12), (1.18), (1.19), (1.22), (1.23), (1.26), (1.28) and following
the similar scheme as used to evaluate I1 , we get
2n
X
I2 = O(log n). (3.30)
k=1
Also Rz
W (t)dt
I3 = Sk (1) R−1
1 .
−1
W (t)dt
2n
X
I3 = O(1). (3.31)
k=1
X
| Uk (z) |= O(1).
k=0,2n+1
Proof.
X
| Uk (z) |=| U0 (z) | + | U2n+1 (z) | .
k=0,2n+1
Z z
1
| U0 (z) |≤ R 1 | W (t)dt | .
| −1 W (t)dt | −1
Similarly, we have
| U2n+1 (z) |= O(1). (3.33)
3.4 Convergence
X 2n
X
HLn (z) = f (zk )Uk (z) + βk Vk (z), (3.35)
k=0,2n+1 k=1
X
| HLn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Uk (z) |
k=0,2n+1
2n
X
+ | βk − Fn0 (zk ) || Vk (z) | + | Fn (z) − f (z) |,
k=1
X 2n
X
| HLn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Uk (z) | + | βk || Vk (z) |
k=0,2n+1 k=1
| {z } | {z }
M1 M2
2n
X
+ | Fn0 (zk ) || Vk (z) | + | Fn (z) − f (z) | . (3.41)
| {z }
k=1 M4
| {z }
M3
A note on Hermite-Fejér
Interpolation on the
non-uniform nodes of the unit
circle
This chapter deals with the Hermite-Fejér interpolation problem on the unit
circle with the nodal system containing the vertically projected zeros of Ja-
cobi’s polynomial with boundary points on the unit circle. Three cases with
different nodal structures is worked upon in this chapter and rate of conver-
gence is obtained for each case. Moreover, a comparison of all the three cases
and provided some important conclusions.
61
Interpolation on non-uniform nodes of the unit circle
Polynomial
Referring to the Case I, II and III mentioned in (4.1) - (4.3), we give following
theorems.
2n
X
Jn (z) = αk Mk (z), (4.4)
k=1
n o
Mk (z) = Lk2 (z) 1 − 2(z − zk )Lk0 (zk ) , (4.5)
Proof. Let
Mk (z) = mk W (z)Lk (z) + Lk2 (z), (4.7)
At z = zk , we have
M0k (zk ) = mk W 0 (zk )Lk (zk ) + mk W (zk )Lk0 (zk ) + 2Lk (zk )Lk0 (zk ).
2Lk0 (zk )
mk = − . (4.8)
W 0 (zk )
2n+1
X
Gn (z) = βk Nk (z), (4.9)
k=0
!
Lk0 (zk ) + L0k (zk )
Nk (z) = − W (z)Lk (z) + Lk (z)Lk (z). (4.10)
W 0 (z )
k
For k = 0, 2n + 1
W (z)Lk (z)
Nk (z) = , (4.11)
W (zk )
Nk0 (z) = nk W 0 (z)Lk (z) + nk W (z)L0k (z) + Lk0 (z)Lk (z) + Lk (z)L0k (z).
At z = zk , we have
Nk0 (zk ) = nk W 0 (zk )Lk (zk ) + nk W (zk )L0k (zk ) + Lk0 (zk )Lk (zk ) + Lk (zk )L0k (zk ).
L 0 (z ) + L0 (z )
k k k k
nk = − . (4.14)
W (zk )
0
For k = 0, 2n + 1’
let
Nk (z) = n1k W (z)Lk (z), (4.15)
1
n1k = , as W (zk ) 6= 0. (4.16)
W (zk )
2n+1
X
Hn (z) = γk Pk (z), (4.17)
k=0
n o
Pk (z) = L2k (z) 1 − 2(z − zk )L0k (zk ) , (4.18)
2n
X 2n
X
| Mk (z) | ≤ max | Lk (z) | | Lk (z) |
k=1 k=1
2n
X | L 0 (zk ) || Lk (z) || W (z) |
k
+ 2 . (4.21)
k=1
| W 0 (zk ) |
2n 2n 2n
X X 1 X 1
| Mk (z) |≤ C1 −α+ 12
+ C2 , (4.22)
k=1 k=1 k k=1
k −2α+1
2n+1
X 2n
X X
Nk (z) = Nk (z) + Nk (z).
k=0 k=1 k=0,2n+1
2n+1
X 2n
X X
| Nk (z)| ≤ | Nk (z)| + | Nk (z) | . (4.24)
k=0 k=1 k=0,2n+1
| {z } | {z }
I1 I2
Now,
2n 2n
X | L 0 (zk ) || Lk (z) || W (z) |
k
X | L0 (zk ) || Lk (z) || W (z) |
k
I1 ≤ +
k=1
|W 0 (z )
k | k=1
| W 0 (zk ) |
X2n
+ | Lk (z) || Lk (z) | . (4.25)
k=1
Now, we have
X | Lk (z) || W (z) |
I2 ≤ . (4.27)
k=0,2n+1
| W (zk ) |
Combining (4.26) and (4.28), after a little computations for different ranges
of α, the lemma follows.
4.4 Convergence
3
ωr (f, n−1 ) = O(n−r+1−q ), {q > + α − r}, (4.30)
2
where ωr (f, n−1 ) be the rth modulus of continuity of f (z), such that
h 3 1i
r− ≤α<r− , (4.31)
2 2
where r = 2, 3, ... .
To prove the theorems , we shall need following.
Let f (z) be a function continuous on T ∪ D and analytic on D. Then, there
exists a polynomial Fn (z) of degree ≤ 4n − 1 satisfying Jackson’s inequality
[62].
1
O(ω(f, n−1 )),
−1 < α ≤ − ,
2
1
| f (z) − Fn (z) |= O(ω2 (f, n−1 )), − < α ≤ ,1 (4.32)
2 2
1
O(ωr (f, n−1 )), α > ,
2
where C is a constant independent of n and z.
2n
X
Jn (z) = f (zk )Mk (z), (4.33)
k=1
2n
X
Fn (z) = Fn (zk )Mk (z), (4.35)
k=1
Then,
2n
X
| Jn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Mk (z) | + | Fn (z) − f (z) | . (4.37)
k=1
2n+1
X
Gn (z) = f (zk )Nk (z), (4.38)
k=0
2n+1
X
Hn (z) = f (zk )Pk (z), (4.40)
k=0
4.5 Conclusion
Extension of Pál-type
Hermite-Fejér Interpolation
onto the unit circle
of Interpolatory Polynomial
73
Interpolation on non-uniform nodes of the unit circle
Rn (zk ) = αk ; k = 1, 2, ..., 2n,
R0n (tk ) = βk ; k = 1, 2, ..., (2n − 2), (5.1)
R(m)
n (±1) = 0 ; m = 0, 1, ..., r ,
2n
X 2n−2
X
Rn (z) = αk Ak (z) + βk Bk (z). (5.2)
k=1 k=1
Here, Ak (z) and Bk (z) are the first and second fundamental polynomials each
of degree ≤ 4n + 2r − 1 satisfying (5.3) and (5.4) respectively.
For k = 1, 2, ..., 2n
Ak (zj ) = δkj ; j = 1, 2, ..., 2n,
0
Ak (tj ) = 0 ; j = 1, 2, ..., 2n − 2, (5.3)
A(m)
k (±1) = 0 ; m = 0, 1, ..., r ,
Z z
Jk (z) = z n+1 (z 2 − 1)r Lk∗ (z) dz (5.5)
0
and
Z z
J1j (z) = z n−j (z 2 − 1)r W ∗ (z) dz ; j = 0, 1 , (5.6)
0
Bk (z) = z −n
W (z) bk Jk (z) + b0k J10 (z) + b1k J11 (z) , (5.7)
where
1
bk = , (5.8)
W (tk )tk (t2k − 1)r
and
bk (Jk (−1) − Jk (1))
b0k = . (5.10)
2J10 (1)
At z = zj , j = 1, 2, ..., n
Bk (zj ) = zj−n W (zj ) bk Jk (zj ) + b0k J10 (zj ) + b1k J11 (zj ) .
Bk0 (z) = − nz −n−1 W (z) + z −n W 0 (z) bk Jk (z) + b0k J10 (z) + b1k J11 (z)
Since tj ’s are the zeros of the polynomial W ∗ (z), so Bk0 (tj ) = t−n
j W (tj )bk Jk (tj ).
0
j W (tj )bk tj
Bk0 (tj ) = t−n (tj − 1)r Lk∗ (tj )
n+1 2
0
Using condition Bk (tj ) = δkj given in (5.4), at j = k, we get (5.8).
One can verify the results for j 6= k.
(m)
Also, from Bk (±1) = 0 for m = 0, 1, ..., r, we get (5.9) and (5.10).
Hence, the theorem 5.1.1 follows.
by (5.11)
0
z
z n (z 2 − 1)r (z − 1)W ∗ (z) + ck W ∗ (z)
Z 2
Sk (z) = − dz,
0 W 0 (zk )(zk2 − 1)r+1 W ∗ (zk ) (z − zk )
(5.12)
Sk (1) − Sk (−1)
a0k = , (5.13)
2J10 (1)
and
0
(1 − zk 2 )W ∗ (zk )
ck = . (5.15)
W ∗ (zk )
1 h
A0k (z) = W ∗
2z(r + 1)(z 2 − 1)r Lk (z) +
(z)
(zk2 − 1)r+1 W ∗ (zk )
0
i
(z 2 − 1)r+1 Lk0 (z) + (z 2 − 1)r+1 Lk (z)W ∗ (z)
n o
+ − nz −n−1 W (z) + z −n W 0 (z) Sk (z) + a0k J10 (z) + a1k J11 (z)
+z W (z) Sk (z) + a0k J10 (z) + a1k J11 (z) .
−n 0 0 0
0
(t2j − 1)r+1 Lk (tj )W ∗ (tj ) n o
A0k (tj ) = + − nt−n−1
W (t j ) + t−n
W 0
(t j )
(zk2 − 1)r+1 W ∗ (zk ) j j
j W (tj )Sk (tj ).
Sk (tj ) + a0k J10 (tj ) + a1k J11 (tj ) + t−n 0
0
(t2j − 1)r+1 Lk (tj )W ∗ (tj ) n o
0 = + − nt−n−1
W (t j ) + t−n
W 0
(t j )
(zk2 − 1)r+1 W ∗ (zk ) j j
j W (tj )Sk (tj ),
Sk (tj ) + a0k J10 (tj ) + a1k J11 (tj ) + t−n 0
0
(t2j − 1)r+1 Lk (tj )W1∗ (tj ) n −n−1 o
j W
t−n W W
−n
(tj )Sk0 (tj ) = − + nt (tj ) − t 0
(t j )
(zk2 − 1)r+1 W ∗ (zk ) j j
× Sk (tj ) + a0k J10 (tj ) + a1k J11 (tj )
0
( !
(t2j − 1)r+1 Lk (tj )W ∗ (tj ) 1 + t2
j
= − + nt−n−1 Kn Pn(α,β) tnj
(zk2 − 1)r+1 W ∗ (zk ) j
2tj
( ! ! ))
2 2
0 1 + tj 1 + tj
− t−n
j Kn Pn(α,β) tnj + nKn Pn(α,β) tn−1
j
2tj 2tj
!
× Sk (tj ) + a0k J10 (tj ) + a1k J11 (tj ) .
0
(t2j − 1)r+1 Lk (tj )W ∗ (tj )
Sk0 (tj ) = − −n .
tj W (tj )(zk2 − 1)r+1 W ∗ (zk )
0
z n (z 2 − 1)r (z − 1)W ∗ (z) + ck W ∗ (z)
2
Sk0 (z) =− 0 . (5.16)
W (zk )(zk2 − 1)r+1 W ∗ (zk ) (z − zk )
0
[(z 2 − 1)W ∗ (z) + ck W ∗ (z)]|z=zk = 0,
which in turn gives (5.15). Similarly, the constants a0k and a1k can be found
(m)
out by satisfying the condition Ak (±1) = 0 for m = 0, 1, ..., r.
Hence, the theorem 5.1.2 follows.
2n
X
2 r/2 r+1
| Ak (z) | = O (1 − x ) n log n , (5.17)
k=1
r
where −1 < α ≤ .
2
2n 2n 2n
X X (z 2 − 1)r+1 Lk (z)W ∗ (z) X
| Ak (z) | ≤ 2 r+1 W ∗ (z )
+ z −n W (z)Sk (z)
k=1 k=1
(zk − 1) k
k=1
| {z } | {z }
I1 I2
2n
X
+ z −n W (z)(a0k J10 (z) + a1k J11 (z)) .
k=1
| {z }
I3
We can write as
2n
X
| Ak (z) | ≤ I1 + I2 + I3 (5.18)
k=1
2n
X (z 2 − 1)r+1 W (z)W ∗ (z)
I1 =
k=1
(zk2 − 1)r+1 (z − zk )W 0 (zk )W ∗ (zk )
n o
2 r+1 (α,β) 1+z 2
X2n (z − 1) K n Pn 2z
zn
= n o0 ×
2 (α,β) 1+z 2
k=1 (zk − 1)r+1 (z − zk ) Kn Pn 2z
z n
|z=zk
n 0
2
o
(α,β) 1+z
Kn∗ Pn 2z
z n−1
n o
(α,β)0 1+zk2 n−1
Kn∗ Pn 2zk
z k
(α,β)
Since Pn (xk ) = 0 and | zk |= 1, we get
2n (α,β) (α,β)0
X | (z 2 − 1)r+1 || Pn (x) || Pn (x) |
I1 = 2 (α,β)0
.
k=1 | (zk2 − 1)r+2 || z − zk || Pn (xk ) |2
1
For | x − xk | ≥ 2
| 1 − x2k |, we have
2n
!
X 1
I1 = O (1 − x2 )r/2 nr+1 .
k=1
k r−2α+1
From r − 2α + 1 ≥ 1, we get
n ro
2 r/2 r+1
I1 = O (1 − x ) n log n −1<α≤ . (5.19)
2
The reader can verify that estimate remains the same in the case, where
1
| x − xk | < 2
| 1 − x2k |.
Following similar scheme as above gives
n ro
2 r/2 r
I2 = O (1 − x ) n log n −1<α≤ (5.20)
2
and
n ro
2 r/2 r
I3 = O (1 − x ) n log n −1<α≤ . (5.21)
2
2n−2
X
2 r/2 r
| Bk (z) | = O (1 − x ) n log n , (5.22)
k=1
r−1
where −1 < α ≤ .
2
2n−2
X 2n−2
X 2n−2
X
| Bk (z) |≤ z −n
W (z)bk Jk (z) + z −n W (z)((b0k J10 (z) + b1k J11 (z)) .
k=1 k=1 k=1
| {z } | {z }
M1 M2
(5.24)
2n−2
X
| Bk (z) | ≤ M1 + M2 . (5.25)
k=1
2n−2
( )Z
z
X 1
M1 = z −n W (z) z n+1 (z 2 − 1)r Lk∗ (z)dz .
k=1
W (tk )tk (t2k − 1)r 0
2n−2
( )
X 1
M1 = z −n W (z) ×
k=1
W (tk )tk (t2k − 1)r
( )
W
Z z ∗
(z)
z n+1 (z 2 − 1)r dz
0 (z − tk )W ∗0 (tk )
2n−2
( )
z −n W (z) z
W ∗ (z)
X Z
≤ max z n+1 (z 2 − 1)r dz
k=1
W (tk )tk (t2k − 1)r |z|=1 0 (z − tk )W ∗0 (tk )
2n−2
z −n W (z)(z 2 − 1)r W ∗ (z) z
X Z
≤ z n+1 dz .
k=1
W (tk )tk (t2k − 1)r (z − tk )W ∗0 (tk ) 0
(α,β)0
Since Pn (x∗k ) = 0 and | tk |= 1, using (1.12) and (1.13) we get
0
X (1 − x2 )r/2 | Pn(α,β) (x) || Pn(α,β) (x) | 1 − xx∗k
2n−2 p
1
M1 ≤ r+1 (α,β) ∗ (α,β)00
.
(n + 2) k=1 (1 − x∗2
k ) 2 | P n (x k ) | | x − x ∗
k || P n (x∗
k ) |
For | x − x∗k | ≥ 1
2
| 1 − x∗2
k |, we have
2n
!
X 1
M1 = O (1 − x2 )r/2 nr .
k=1
k r−2α
From r − 2α ≥ 1, we get
n r − 1o
2 r/2 r
M1 = O (1 − x ) n log n −1<α≤ . (5.27)
2
The estimate remains the same in the case, where | x − x∗k | < 1
2
| 1 − x∗2
k |.
Similarly, we have
n r − 1o
2 r/2 r
M2 = O (1 − x ) n log n −1<α≤ . (5.28)
2
5.3 Convergence
Remark 5.1 :To prove the Theorem 5.3.1, we shall need following.
Let f (z) be a function continuous on T ∪ D and analytic on D. Then, there
exists a polynomial Fn (z) of degree ≤ 2n+1 satisfying Jackson’s inequality
[62].
| f (z) − Fn (z) | ≤ C ωr+1 (f, n−1 ), r≥0 (5.29)
where ωr+1 (f, n−1 ) denotes the (r + 1)th modulus of continuity of f (z) and
C is a constant independent of n and z.
2n
X 2n−2
X
Rn (z) = f (zk )Ak (z) + βk Bk (z), (5.32)
k=1 k=1
2 r/2 r+1 −1
| Rn (z) − f (z) | = O (1 − x ) n ωr+1 (f, n ) log n , (5.33)
2n 2n−2
0
X X
Fn (z) = Fn (zk )Ak (z) + Fn (zk )Bk (z). (5.34)
k=1 k=1
2n
X
| Rn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Ak (z) |
k=1
2n−2
X
+ | βk − Fn0 (zk ) || Bk (z) | + | Fn (z) − f (z) |
k=1
X2n 2n−2
X
≤ | f (zk ) − Fn (zk ) || Ak (z) | + | βk || Bk (z) |
k=1 k=1
| {z } | {z }
N1 N2
2n−2
X
+ | Fn0 (zk ) || Bk (z) | + | Fn (z) − f (z) | .
| {z }
k=1 N4
| {z }
N3
where
2n
X
N1 = | f (zk ) − Fn (zk ) || Ak (z) | .
k=1
2 r/2 −1 r+1
N1 = O (1 − x ) ωr+1 (f, n )n log n . (5.37)
2n−2
X
N2 = | βk || Bk (z) | .
k=1
2 r/2 −1 r
N2 = O (1 − x ) n ωr+1 (f, n ))n log n . (5.38)
2n−2
X
N3 = | Fn0 (zk ) || Bk (z) | .
k=1
2 r/2 −1 r
N3 = O (1 − x ) n ωr+1 (f, n )n log n . (5.39)
N4 = | Fn (z) − f (z) | .
−1
N4 = O ωr+1 (f, n ) . (5.40)
2 r/2 r+1 −1
| Rn (z) − f (z) | = O (1 − x ) n ωr+1 (f, n ) log n .
5.4 Conclusion
88
Interpolation on non-uniform nodes of the unit circle
Fn (tk ) = ηk
; k = 1(1)2n − 2,
(1)
Fn (zk ) = ϕk ; k = 1(1)2n, (6.1)
Fn(m) (±1) = 0
; m = 0(1)r,
(m)
where, ηk and ϕk are complex constants. Also, Fn (z) denotes mth
derivative of Fn (z). We shall write Fn (z) satisfying (6.1)
2n−2
X 2n
X
Fn (z) = ηk Ck (z) + ϕk Dk (z), (6.2)
k=1 k=1
For k = 1(1)2n
Dk (tj ) = 0
; j = 1(1)2n − 2,
0
Remark: The equations (6.5) and (6.6) have been developed while de-
riving out explicit representation of the interpolatory polynomial. Readers
can get the motivation to form such expression from the idea to maintain
the degree of the polynomial as well as simultaneously satisfy the conditions
required to form the fundamental polynomial.
Z z
Mk (z) = z n+2α (z 2 − 1)r−2α Lk∗ (z)dz. (6.5)
0
Z z
M1j (z) = z n+2α+j (z 2 − 1)r−2α−2 W (z)dz ; j = 0, 1 , (6.6)
0
Dk (z) = (z − 1)
2 2+2α −n−2α−1
z W (z) gk Mk (z) + g0k M10 (z) + g1k M11 (z)
∗
(6.7)
zk
gk = , (6.8)
(zk2 − 1)2+r W ∗ (zk )
and
gk [Mk (−1) − Mk (1)]
g0k = . (6.10)
2M10 (1)
where
0
z
z n+2α+1 (z 2 − 1)r−2−2α (z 2 − 1)W (z) + dk W (z)
Z
Nk (z) = − dz,
0 W ∗0 (tk )(t2k − 1)r+1 W (tk ) (z − tk )
(6.12)
Nk (1) − Nk (−1)
h0k = , (6.13)
2M10 (1)
(1 − tk 2 )W 0 (tk )
dk = . (6.15)
W (tk )
2n−2
2 r−2
X
| Ck (z) | = O (1 − x ) n log n ,
2
r
(6.16)
k=1
r−1
where −1 < α ≤ .
2
2n−2 2n−2
X X (z 2 − 1)r+1 Lk∗ (z)W (z)
| Ck (z) | ≤ (6.17)
k=1 k=1
(t2k − 1)r+1 W (tk )
| {z }
I1
2n−2
X
+ (z 2 − 1)2+2α z −n−2α−1 W ∗ (z)Nk (z)
k=1
| {z }
I2
2n−2
X
+ (z 2 − 1)2+2α z −n−2α−1 W ∗ (z)(h0k M10 (z) + h1k M11 (z)) ,
k=1
| {z }
I3
2n−2
X
| Ck (z) | ≤ I1 + I2 + I3 . (6.18)
k=1
We have
2n−2
X (z 2 − 1)r+1 Lk∗ (z)W (z)
I1 = .
k=1
(t2k − 1)r+1 W (tk )
2n−2
X (z 2 − 1)r+1 W ∗ (z)W (z)
I1 = .
k=1
(t2k − 1)r+1 (z − tk )W ∗0 (tk )W (tk )
2n
!
X 1
I1 = O (1 − x2 )r−2/2 nr−1 .
k=1
k r−2α
From r − 2α ≥ 1, we get
n r − 1o
2 r−2 r−1
I1 = O (1 − x ) 2 n log n −1<α≤ . (6.19)
2
The reader can verify that estimate remains the same in the case, where
| x − x∗k | < 1
2
| 1 − x∗2
k |.
2 r−2 r
I2 = O (1 − x ) n log n .
2 (6.20)
2 r−2 r
I3 = O (1 − x ) 2 n log n . (6.21)
2n
2 r−2
X
| Dk (z) | = O (1 − x ) 2 n log n ,
r−1
(6.22)
k=1
r−1
where, −1 < α ≤ .
2
| Dk (z) |= (z 2 − 1)2+2α z −n−2α−1 W ∗ (z) [gk Mk (z) + g0k M10 (z) + g1k M11 (z)] .
(6.23)
2n
X 2n
X
| Dk (z) | ≤ (z 2 − 1)2+2α z −n−2α−1 W ∗ (z)gk Mk (z) (6.24)
k=1 k=1
X2n
+ (z 2 − 1)2+2α z −n−2α−1 W ∗ (z)[g0k M10 (z) + g1k M11 (z)] ,
k=1
2n
X
| Dk (z) | ≤ I1 + I2 , (6.25)
k=1
where
2 r−2 r−1
I1 = O (1 − x ) n log n .
2 (6.26)
r−2
2 r−1
I2 = O (1 − x ) 2 n log n . (6.27)
6.3 Convergence
2n−2
X 2n
X
Fn (z) = f (zk )Ck (z) + ϕk Dk (z), (6.32)
k=1 k=1
r−2
| Fn (z) − f (z) | = O((1 − x2 ) 2 nr ωr+1 (f, n−1 ) log n), (6.33)
2n−2 2n
0
X X
Fn (z) = Fn (zk )Ck (z) + Fn (zk )Dk (z) (6.34)
k=1 k=1
Then
2n−2
X
| Fn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Ck (z) |
k=1
2n
X
+ {| ϕk | + | Fn0 (zk ) |} | Dk (z) | + | Fn (z) − f (z) |
k=1
Using equations (6.29), (6.30), (6.31), Lemma 6.2.1 and Lemma 6.2.2, we get
(6.33).
98
Interpolation on non-uniform nodes of the unit circle
(r)
where αk0 s are complex constants and Qn (zk ) denotes rth derivative of Qn (z)
at z = zk .
Polynomial
2n+1
X
Qn (z) = f (zk )A0k (z), (7.2)
k=0
where A0k (z) is a polynomial of degree less than (2n + 2)(r + 1) satisfying
the conditions given in (7.3).
For j, k = 0, 1, ..., 2n + 1 ,
A0k (zj ) = δkj ,
(7.3)
A(r)
0k (zj ) = 0 ; r = 1, 2, 3, 4 ,
such that
4
X
5
A0k (z) = [Lk (z)] + cpk Apk (z), (7.4)
p=1
5 L0k (zk )
c1k = − , (7.6)
R 0 (zk )
5 00 0 2
c2k = − Lk (zk ) + 10 (Lk (zk )) , (7.7)
2! [R 0 (zk )]2
5 00 0 000 0 3
c3k = − − 18 Lk (zk ) Lk (z k ) + Lk (zk ) − 198(Lk (zk )) , (7.8)
3! [R 0 (zk )]3
4
X
5
A0k (z) = [Lk (z)] + cpk [R(z)]p (Lk (z))5−p . (7.10)
p=1
At z = zj , where j = 0, 1, ..., 2n + 1 ,
4
X
5
A0k (zj ) = [Lk (zj )] + cpk [R(zj )]p (Lk (zj ))5−p .
p=1
4
hX i0
A00k (z) = 5L0k (z)[Lk (z)]4 + c1k [R(z)(Lk (z)) ] + 4 0
cpk [R(z)]p (Lk (z))5−p .
p=2
(7.12)
Clearly, at z = zj (j 6= k), we have A00k (z) = 0.
At z = zk , A00k (z) must be equal to zero. We have
als
p
[ 2 ] p−s
5 X X (s) (p−sr)
cpk = epsr [Lk (zk )]r Lk (zk ), (7.13)
p![R 0 (zk )]p s=0 r=s
p
where epsr are the constants independent of n and z and denotes
2
(s)
greatest integer function Also, Lk (zk ) denotes sth derivative of Lk (z) at
z = zk .
1
| Lk (z) |= O 3 (7.14)
k −α+ 2
R(z)
Lk (z) = . (7.15)
(z − zk )R 0 (zk )
| R(z) |
| Lk (z) |= .
| z − zk || R 0 (zk ) |
| (z 2 − 1)W (z) |
| Lk (z) |= .
| z − zk || {2zW (z) + (z 2 − 1)W 0 (z))}z=zk |
(α,β) 1+z 2
(z 2 − 1)Kn Pn 2z
zn
| Lk (z) | = n o0
(α,β) 1+z 2
| z − zk | (zk2 − 1) Kn Pn 2z
zn
2 z=zk
(α,β) 1+z
(z 2 − 1)Pn 2z
zn
= n 2 o
(α,β) (α,β)0 zk −1
| z − zk | (zk2 − 1) nzkn−1 Pn (xk ) + zkn Pn (xk ) 2z 2
k
(α,β)
2|z − 1||Pn (x)||z|n
2
= (α,β)0
.
|z− zk | |zk |n−2 |zk2 − 1|2 |Pn (xk )|
√ (α,β)
2.2 1 − x2 |Pn (x)|
| Lk (z) | = √ q √ (α,β)0
p
4(1 − x2k ) 2 1 − xxk − 1 − x2 1 − x2k |Pn (xk )|
√ (α,β)
q √ p
1 − x2 |Pn (x)| 1 − xxk + 1 − x2 1 − x2k
= √ p (α,β)0
2(1 − x2k ) (1 − xxk )2 − (1 − x2 )(1 − x2k )|Pn (xk )|
√ (α,β)
q p
1 − x2 |Pn (x)| 1 − xxk + 1 − x2 − x2k + x2 x2k
= √ (α,β)0
2(1 − x2k ) | x − xk | |Pn (xk )|
√ (α,β)
q p
1 − x2 |Pn (x)| 1 − xxk + (1 − xxk )2 − (x − xk )2
= √ (α,β)0
2(1 − x2k ) | x − xk | |Pn (xk )|
√ (α,β) √
1 − x2 |Pn (x)| 1 − xxk
≤ (α,β)0
.
(1 − x2k ) | x − xk | |Pn (xk )|
1
For | x − xk |≥ 2
| 1 − x2k |, we get
√ (α,β)
1 − x2 |Pn (x)|
| Lk (z) | ≤ C (α,β)0
,
(1 − x2k )3/2 |Pn (xk )|
1
Similarly, for | x − xk |≤ 2
| 1 − x2k |, we get the same result as (7.16).
For k = 0 and k = 2n + 1, we have
!
1
| cpk |= O . (7.18)
Knp np(α−1) k p/2−pα
!
Kn n−2 2 0
R (zk ) =
0
zk (zk − 1)2 Pn(α,β) (xk ).
2
1
| R 0 (zk ) |= O(Kn k −α+ 2 nα ). (7.19)
(s)
| Lk (zk ) |= O(ns ). (7.20)
Lemma 7.3.3. Let A0k (z) be given by (7.4) and cpk given by (7.13), then
for z ∈ T ∪ D,
2n+1
X
| A0k (z) | = O(log n), (7.21)
k=0
1
where −1 < α ≤ .
2
!
1 + z2 n
R(z) = Kn (z 2 − 1)Pn(α,β) z . (7.22)
2z
Taking modulus on both the sides and using (1.18) and (1.20), we have
1
For | xk − x | ≥ | 1 − x2k | and from (7.4) and (7.5), we have
2
2n+1
X 2n+1
X 2n+1
XX 4
| A0k (z) | = 5
| Lk (z) | + | cpk | | R(z) |p | Lk (z) |5−p .
k=0 k=0 k=0 p=1
Now, we get
2n+1
( )
X 13 p
| A0k (z) |= O(log n), −1<α≤ − . (7.24)
k=0
10 5
1
Similarly, for | xk − x |≤ 2
| 1 − x2k |, we get the same result.
Since range of α with p = 4 lies in the intersection of all the cases. Hence,
the lemma follows.
7.4 Convergence
where ω(f, n−1 ) represents the modulus of continuity of the function f (z).
where ω(f, n−1 ) denotes the first modulus of continuity of f (z) as well as C
is a constant independent of n and z.
Then,
2n+1
X
| Qn (z) − f (z) | ≤ | f (zk ) − Fn (zk ) || Ak (z) | + | Fn (z) − f (z) | .
| {z }
k=0 N2
| {z }
N1
N2 = O ω(f, n−1 ) .
(7.30)
Convergence of Lagrange
Interpolation on the
non-uniform nodes on the unit
circle
107
Interpolation on non-uniform nodes of the unit circle
2n+1
X
Ln (z) = f (zk )Lk (z) (8.2)
k=0
mial
2n+1
X
|Lk (z)| = O(log n), (8.4)
k=0
1
where −1 < α ≤ .
2
R(z)
Lk (z) = . (8.5)
(z − zk )R 0 (zk )
| R(z) |
| Lk (z) | = ,
| z − zk || R 0 (zk ) |
| (z 2 − 1)W (z) |
= .
| z − zk || {2zW (z) + (z 2 − 1)W 0 (z)}z=zk |
(α,β) 1+z 2
(z 2 − 1)Kn Pn 2z
zn
| Lk (z) | = n o0
(α,β) 1+z 2
| z − zk | (zk2 − 1) Kn Pn 2z
zn
2 z=zk
(α,β) 1+z
(z 2 − 1)Pn 2z
zn
= n 2 o
(α,β) (α,β)0 zk −1
| z − zk | (zk2 − 1) nzkn−1 Pn (xk ) + zkn Pn (xk ) 2z 2
k
(α,β)
2|z − 1||Pn (x)||z|n
2
= (α,β)0
.
|z− zk | |zk |n−2 |zk2 − 1|2 |Pn (xk )|
√ (α,β)
2.2 1 − x2 |Pn (x)|
| Lk (z) | = √ q √ (α,β)0
p
4(1 − x2k ) 2 1 − xxk − 1 − x2 1 − x2k |Pn (xk )|
√ (α,β)
q √ p
1 − x2 |Pn (x)| 1 − xxk + 1 − x2 1 − x2k
= √ p (α,β)0
2(1 − x2k ) (1 − xxk )2 − (1 − x2 )(1 − x2k )|Pn (xk )|
√ (α,β)
q p
1 − x2 |Pn (x)| 1 − xxk + 1 − x2 − x2k + x2 x2k
= √ (α,β)0
2(1 − x2k ) | x − xk | |Pn (xk )|
√ (α,β)
q p
1 − x2 |Pn (x)| 1 − xxk + 1 + x2 x2k − 2xxk + 2xxk − x2 − x2k
= √ (α,β)0
2(1 − x2k ) | x − xk | |Pn (xk )|
√ (α,β)
q p
1 − x2 |Pn (x)| 1 − xxk + (1 − xxk )2 − (x − xk )2
= √ (α,β)0
2(1 − x2k ) | x − xk | |Pn (xk )|
√ (α,β) √
1 − x2 |Pn (x)| 1 − xxk
≤ (α,β)0
.
(1 − x2k ) | x − xk | |Pn (xk )|
1
For | x − xk |≥ 2
| 1 − x2k |,
√ (α,β)
1 − x2 |Pn (x)|
| Lk (z) | ≤ C (α,β)0
,
(1 − x2k )3/2 |Pn (xk )|
1
Similarly, for | x − xk |≤| 1 − x2k |, we get the same result as (8.6).
2
For k = 0 and k = 2n + 1, we have
2n+1 2n+1
!
X X 1
|Lk (z)| = O 3 , (8.8)
k=0 k=0 k −α+ 2
2n+1
X 1
|Lk (z)| = O(log n), −1 < α ≤ .
k=0
2
8.3 Convergence
where ω(f, n−1 ) represents the modulus of continuity of the function f (z).
inequality[62].
| f (z) − Fn (z) | ≤ C ω(f, n−1 ) (8.10)
2n+1
X
Fn (z) = Fn (zk )Lk (z). (8.11)
k=0
Therefore,
2n+1
X
|Ln (z) − f (z)| ≤ |f (zk ) − Fn (zk )||Lk (z)| + |Fn (z) − f (z)| .
| {z }
k=0
| {z }
We have,
|Ln (z) − f (z)| ≤ N1 + N2 . (8.12)
2n+1 2n+1
!
X X
N1 = |f (zk ) − Fn (zk )||Lk (z)| = O ω(f, n−1 ) |Lk (z)|
k=0 k=0
1
= O(ω(f, n−1 ) log n), −1 < α ≤ . (8.13)
2
and
N2 = |Fn (z) − f (z)| = O(ω(f, n−1 )). (8.14)
• We detail a test function f (z) in each example and the Lagrange inter-
polating polynomial Ln (z) that interpolates the test function is always
plot in black.
2z
f (z) = .
z2+2
1.8 0.4
1.6
0.3
1.4
1.2 0.2
1.0
0.1
0.8
1 2 3 4 5 6 0 1 2 3 4 5 6
Figure 8.1: |f (eiθ )| and |L2 (eiθ )|; |L2 (z) − f (z)| for θ ∈ [0, 2π) and n = 2.
0.14
2.0
0.12
1.8
0.10
1.6
0.08
1.4
0.06
1.2
0.04
1.0
0.02
0.8
0.6 1 2 3 4 5 6 0 1 2 3 4 5 6
Figure 8.2: |f (eiθ )| and |L4 (eiθ )|; |L4 (z) − f (z)| for θ ∈ [0, 2π) and n = 4.
0.008
1.8
0.006
1.6
1.4
0.004
1.2
1.0 0.002
0.8
0.6 1 2 3 4 5 6 0 1 2 3 4 5 6
Figure 8.3: |f (eiθ )| and |L8 (eiθ )|; |L8 (z) − f (z)| for θ ∈ [0, 2π) and n = 8.
!
−1 1
Figure 8.4: |f (x+iy)| and |L16 (x+iy)|; |L16 (z)−f (z)| for x, y ∈ √ ,√
2 2
and n = 16.
f (z) = ez .
Case 1: Here z ∈ T i.e z = eiθ , θ ∈ [0, 2π) and we choose Jacobi’s poly-
1
nomial parameter as α = β = . We have plotted the graphic of absolute
2
values of f and its interpolant L8 (z) constructed with the 18 roots of T as
nodes. The results can be seen in the following Fig. 8.5. In right of the Fig.
8.5, graphic of the error term |L8 (z) − f (z)| can be seen with error less than
8 × 10−10
8. × 10-10
2.5
6. × 10-10
2.0
4. × 10-10
1.5
1.0 2. × 10-10
0.5
1 2 3 4 5 6
1 2 3 4 5 6
Figure 8.5: |f (eiθ )| and |L8 (eiθ )|; |L8 (z) − f (z)| for θ ∈ [0, 2π) and n = 8.
!
−1 1
Figure 8.6: |f (x + iy)| and |L8 (x + iy)|; |L8 (z) − f (z)| for x, y ∈ √ ,√
2 2
and n = 8.
We also took random values except this rectangular domain within the
unit circle and results were found to be same indicating Lagrange interpola-
tory polynomials on the projected nodes delivers a very good approximation
of the analytic functions within unit circle.
8.5 Discussions
Taking into account the Theorem 8.3.1, results obtained via numerical exam-
ples are similar to the expected ones as per the theory developed throughout
the chapter. Polynomials Ln (z) converges to the test function f (z), which
clearly indicates that Lagrange interpolatory polynomial formed using the
projected nodes of the zeros of Jacobi polynomials does a pretty good job in
approximating suitable class functions on the unit circle.
The experiments given in the section 8.4 were obtained by using personal
codes elaborated with software Wolfram Mathematica [134] without using
any extra precision. All computations were performed on the Macbook Pro
13 inch 2019 edition (1.4 GHz Quad-Core Intel Core i5 and 8 GB 2133 MHz
LPDDR3).
117
Interpolation on non-uniform nodes of the unit circle
[14] Bahadur S., (2012), Pál - type (0;0,2)-Interpolation on the unit circle;
Int. J. Math. Sci., 11(1-2), pp. 103 - 111.
[17] Bahadur S. and Shukla M., (2014), A New Kind of Hermite In-
terpolation, Adv. Inequal. Appl., http://scik.org/index.php/aia/
article/view/1265
[18] Bahadur S. and Shukla M., (2014), On the Derivative of Jacobi’s Poly-
nomial, Asian Journal Mathematics and Applications, Article ID ama
0163, pp. 1-6. https://scienceasia.asia/files/163.pdf
[21] Bahadur S. and Shukla M., (2015), Weighted Pál - type Interpolation
On The Unit Circle, Int. J. Emer. Tech. Comp. App. Sci., 12(2), pp.
135-139 http://iasir.net/IJETCASpapers/IJETCAS15-350.pdf
[22] Bahadur S. and Shukla M., (2014), (0,2)-Interpolation On The Unit Cir-
cle, International Journal of Advancements in Research and Technology,
3(1).
[25] Bahadur S., (2017) and Bano S., Modified Hermite Interpolation on
The Unit Circle, Annals of Pure and Applied Mathematics, 14(3), pp.
531-537.
[26] Burkett, J.; Varma, A.K., (1995), On Birkhoff interpolation (0;2) case.,
11(2), pp. 59-66. https://zbmath.org/?q=an%3A0930.41004
[28] Balázs J., Turan P., (1958), Notes On Interpolation III, Acta. Math.
Acad. Sci. Hung., 9, pp. 195-214. https://link.springer.com/
article/10.1007/BF02023872
[31] Balázs K., Kilgore T., (1993), Interpolation and Simultaneous mean
convergence of derivatives, Acta. Math. Hungar, 61, pp. 63-72. https:
//link.springer.com/article/10.1007/BF01872098
[32] Berriochoa E., Cachafeiro A., Martinez Brey E., (2011), Some im-
provements to the Hermite - Féjer interpolation on the unit cir-
cle and bounded interval, Computers and Mathematics with Ap-
plications, 61(4), pp. 1228–1240. https://www.sciencedirect.com/
science/article/pii/S0898122110009909
[34] Berriochoa E., Cachafeiro A., Dı́az J., Martı́nez Brey E., (2013), Rate of
convergence of Hermite- Fejér interpolation on the unit circle, Journal
of Applied Mathematics. Article ID 407128, 8 pages https://doi.org/
10.1155/2013/407128.
[35] Berriochoa E., Cachafeiro A., Dı́az J., (2014), Hermite interpolation on
the Unit Circle Considering up to the Second Derivative, ISRN Mathe-
matical Analysis, http://dx.doi.org/10.1155/2014/808519.
[36] Berriochoa E., Cachafeiro A., Garcı́a-Amor J.M., (2020),
Gibbs–Wilbraham phenomenon on Lagrange interpolation based
on analytic weights on the unit circle, J. Comput. Appl. Math., 365.
https://doi.org/10.1016/j.cam.2019.112376.
[37] Berriochoa E., Cachafeiro A., Rabade H.G. and Garcı́a-Amor J.M.,
(2022), Mechanical Models for Hermite Interpolation on the Unit Circle,
Mathematics, , 9(9). https://www.mdpi.com/2227-7390/9/9/1043
[38] Brück R., (1996), Lagrange interpolation in non–uniformly distributed
zeros on unit circle, Analysis, 16(3), pp. 273–282 https://www.
degruyter.com/document/doi/10.1524/anly.1996.16.3.273/html.
[39] Bruin Marcel G. de and Sharma A., (2003), Lacunary Pál - type
Interpolation and over – convergence, Computational and Function
Theory, 3(1), pp. 305-323. https://sci-hub.hkvisa.net/10.1007/
BF03321040.
[40] Cavaretta A.S., Sharma A. and Varga R. S., (1980), Interpolation in
the roots of unity : An extension of a theorem of J. L. Walsh, Re-
sults in Mathematics, 3(2), pp. 155-191. https://link.springer.com/
article/10.1007/BF03323355.
[41] Chen W. and Sharma A., (2004), Lacunary interpolation on some non-
uniformly distributed nodes on the unit circle, Annales Universitatis
Scientiarum Budapestinensis, 16, pp. 69-82. http://ac.inf.elte.hu/
Vol_016_1996/069.pdf
[42] Cheney E.W., (1966), Introduction to Approximation Theory, McGraw-
Hill Book Co., New York. https://www.worldcat.org/title/
introduction-to-approximation-theory/oclc/243547
[43] Daruis L. and González-vera P., (2001), A Note on Hermite -Fejér
Interpolation for the Unit Circle, Applied Mathematic Letters, 14,
pp. 997-1003. https://www.sciencedirect.com/science/article/
pii/S0893965901000787
[44] Davis P.J., (1963), Interpolation and Approximation, Blaise-
lell Pub. Co., Waltham Mass. https://www.worldcat.org/title/
interpolation-and-approximation/oclc/228351
[49] Fejér L., (1916), Über Interpolation, Gött. Nachr, 6, pp. 66–91. https:
//www.math.technion.ac.il/hat/people/fejer.html
[51] Freud G., (1958), Bemerkungen Uber die Kanverg enzeins interpolation
sverfahrens Von. P. Turan; Acta. Math. Acad. Sci. Hung., 9, pp. 337-341.
https://link.springer.com/article/10.1007/BF02020263
[53] Freud G. and Vértesi P., (1967), A new proof of A. F. Timan Ap-
proximation Theorem; Studia Sci. Math. Hung., 2, pp. 403-414. https:
//mathscinet.ams.org/mathscinet-getitem?mr=221165
[55] Freud G. and Sharma A., (1974), Some good sequences of Interpolatory
polynomials, Cand. J. of Maths. 26(1), pp. 233-246. https://www.
cambridge.org/core/services/aop-cambridge-core/content/
view/BCE1E28A193914D480419ED3F1C53A5C/S0008414X0004877Xa.
pdf/some-good-sequences-of-interpolatory-polynomials.pdf
[68] Joó I. and Pál L.G., (1996), Lacunary (0;0,1)-Interpolation on the roots
of Jacobi polynomials and their derivatives respectively I (Existence,
Explicit Formulae, Unicity); Analysis Math., 22, pp. 289-298.
[69] Kilgore T. and Prestin J., (1991), On the order of simultaneous ap-
proximation by Lagrange - Hermite Interpolation; Math. Nachr., 150,
pp. 143-150. https://onlinelibrary.wiley.com/doi/abs/10.1002/
mana.19911500111
[70] Kiš O., (1960), Remarks on Interpolation (Russian), Acta Math. Acad.
Sci. Hungar, 11, pp. 49-64.
[71] Kiš O. and Vértesi P., (1967), On a new Interpolation process; Ann.
Univ. Sci. Budapest Etövos Sect. Maths, 10, pp. 117-128.
[73] Krebsz A., (2004), Weighted (0,1,3)- Interpolation on the roots of the
classical orthogonal polynomials, Mathematica Pannonica, 15(1), pp.
21-35.
[74] Knoop H., Locher F., (1990), Hermite-Fejér type interpolation and
Korovkin’s theorem, Bulletin of the Australian Mathematical Society,
42(3), pp. 383-390. https://doi:10.1017/S0004972700028549.
[79] Mastroianni G., Milovanovic G.V. and Notarangelo I., (2012), An In-
terpolation Process of Lagrange - Hermite Type, De L’ Instinct Math.,
Nouvelleserie Tome, 91(105), pp. 163-175.
[80] Milovanovic G.V., Notarangelo I. and Paster P., (2016), Lagrange - Her-
mite Interpolation on the real axis, Calcolo, 53(2), pp. 235-261.
[81] Mathur K.K. and Sharma A., (1961), Some interpolatory properties of
Hermite Polynomials, Acta. Math. Acad. Sci. Hung., 12, pp. 193-207.
[83] Mathur K.K. and Srivastava R., (1995), Pál – type interpolation on
infinite interval, J. of Math. Anal. and Appl., 192, pp. 346-359.
[90] Min Guo Huo, (1992), Mean convergence of derivatives of the Hermite
Interpolation operator, Adv. in Math (China), 21(3), pp. 329-333.
[92] Pál L.G., (1975), A new modification of H-F Interpolation; Anals. Math.,
1, pp. 197-205.
[95] Prasad J., (1992), On the degree of approximation of the Hermite and
Hermite- Féjer Interpolation, International J. Math. Sci., 15(1) pp. 47-
56.
[96] Prasad J., (1992), On the order of convergence and mean convergence
of interpolatory process, Ganita, 43(1-2) pp. 181-204.
[98] Rainville R.D., (1963), Special Functions, Mac Millan Co., New York.
[99] Sakai R., (1991), Hermite - Féjer Interpolation, Approx. Theory, Collog.
Math. Soc. Janos Bolyai, North – Holland, Amsterdam, pp. 591-601.
[100] Saxena R.B. and Sharma A., (1958), On Some Interpolatory properties
of Legendre polynomial; Acta. Math. Acad. Sci. Hung., 9, pp. 345-358.
[104] Sharma A. and Leeming D., (1969), Lacunary Interpolation (0, n−1, n)-
case, Mathematical (cluj), 11(34), pp. 155-162.
[106] Shen Xie. Chang and Shuai Bin Peny, (1991), The order of approxima-
tion of Lagrange Interpolating polynomials based on asymptotic roots
of unity; J. Math., 11(3), pp. 287-297.
[109] Srivastava V., Mathur P. and Mathur N., (2011), A new kind of Pál-
type Interpolation II, Int. J. Contemp. Math. Sciences, 6(4-5), pp. 2237-
2246.
[110] Srivastava V., Mathur P. and Mathur N., (2013), A new kind of Pál -
type (0;0,1)-Interpolation; Int. J. of Math. and soft computing, 3, pp.
21-32.
[112] Surányi J. and Turán P., (1995), Notes On Interpolation I, Acta. Math.
Sci. Hung., 6, pp. 67-79.
[114] Szabadós J. and Varma A.K., (1985), On the almost everywhere diver-
gence of Lagrange Interpolation on infinite interval; Acad. Sci. Math.
Hung., 48(14), pp. 443-450.
[117] Szabadós J. and Varma A.K., (1995), On a convergent Pál - type (0,2)-
Interpolation process; Acta. Math. Hung., 66(4), pp. 301-326.
[118] Szegő, G., (1975), Orthogonal Polynomials, Amer. Math. Soc. Coll.,
23.
[120] Szilli L., (1985), A convergence theorem for the Pál method of Inter-
polation on the roots of Hermite polynomials; Analysis Math., 11, pp.
75-84.
[121] Szilli L., (1993), Weighted (0,2) Interpolation on the roots of Classical
Orthogonal Polynomials; Bull. of Allahabad Math. Soc., 8, pp. 1-10.
[127] Varma A.K. and Sharma A., (1961), Some interpolatory properties of
Tchebicheff polynomials, (0,2) case modified, Pub. Math. Debreen, 8,
pp. 336-349.
[130] Varma A.K. and Prasad J., (1990), Rate of convergence for Lacunary
Interpolation processes based on the roots of unity; Acta Math. Hung.,
55(3), pp. 211-217.
[131] Vecchia B.D. and Mastroianni G., (1992), The Hermite Interpolations,
Calcolo, 30(4), pp. 371-394.
[135] Xiang, S., He, G., (2015), The Fast Implementation of Higher Order
Hermite–Fejér Interpolation, SIAM Journal on Scientific Computing, 37,
pp. A1727–A1751. https://doi.org/10.1137/140971488
[138] Xie S., (1996), (0,1,3)*- Interpolation on the unit circle, Acta Math.
Sinica, 39(5), pp. 690-700.
[139] Xie S., (2001) Weighted least square convergence of Lagrange interpo-
lation on the unit circle, Approximation Theory and its Applications,
17, pp. 60–68. https://doi.org/10.1023/A:1015502428520
?????