Weighted Least Square Convergence of Lagrange Interpolation on the Unit Circle

In the paper, a result of Walsh and Sharma on least square convergence of Lagrange interpolation polynomials based on the n-th roots of unity is extended to Lagrange interpolation on the sets obtained by projecting vertically the zeros of (1-x)²=P ^(a,) _n(x),a>0,>0,(1-x)P^(a,) _n(x),a>0,>-1,(1+x)P P^(a,) _n(x),a>-1,0 and P^(a,) _n(x),a>-1,>-1, respectively, onto the unit circle, where P^(a,) _n(x),a>-1,>-1, stands for the n-th Jacobi polynomial. Moreover, a result of Saff and Walsh is also extended.     

Uniform and mean approximation by certain weighted polynomials,with applications

LetW(x):= exp(-{tiQ(x})), where, for example, Q(x) is even and convex onR, and Q(x)/logx → ∞ asx → ∞. A result of Mhaskar and Saff asserts that ifa _n =a _n (W) is the positive root of the equation $$n = ({2 \mathord{\left/ {\vphantom {2 \pi }} \right. \kern-\nulldelimiterspace} \pi })\int_0^1 {{{a_n xQ'(a_n x)} \mathord{\left/ {\vphantom {{a_n xQ'(a_n x)} {\sqrt {1 - x^2 } }}} \right. \kern-\nulldelimiterspace} {\sqrt {1 - x^2 } }}dx,}$$ then, given any polynomialP _n(x) of degree at mostn, the sup norm ofP _n(x)W(a _n x) overR is attained on [-1, 1]. In addition, any sequence of weighted polynomials {p _n (x)W(a _n x)} ₁ ^∞ that is uniformly bounded onR will converge to 0, for ¦x¦>1. In this paper we show that under certain conditions onW, a function g(x) continuous inR can be approximated in the uniform norm by such a sequence {p _n (x)W(a _n x)} ₁ ^∞ if and only if g(x)=0 for ¦x¦? 1. We also prove anL _p analogue for 0W(x)=exp(?|x| ^α ), when α >1. Further applications of our results are upper bounds for Christoffel functions, and asymptotic behavior of the largest zeros of orthogonal polynomials. A final application is an approximation theorem that will be used in a forthcoming proof of Freud's conjecture for |x| ^p exp(?|x| ^α ),α > 0,p > ?1.     

An even-order three-point boundary value problem on time scales

We study the even-order dynamic equation (−1)ⁿx^(Δ∇)n(t)=λh(t)f(x(t)), t∈[a,c] satisfying the boundary conditions x^(Δ∇)i(a)=0 and x^(Δ∇)i(c)=βx^(Δ∇)i(b) for 0?i?n−1. The three points a,b,c are from a time scale , where 0<β(b−a)<c−a for b∈(a,c), β>0, f is a positive function, and h is a nonnegative function that is allowed to vanish on some subintervals of [a,c] of the time scale.     

Asymptotic Properties of Orthogonal Polynomials with Respect to a Non-discrete Jacobi-Sobolev Inner Product

Let {Q _n ^(α,β) (x)} _n=0 ^∞ denote the sequence of polynomials orthogonal with respect to the non-discrete Sobolev inner product

$\langle f,g\rangle=\int_{-1}^{1}f(x)g(x)d\mu_{\alpha,\beta}(x)+\lambda\int_{-1}^{1}f'(x)g'(x)d\nu_{\alpha,\beta}(x)$

where λ>0 and d μ _α,β(x)=(x?a)(1?x)^α?1(1+x)^β?1 dx, d ν _α,β(x)=(1?x) ^α (1+x) ^β dx with aα,β>0. Their inner strong asymptotics on (?1,1), a Mehler-Heine type formula as well as some estimates of the Sobolev norms of Q _n ^(α,β) are obtained.

     

Finite Element Matrix Sequences: the Case of Rectangular Domains

