The umbral calculus and orthogonal polynomials

We determine all orthogonal polynomials having Boas-Buck generating functions g(t)(xf(t)), where% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacqqHOo% qwcaGGOaGaamiDaiaacMcacqGH9aqpruqqYLwySbacfaGaa8hiamaa% BeaaleaacaaIWaaabeaakiaadAeacaqGGaWaaSbaaSqaaiaabgdaae% qaaOGaaeikaiaadggacaGGSaGaa8hiaiaadshacaqGPaGaaeilaiaa% bccacaqGGaGaaeiiaiaadggacqGHGjsUcaaIWaGaaiilaiaa-bcacq% GHsislcaaIXaGaaiilaiaa-bcacqGHsislcaaIYaGaaiilaiablAci% ljaacUdaaeaacqqHOoqwcaGGOaGaamiDaiaacMcacqGH9aqpcaWFGa% WaaSraaSqaaiaaicdaaeqaaOGaamOraiaabccadaWgaaWcbaGaaeOm% aaqabaGccaGGOaWaaSqaaSqaaiaaigdaaeaacaaIZaaaaOGaaiilai% aa-bcadaWcbaWcbaGaaGOmaaqaaiaaiodaaaGccaGGSaGaa8hiaiaa% dshacaGGPaGaa8hiamaaBeaaleaacaaIWaaabeaakiaadAeacaqGGa% WaaSbaaSqaaiaabkdaaeqaaOGaaeikamaaleaaleaacaaIYaaabaGa% aG4maaaakiaacYcacaWFGaWaaSqaaSqaaiaaisdaaeaacaaIZaaaaO% Gaaiilaiaa-bcacaWG0bGaaiykaiaacYcacaWFGaWaaSraaSqaaiaa% icdaaeqaaOGaamOraiaabccadaWgaaWcbaGaaeOmaaqabaGccaGGOa% WaaSqaaSqaaiaaisdaaeaacaaIZaaaaOGaaiilaiaa-bcadaWcbaWc% baGaaGynaaqaaiaaiodaaaGccaGGSaGaa8hiaiaadshacaGGPaGaai% 4oaaqaaiabfI6azjaacIcacaWG0bGaaiykaiabg2da9iaa-bcadaWg% baWcbaGaaGimaaqabaGccaWGgbGaaeiiamaaBaaaleaacaqGZaaabe% aakiaacIcadaWcbaWcbaGaaGymaaqaaiaaisdaaaGccaGGSaGaa8hi% amaaleaaleaacaaIYaaabaGaaGinaaaakiaacYcacaWFGaWaaSqaaS% qaaiaaiodaaeaacaaI0aaaaOGaaiilaiaa-bcacaWG0bGaaiykaiaa% -bcadaWgbaWcbaGaaGimaaqabaGccaWGgbGaaeiiamaaBaaaleaaca% qGZaaabeaakiaabIcadaWcbaWcbaGaaGOmaaqaaiaaisdaaaGccaGG% SaGaa8hiamaaleaaleaacaaIZaaabaGaaGinaaaakiaacYcacaWFGa% WaaSqaaSqaaiaaiwdaaeaacaaI0aaaaOGaaiilaiaa-bcacaWG0bGa% aiykaiaacYcaaeaadaWgbaWcbaGaaGimaaqabaGccaWGgbGaaeiiam% aaBaaaleaacaqGZaaabeaakiaacIcadaWcbaWcbaGaaG4maaqaaiaa% isdaaaGccaGGSaGaa8hiamaaleaaleaacaaI1aaabaGaaGinaaaaki% aacYcacaWFGaWaaSqaaSqaaiaaiAdaaeaacaaI0aaaaOGaaiilaiaa% -bcacaWG0bGaaiykaiaacYcacaGGUaGaa8hiamaaBeaaleaacaaIWa% aabeaakiaadAeacaqGGaWaaSbaaSqaaiaabodaaeqaaOGaaeikamaa% leaaleaacaaI1aaabaGaaGinaaaakiaacYcacaWFGaWaaSqaaSqaai% aaiAdaaeaacaaI0aaaaOGaaiilaiaa-bcadaWcbaWcbaGaaG4naaqa% aiaaisdaaaGccaGGSaGaa8hiaiaadshacaGGPaGaaiOlaaaaaa!C1F3!\[\begin{gathered}\Psi (t) = {}_0F{\text{ }}_{\text{1}} {\text{(}}a, t{\text{), }}a \ne 0, - 1, - 2, \ldots ; \hfill \\\Psi (t) = {}_0F{\text{ }}_{\text{2}} (\tfrac{1}{3}, \tfrac{2}{3}, t) {}_0F{\text{ }}_{\text{2}} {\text{(}}\tfrac{2}{3}, \tfrac{4}{3}, t), {}_0F{\text{ }}_{\text{2}} (\tfrac{4}{3}, \tfrac{5}{3}, t); \hfill \\\Psi (t) = {}_0F{\text{ }}_{\text{3}} (\tfrac{1}{4}, \tfrac{2}{4}, \tfrac{3}{4}, t) {}_0F{\text{ }}_{\text{3}} {\text{(}}\tfrac{2}{4}, \tfrac{3}{4}, \tfrac{5}{4}, t), \hfill \\{}_0F{\text{ }}_{\text{3}} (\tfrac{3}{4}, \tfrac{5}{4}, \tfrac{6}{4}, t),. {}_0F{\text{ }}_{\text{3}} {\text{(}}\tfrac{5}{4}, \tfrac{6}{4}, \tfrac{7}{4}, t). \hfill \\\end{gathered}\]We also determine all Sheffer polynomials which are orthogonal on the unit circle. The formula for the product of polynomials of the Boas-Buck type is obtained.     