In the present paper, we consider a preconditioning strategy for Finite Element (FE) matrix sequences {A _n (a)} _n discretizing the elliptic problem $$\left\{ \begin{gathered} A_a u \equiv ( - )^k \nabla ^k [a(x,y)\nabla ^k u(x,y)] = f(x,y),{ }(x,y) \in \Omega = (0,1)^2 , \hfill \\ \left. {\left( {\frac{{\partial ^s }}{{\partial v^s }}u(x,y)} \right)} \right|_{\partial \Omega } \equiv 0,{ }s = 0,...,k - 1,{ }^{^{^{^{^{^{(1)} } } } } } \hfill \\ \end{gathered} \right.$$ with a(x,y) being a uniformly positive function and ν denoting the unit outward normal direction. More precisely, in connection with preconditioned conjugate gradient (PCG) like methods, we define the preconditioning sequence: {P _n (a)} _n , P _n (a):= $$\widetilde D$$ _n ^1/2(a)A _n (1) $$\widetilde D$$ _n ^1/2(a), where $$\widetilde D$$ _n (a) is the suitable scaled main diagonal of A _n (a). In fact, under the mild assumption of Lebesgue integrability of a(x), the weak clustering at the unity of the corresponding preconditioned sequence is proved. Moreover, if a(x,y) is regular enough and if a uniform triangulation is considered, then the preconditioned sequence shows a strong clustering at the unity so that the sequence {P _n (a)} _n turns out to be a superlinear preconditioning sequence for {A _n (a)} _n . The computational interest is due to the fact that the computation with A _n (a) is reduced to computations involving diagonals and two-level Toeplitz structures {A _n (1)} _n with banded pattern. Some numerical experimentations confirm the efficiency of the discussed proposal.     

A method of constructing Krall's polynomials

We propose a method of constructing orthogonal polynomials P_n(x) (Krall's polynomials) that are eigenfunctions of higher-order differential operators. Using this method we show that recurrence coefficients of Krall's polynomials P_n(x) are rational functions of n. Let P_n^(a,b;M)(x) be polynomials obtained from the Jacobi polynomials P_n^(a,b)(x) by the following procedure. We add an arbitrary concentrated mass M at the endpoint of the orthogonality interval with respect to the weight function of the ordinary Jacobi polynomials. We find necessary conditions for the parameters a,b in order for the polynomials P_n^(a,b;M)(x) to obey a higher-order differential equation. The main result of the paper is the following. Let a be a positive integer and b⩾−1/2 an arbitrary real parameter. Then the polynomials P_n^(a,b;M)(x) are Krall's polynomials satisfying a differential equation of order 2a+4.     

A weight space invariant with respect to a singular linear operator

For the singular operator $$Su = \int_a^b {\frac{{K(x, s) u (s)}}{{s - x}}} ds$$ invariant weight spacesλ _α ^β , _p (u(x)∈λ _α ^β , _p if 1⁰,u (x) ρ (x)∈ H _β ⁰ , 2⁰.‖u‖L _p(ρ₀)<∞, ρ (x) = (x?a) (b ?x)^1+β, ρ₀(x)=(b?x)^α(p?1), 0<α, β<1,p>1H ₀ ^β is a Hölder space. Multiplicative inequalities of the type of Kh. Sh. Mukhtarov are also obtained.     

On two-sided and asymptotic estimates for the norms of embedding operators of into L q (dμ)

Explicit upper and lower estimates are given for the norms of the operators of embedding of , n ∈ ?, in L _q (dµ), 0 < q < ∞. Conditions on the measure µ are obtained under which the ratio of the above estimates tends to 1 as n → ∞, and asymptotic formulas are presented for these norms in regular cases. As a corollary, an asymptotic formula (as n → ∞) is established for the minimum eigenvalues λ_{1, n, β} , β > 0, of the boundary value problems (?d ²/dx ²) ⁿ u(x) = λ|x| ^β?1, x ∈ (?1, 1), u ^(k)(±1) = 0, k ∈ {0, 1, ..., n ? 1}.     

A polynomial dual of partitions

Let P(x) = Σ_i=0ⁿa_ixⁱ be a nonnegative integral polynomial. The polynomial P(x) is m-graphical, and a multi-graph G a realization of P(x), provided there exists a multi-graph G containing exactly P(1) points where a_i of these points have degree i for 0≤i≤n. For multigraphs G, H having polynomials P(x), Q(x) and number-theoretic partitions (degree sequences) π, ?, the usual product P(x)Q(x) is shown to be the polynomial of the Cartesian product G × H, thus inducing a natural product π? which extends that of juxtaposing integral multiple copies of ?. Skeletal results are given on synthesizing a multi-graph G via a natural Cartesian product G₁ × … × G_k having the same polynomial (partition) as G. Other results include an elementary sufficient condition for arbitrary nonnegative integral polynomials to be graphical.     

On a problem of Chebyshev and Markov

п. л. ЧЕБышЕВыМ БылА пОс тАВлЕНА И РЕшЕНА жАДА ЧА: пРИ пРОИжВОльНО жАДАННО М НА [?1,1] пОлОжИтЕльНОМ МНОг ОЧлЕНЕP _l (x) жАДАННОИ ст ЕпЕНИl НАИтИ пРИ кАжДОМn≧1 МНОгОЧлЕНР _n ^* (x) стЕпЕН Ип с кОЁФФИцИЕНтОМ п РИx ⁿ , РАВНыМ 1, кОтОРыИ ьВл ьлсь Бы МНОгО-ЧлЕНОМ, НАИМЕНЕЕ УклОНьУЩИМ сь От НУль с ВЕсОМP _l ^?1 (x) В МЕтРИкЕC[? 1,1]. А. А. МАРкОВ ОБОБЩИл ЁтО т РЕжУльтАт И пРИ тЕх ж Е УслОВИьх НАP _l (x) пОстРО Ил Дль≧1/2 МНОгОЧлЕНыP _n ^* (x) стЕп ЕНИ п, НАИМЕНЕЕ УклОНьУЩИЕсь От НУль с ВЕсОМP _l ^-1/2 (x). В ДАННОИ стАтьЕ УкАжы ВАЕтсь БОлЕЕ пРОстОИ, ЧЕМ В [3], спОсОБ пОстРОЕНИь МН ОгОЧлЕНОВP _n ^* (x), ДАУЩИх РЕшЕНИЕ жАД АЧИ МАРкОВА, И пРИ ЁтОМ, ВО-пЕРВых, УстАНОВлЕН О, ЧтО ВськИИ тАкОИ МНОгОЧлЕН МОжН О пРЕДстАВИть В ВИДЕ л ИНЕИНОИ кОМБИНАцИИ НЕ БОлЕЕ Ч ЕМ Ижl+1 МНОгОЧлЕНОВ ЧЕБышЕВ АT _j (x)=cos (jarc cosx): $$P_n^* (x) = \mathop \Sigma \limits_{k = 0}^l \gamma _k T_{|n - l + k|} (x)$$ , В кОтОРОИ кОЁФФИцИЕН тыγ _k НЕ жАВИсьт Отп И, ВО-ВтОРых, УкАжАН спОс ОБ ЁФФЕктИВНОгО РАжы с-кАНИь кОЁФФИцИЕНтОВγ _k пО М НОгОЧлЕНУР _l (х).     

Copositive pointwise approximation

We prove that if a functionf ∈C ⁽¹⁾ (I),I: = [?1, 1], changes its signs times (s ∈ ?) within the intervalI, then, for everyn > C, whereC is a constant which depends only on the set of points at which the function changes its sign, andk ∈ ?, there exists an algebraic polynomialP _n =P _n (x) of degree ≤n which locally inherits the sign off(x) and satisfies the inequality $$\left| {f\left( x \right) - P_n \left( x \right)} \right| \leqslant c\left( {s,k} \right)\left( {\frac{1}{{n^2 }} + \frac{{\sqrt {1 - x^2 } }}{n}} \right)\omega _k \left( {f'; \frac{1}{{n^2 }} + \frac{{\sqrt {1 - x^2 } }}{n}} \right), x \in I$$ , where ω _k (f′;t) is thekth modulus of continuity of the functionf’. It is also shown that iff ∈C (I) andf(x) ≥ 0,x ∈I then, for anyn ≥k ? 1, there exists a polynomialP _n =P _n (x) of degree ≤n such thatP _n (x) ≥ 0,x ∈I, and |f(x) ?P _n (x)| ≤c(k)ω _k (f;n ^?2 +n ^?1 √1 ?x ²),x ∈I.     

On convergence of a generalized Pál interpolation process

Let f∈C _[−1,1] ^″ (r≥1) and R_n(f,α,β,x) be the generalized Pál interpolation polynomials satisfying the conditions R_n(f,α,β,x_k)=f(x_k),R_n ^′(f,α,β,x_k)=f′(x_k)(k=1,2,…,n), where {x_k} are the roots of n-th Jacobi polynomial P_n(α,β,x),α,β>−1 and {x _k ^″ } are the roots of (1−x²)P_n″(α,β,x). In this paper, we prove that holds uniformly on [0,1]. In Memory of Professor M. T. Cheng Supported by the Science Foundation of CSBTB and the Natural Science Foundatioin of Zhejiang.     

A note on the Jacobian Conjecture

In this note, we show that, if the Druzkowski mappings F(X)=X+(AX)^∗3, i.e. F(X)=(x₁+(a₁₁x₁+?+a_1nx_n)³,…,x_n+(a_n1x₁+?+a_nnx_n)³), satisfies TrJ((AX)^∗3)=0, then where δ is the number of diagonal elements of A which are equal to zero. Furthermore, we show the Jacobian Conjecture is true for the Druzkowski mappings in dimension ?9 in the case .     

Preconditioned HSS methods for the solution of non-Hermitian positive definite linear systems and applications to the discrete convection-diffusion equation

Summary. We study the role of preconditioning strategies recently developed for coercive problems in connection with a two-step iterative method, based on the Hermitian skew-Hermitian splitting (HSS) of the coefficient matrix, proposed by Bai, Golub and Ng for the solution of nonsymmetric linear systems whose real part is coercive. As a model problem we consider Finite Differences (FD) matrix sequences {A_n(a,p)}_n discretizing the elliptic (convection-diffusion) problem with being a plurirectangle of R^d with a(x) being a uniformly positive function and p(x) denoting the Reynolds function: here for plurirectangle we mean a connected union of rectangles in d dimensions with edges parallel to the axes. More precisely, in connection with preconditioned HSS/GMRES like methods, we consider the preconditioning sequence {P_n(a)}_n, P_n(a):= D_n^1/2(a)A_n(1,0) D_n^1/2(a) where D_n(a) is the suitably scaled main diagonal of A_n(a,0). If a(x) is positive and regular enough, then the preconditioned sequence shows a strong clustering at unity so that the sequence {P_n(a)}_n turns out to be a superlinear preconditioning sequence for {A_n(a,0)}_n where A_n(a,0) represents a good approximation of Re(A_n(a,p)) namely the real part of A_n(a,p). The computational interest is due to the fact that the preconditioned HSS method has a convergence behavior depending on the spectral properties of {P_n^-1(a)Re(A_n(a,p))}_n {P_n^-1(a)A_n(a,0)}_n: therefore the solution of a linear system with coefficient matrix A_n(a,p) is reduced to computations involving diagonals and to the use of fast Poisson solvers for {A_n(1,0)}_n.Some numerical experimentations confirm the optimality of the discussed proposal and its superiority with respect to existing techniques.Mathematics Subject Classification (1991): 65F10, 65N22, 15A18, 15A12, 47B65     

Linear equations with unknowns from a multiplicative group whose solutions lie in a small number of subspaces

Let K be a field of characteristic 0 and let (K^*)ⁿ denote the n-fold Cartesian product of K^*, endowed with coordinatewise multiplication. Let Γ be a subgroup of (K^*)ⁿ of finite rank. We consider equations (*) a₁x₁ + … + a_nx_n = 1 in x = (x₁x_n)Γ, where a = (a₁,a_n)(K^*)ⁿ. Two tuples a, b(K^*)ⁿ are called Γ-equivalent if there is a uΓ such that b = u · a. Gy?ry and the author [Compositio Math. 66 (1988) 329-354] showed that for all but finitely many Γ-equivalence classes of tuples a(K^*)ⁿ, the set of solutions of (*) is contained in the union of not more than 2⁽ⁿ⁺¹! proper linear subspaces of Kⁿ. Later, this was improved by the author [J. reine angew. Math. 432 (1992) 177-217] to (n!)²ⁿ⁺². In the present paper we will show that for all but finitely many Γ-equivalence classes of tuples of coefficients, the set of non-degenerate solutions of (*) (i.e., with non-vanishing subsums) is contained in the union of not more than 2ⁿ proper linear subspaces of Kⁿ. Further we give an example showing that 2ⁿ cannot be replaced by a quantity smaller than n.     

On Extraction of Smooth Solutions of a Class of Almost Hypoelliptic Equations with Constant Power

A linear differential operator P(x, D) = P(x₁,... x _n , D₁,..., D _n ) = ∑_αγ_α(x)D^α with coefficients γ_α(x) defined in E ⁿ is called formally almost hypoelliptic in E ⁿ if all the derivatives D^ν_ξP(x, ξ) can be estimated by P(x, ξ), and the operator P(x, D) has uniformly constant power in Eⁿ. In the present paper, we prove that if P(x, D) is a formally almost hypoelliptic operator, then all solutions of equation P(x, D)u = 0, which together with some of their derivatives are square integrable with a specified exponential weight, are infinitely differentiable functions.     

Solution of the Szili type conjugate problem on the roots of the laguerre polynomials

Let ?1<α≤0 and let $$L_n^{(\alpha )} (x) = \frac{1}{{n!}}x^{ - \alpha } e^x \frac{{d^n }}{{dx^n }}(x^{\alpha + n} e^{ - x} )$$ be the generalizednth Laguerre polynomial,n=1,2,… Letx ₁,x ₂,…,x _n andx*₁,x*₂,…,x* _n?1 denote the roots ofL _n ^(α) (x) andL _n ^(α)′ (x) respectively and putx*₀=0. In this paper we prove the following theorem: Ify ₀,y ₁,…,y _n ?1 andy ₁ ^′ ,…,y _n ^′ are two systems of arbitrary real numbers, then there exists a unique polynomialP(x) of degree 2n?1 satisfying the conditions $$\begin{gathered} P\left( {x_k^* } \right) = y_k (k = 0,...,n - 1) \hfill \\ P'\left( {x_k } \right) = y_k^\prime (k = 1,...,n). \hfill \\ \end{gathered} $$ .     

Zeros of Jacobi polynomials Pn( α, β)$ P_{n}^{(\alpha ,\beta)} $, − 2 < α, β < − 1

The sequence of Jacobi polynomials \(\{P_{n}^{(\alpha ,\beta )}\}_{n = 0}^{\infty }\) is orthogonal on (??1,1) with respect to the weight function (1 ? x)^α(1 + x)^β provided α > ??1,β > ??1. When the parameters α and β lie in the narrow range ??2 < α, β < ??1, the sequence of Jacobi polynomials \(\{P_{n}^{(\alpha ,\beta )}\}_{n = 0}^{\infty }\) is quasi-orthogonal of order 2 with respect to the weight function (1 ? x)^α+?1(1 + x)^β+?1 and each polynomial of degree n,n ≥?2, in such a Jacobi sequence has n real zeros. We show that any sequence of Jacobi polynomials \(\{P_{n}^{(\alpha ,\beta )}\}_{n = 0}^{\infty }\) with ??2 < α, β < ??1, cannot be orthogonal with respect to any positive measure by proving that the zeros of Pn??1(α,β) do not interlace with the zeros of Pn(α,β) for any \(n \in \mathbb {N},\)n ≥?2, and any α,β lying in the range ??2 < α, β < ??1. We also investigate interlacing properties satisfied by the zeros of equal degree Jacobi polynomials Pn(α,β),Pn(α,β+?1) and Pn(α+?1,β+?1) where ??2 < α, β < ??1. Upper and lower bounds for the extreme zeros of quasi-orthogonal order 2 Jacobi polynomials Pn(α,β) with ??2 < α, β < ??1 are derived.     

On the equivalence of the Euler-Pommier operators in spaces of analytic functions

In the space A (θ) of all one-valued functions f(z) analytic in an arbitrary region G ? ? (0 ∈ G) with the topology of compact convergence, we establish necessary and sufficient conditions for the equivalence of the operators L ₁=α _n z ⁿ Δ ⁿ + ... + α₁ zΔ+α₀ E and L ₂= z ⁿ a _n (z)Δ ⁿ + ... + za ₁(z)Δ+a ₀(z)E, where δ: (Δ?)(z)=(f(z)-?(0))/z is the Pommier operator in A(G), n ∈ ?, α _n ∈ ?, a _k (z) ∈ A(G), 0≤k≤n, and the following condition is satisfied: Σ _j=s ^n?1 α _j+1 ∈ 0, s=0,1,...,n?1. We also prove that the operators z ^s+1Δ+β(z)E, β(z) ∈ A _R , s ∈ ?, and z ^s+1 are equivalent in the spaces A _R, 0?R?-∞, if and only if β(z) = 0.