On the error of compound quadrature formulas forr-convex functions

LetC _m be a compound quadrature formula, i.e.C _m is obtained by dividing the interval of integration [a, b] intom subintervals of equal length, and applying the same quadrature formulaQ _n to every subinterval. LetR _m be the corresponding error functional. Iff ^(r) > 0 impliesR _m [f] > 0 (orR _m [f] < 0),=" then=" we=" say=">C _m is positive definite (or negative definite, respectively) of orderr. This is the case for most of the well-known quadrature formulas. The assumption thatf ^(r) > 0 may be weakened to the requirement that all divided differences of orderr off are non-negative. Thenf is calledr-convex. Now letC _m be positive definite or negative definite of orderr, and letf be continuous andr-convex. We prove the following direct and inverse theorems for the errorR _m [f], where , denotes the modulus of continuity of orderr:

Integrals involving Gegenbauer and Hermite polynomials on the imaginary axis

Summary The integral _- [C _2n (it)]^–2(1+t ²)^–-1/2 dt is evaluated for > –1/2 whereC _2n is the Gegenbauer polynomial of degree 2n. Letting gives the value _- [H _2n (it)]^–2 e ^{1-1/2t 2} dt involving the Hermite polynomialH _2n of degree 2n. The result is obtained using Gegenbauer functions of the second kind.     

Numerical solution of a singularly perturbed problem via exponential splines

It is proved that a spline difference scheme for a singularly perturbed self-adjoint problem, derived by using exponential cubic splines at mid-points, has second order uniform convergence in a small parameter . Numerical experiments are presented to confirm the theoretical predictions.     

Quadratic spline interpolation on uniform meshes

The interpolation problem at uniform mesh points of a quadratic splines(x _i)=f _i,i=0, 1,...,N ands(x ₀)=f₀ is considered. It is known that s–f=O(h ³) and s–f=O(h ²), whereh is the step size, and that these orders cannot be improved. Contrary to recently published results we prove that superconvergence cannot occur for any particular point independent off other than mesh points wheres=f by assumption. Best error bounds for some compound formulae approximatingf _i andf _i ⁽³⁾ are also derived.     

Comparison theorems for weak splittings of bounded operators

Summary Comparison theorems for weak splittings of bounded operators are presented. These theorems extend the classical comparison theorem for regular splittings of matrices by Varga, the less known result by Wonicki, and the recent results for regular and weak regular splittings of matrices by Neumann and Plemmons, Elsner, and Lanzkron, Rose and Szyld. The hypotheses of the theorems presented here are weaker and the theorems hold for general Banach spaces and rather general cones. Hypotheses are given which provide strict inequalities for the comparisons. It is also shown that the comparison theorem by Alefeld and Volkmann applies exclusively to monotone sequences of iterates and is not equivalent to the comparison of the spectral radius of the iteration operators.This work was supported by the National Science Foundation grants DMS-8807338 and INT-8918502     

On a problem proposed by H. Braß concerning the remainder term in quadrature for convex functions

Summary We prove that the error inn-point Gaussian quadrature, with respect to the standard weight functionw1, is of best possible orderO(n ^–2) for every bounded convex function. This result solves an open problem proposed by H. Braß and published in the problem section of the proceedings of the 2. Conference on Numerical Integration held in 1981 at the Mathematisches Forschungsinstitut Oberwolfach (Hämmerlin 1982; Problem 2). Furthermore, we investigate this problem for positive quadrature rules and for general product quadrature. In particular, for the special class of Jacobian weight functionsw _, (x)=(1–x)(1+x), we show that the above result for Gaussian quadrature is not valid precisely ifw _, is unbounded.Dedicated to Prof. H. Braß on the occasion of his 55th birthday     

Parallel and superfast algorithms for Hankel systems of equations

Summary Utilizing kernel structure properties a unified construction of Hankel matrix inversion algorithms is presented. Three types of algorithms are obtained: 1)O(n ²) complexity Levinson type, 2)O (n) parallel complexity Schur-type, and 3)O(n log² n) complexity asymptotically fast ones. All algorithms work without additional assumption (like strong nonsingularity).     

Nonlinear stability and convergence of finite-difference methods for the “good” Boussinesq equation

Summary The good Boussinesq equationu _tt =–u _xxxx +u _xx +(u ²) _xx has recently been found to possess an interesting soliton-interaction mechanism. In this paper we study the nonlinear stability and the convergence of some simple finite-difference schemes for the numerical solution of problems involving the good Boussinesq equation. Numerical experimentas are also reported.     

Fully discrete Galerkin approximations of parabolic boundary-value problems with nonsmooth boundary data

Summary In this paper we study the convergence properties of a fully discrete Galerkin approximation with a backwark Euler time discretization scheme. An approach based on semigroup theory is used to deal with the nonsmooth Dirichlet boundary data which cannot be handled by standard techniques. This approach gives rise to optimal rates of convergence inL _p[O,T;L ₂()] norms for boundary conditions inL _p[O,T;L ₂()], 1p.