Near-minimax Interpolation by a Polynomial in z and z-1 on a Circular Annulus

A polynomial of degree n in z^–1 and n–1 in z isdefined by an interpolation projection from the space A(N_p) of functions f analytic in the circular annulusp^–1 < <p and continuous on itsboundaries = p^–1, p. The points ofinterpolation are chosen to be spaced at equal angles aroundthe two boundaries, with arguments on the inner boundary midwaybetween those on the outer boundary. By calculating the Lebesgueconstants numerically, is found to be close to a minimax approximation for all p 1and all degrees n in the range 1 n 15. In the limiting casesp = 1 and, it is proved that is asymptotic to 2^–1 log n. More specifically and , where _nis the Lebesgueconstant of Gronwall for equally spaced interpolation on a circleby a polynomial of degree n. It is also demonstrated that is not in general monotonic in p, and that is not everywhere differentiable in p.     

On Mason's conjecture concerning interpolation by polynomials in z and z-1 on an annulus

A polynomial of degree n in z^–1 and n – 1 in z isdefined by an interpolation projection from the space A(N) of functions f analytic in thecircular annulus ^–1<|z| and continuous on its boundaries|z|=^–1, . The points of interpolation are chosen to coincidewith the n roots of zⁿ=–n and the n roots of zⁿ=^–n.We prove Mason's conjecture that the corresponding Lebesguefunction attains its maximal value on the inner circle. We alsoestimate the bound of the Lebesgue constant . It is proved that the following estimate for theoperator norm holds: where _n, is the Lebesgue constant of Gronwall for equally spacedinterpolation on a circle by a polynomial of degree n.     

An Efficient Enthalpy-type Method for the Stefan Problem

The paper presents a new finite-difference method for solvingthe one-dimensional two-phase Stefan problem. Under assumptionson the data which guarantee the temperature u and the movingboundary s to belong to and , respectively, we obtain L₂-errorestimates of order O(h + h^–?) provided the time step is chosen such that Numerical aspects are discussed.     

Finite-Element Approximation of a Plasma Equilibrium Problem

The plasma problem studied is: given R⁺ find (, d, u) R ?R ? H¹() such that Let ₁ < ₂ be the first two eigenvalues of the associatedlinear eigenvalue problem: find $$\left(\lambda ,\phi \right)\in\mathrm{R;}\times {\hbox{ H }}_{0}^{1}\left(\Omega \right)$$such that For ₀(0,₂) it is well known that there exists a unique solution(₀, d₀, u₀) to the above problem. We show that the standard continuous piecewise linear Galerkinfinite-element approximatinon $$\left({\lambda }_{0},{\hbox{d }}_{0}^{k},{u}_{0}^{h}\right)$$, for ₀(0,₂), converges atthe optimal rate in the H¹, L², and L norms as h, the mesh length,tends to 0. In addition, we show that dist (, ^h)Ch² ln 1/h,where $${\Gamma }^{\left(h\right)}=\left\{x\in \Omega :{u}_{0}^{\left(h\right)}\left(x\right)=0\right\}$$.Finally we consider a more practical approximation involvingnumerical integration.     

Generalized compound quadrature formulae for finite-part integrals

Received on 31 July 1995. Revised on 19 August 1996. We investigate the error term of the dth degree compound quadratureformulae for finite-part integrals of the form where and p 1.We are mainly interested in error bounds of the form with best possible constants c. Itis shown that, for and n uniformlydistributed nodes, the error behaves as O(n^p–s–1for , p–1 <s d+1.In a previous paper we have shown that this is not true for As an improvement, we consider the case of non-uniformly distributednodes. Here, we show that for all p I and , an O(n^–s) error estimate can be obtainedin theory by a suitable choice of the nodes. A set of nodeswith this property is staled explicitly. In practice, this gradedmesh causes stability problems which are computationally expensiveto overcome. E-mail address: diethelm{at}informatik.uni-hildesheim.de     

Numerical Solvability of Hammerstein Integral Equations of Mixed Type

We consider a mixed Hammerstein integral equation of the form where –<a<b<, y, f_i and k_i, (1im) are known functionsand x is a solution to be determined. In this paper, we obtainexistence, uniqueness, and numerical solvability of (I) undercertain smoothness assumptions on the known functions y, f_iand k_i.     

An Order of Convergence for Some Radial Basis Functions

The quasi-interpolant to a function f : RⁿR on an infinite regulargrid of spacing h can be defined by where : RⁿR is a function which decays quickly for large argument.In the case of radial basis functions has the form where : R⁺R is known as a radial basis function and, in general,?_j R (j = 1,...,m) and x^j Rⁿ (j = 1,...,m), though here onlythe particular case x^j Zⁿ (j = 1,..., m) is considered. Thispaper concentrates on the case (r) = r, a generalization oflinear interpolation, although some of the analysis is moregeneral. It is proved that, if n is odd, then there is a function such that the maximum difference between a sufficiently smoothfunction and its quasi-interpolant is bounded by a constantmultiple of hⁿ⁺¹. This is done by first showing that such aquasi-interpolation formula can reproduce polynomials of degreen.     

Inverse Spectral Problems for Sturm-Liouville Equations with Eigenparameter Dependent Boundary Conditions

Inverse Sturm–Liouville problems with eigenparameter-dependentboundary conditions are considered. Theorems analogous to thoseof both Hochstadt and Gelfand and Levitan are proved. In particular, let ly = (1/r)(–(py')'+qy), , where det = > 0, c 0, det > 0, t 0 and (cs + dr –au – tb)² < 4(cr – ta)(ds – ub). Denoteby (l; ; ) the eigenvalue problem ly = y with boundary conditionsy(0)cos+y'(0)sin = 0 and (a+b)y(1) = (c+d)(py')(1). Define (; ; ) as above but with l replacedby . Let w_n denote the eigenfunctionof (l; ; ) having eigenvalue n and initial conditions w_n(0)= sin and pw'_n(0) = –cos and let _n = –aw_n(1)+cpw'_n(1).Define _n and _n similarly. As sample results, it is proved that if (l; ; ) and (; ; ) have the same spectrum, and (l;; ) and (; ; ) have the samespectrum or for all n, thenq/r = /.     

Finite-Element Approximation of Elliptic Equations with a Neumann or Robin Condition on a Curved Boundary

This paper considers a finite-element approximation of a second-orderself adjoint elliptic equation in a region Rⁿ (with n=2 or 3)having a curved boundary on which a Neumann or Robin conditionis prescribed. If the finite-element space defined over , a union of elements, has approximation power h^kin the L² norm, and if the region of integration is approximatedby ^h with dist (, ^h)Ch^k, then it is shown that one retains optimalrates of convergence for the error in the H¹ and L² norms, whetherQ^h is fitted or unfitted , provided that the numerical integration scheme has sufficientaccuracy.     

Weighted Integral Inequalities for the Hardy Type Operator and the FractionalMaximal Operator

Let w(x), u(x) and (x) be weight functions. In this paper, underappropriate conditions on Young's functions ₁, ₂ we characterizethe inequality for the Hardy-typeoperator T defined in [1] and the inequality for the fractional maximal operator M_{, ;} definedin [8], as well as the corresponding weak-type inequalities.     

Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: the scalar case

^** Email: Paul.Houston{at}mcs.le.ac.uk^*** Email: Janice.Robson{at}comlab.ox.ac.uk^**** Email: Endre.Suli{at}comlab.ox.ac.uk We develop a one-parameter family of hp-version discontinuousGalerkin finite element methods, parameterised by [–1,1], for the numerical solution of quasilinear elliptic equationsin divergence form on a bounded open set ^d, d 2. In particular,we consider the analysis of the family for the equation –·{µ(x, |u|)u} = f(x) subject to mixed Dirichlet–Neumannboundary conditions on . It is assumed that µ is a real-valuedfunction, µ C( x [0, )), and thereexist positive constants m_µ and M_µ such that m_µ(t– s) µ(x, t)t – µ(x, s)s M_µ(t– s) for t s 0 and all x . Using a result from the theory of monotone operators for any valueof [–1, 1], the corresponding method is shown to havea unique solution u_DG in the finite element space. If u C¹() H^k(), k 2, then with discontinuous piecewise polynomials ofdegree p 1, the error between u and u_DG, measured in the brokenH¹()-norm, is (h^s–1/p^k–3/2), where 1 s min {p+ 1, k}.     

Rate of convergence of numerical approximations to homoclinic bifurcation points

For x=f (x, ), x Rⁿ, R, having a hyperbolic or semihyperbolicequilibrium p(), we study the numerical approximation of parametervalues ^* at which there is an orbit homoclinic to p(). We approximate^* by solving a finite-interval boundary value problem on J=[T_–,T₊], T_–<0<T₊, with boundary conditions that sayx(T_–) and x(T₊) are in approximations to appropriate invariantmanifolds of p(). A phase condition is also necessary to makethe solution unique. Using a lemma of Xiao-Biao Lin, we improve,for certain phase conditions, existing estimates on the rateof convergence of the computed homoclinic bifurcation parametervalue , to the true value ^*. The estimates we obtain agree withthe rates of convergence observed in numerical experiments.Unfortunately, the phase condition most commonly used in numericalwork is not covered by our results.     

Almost Everywhere Convergence of Bochner-Riesz Means On The Heisenberg Group and Fractional Integration On The Dual

Let L denote the sub-Laplacian on the Heisenberg group H_n and the corresponding Bochner-Riesz operator. Let Q denote the homogeneous dimension and D the Euclideandimension of H_n. We prove convergence a.e. of the Bochner-Rieszmeans as r 0 for > 0and for all f L^p(H_n), provided that . Our proof is based on explicit formulas for the operators with a C, defined on the dual ofH_n by , which may be of independent interest. Here is given by for all (z,u) H_n. 2000 Mathematical Subject Classification: 22E30, 43A80.     

Best Approximation with Coefficient Constraints

For given = (₁,..., _n) and ß = (ß₁,...,ß_n), with – _i < ß_i (i = 1, ...,n) and continuous functions u₁,...,u_n, set This paper is concerned with best approximating continuous functions,in the uniform norm, from U(; ß). We exactly characterizethe u₁,..., u_n for which the best approximant to every continuousfunction is unique. We also present a general theorem characterizingall best approximants. When (u₁,..., u_n) is a Descartes, ora weak Descartes, system on [0, 1], explicit characterizationsof the best approximants in terms of equioscillations are given.These results are applied to spline spaces. They are also usedto complete the characterizations in certain specific examplespreviously considered in the literature.     

The Decay of Eddy-currents in Thin Sheets and Related Water-wave Problems

The decay of the eddy-currents that are induced in a thin, uniform,imperfectly-conducting sheet by switching off the source ofan external magnetic field is investigated. For the two-dimensionalproblem of an infinite strip the (non-dimensional) decay constants_n and eddy-current distributions i_n(x) are the eigenvalues andeigenfunctions of the integral equation with the constraint. For the circular disc the corresponding equation is where and K and E are complete elliptic integrals. For both problemsthe initial eddy-currents have inverse-square-root singularitiesat the edges but during their decay the eddy currents are finiteat the edges and the normal magnetic fields have logarithmicsingularities there. Numerical results are given for variousinitial-value problems. The eddy current problems are closely related to water-waveproblems in which there is a strip-shaped or circular aperturein a horizontal rigid dock. If _n and _n are the decay constantsand magnetic scalar potentials for the strip and _n and _n theangular frequencies and velocity potentials for the normal modesin the strip-shaped aperture, then _n =_n² and _n and _n are thereal and imaginary parts respectively of a holomorphic function.The velocities in the normal modes are deduced from the solutionof the eddy-current problem and are found to agree with resultsgiven in Miles (1972). For circular geometries the eigenvaluesand eigenfunctions of the axisymmetric eddy-current problemare the same as those of the water-wave problem that has angularvariation eⁱ; where (, , z) are cylindrical polar co-ordinateslocated at the centre of the basin.     

On Artin's Braid Group And Polyconvexity In The Calculus Of Variations

Let R² be a bounded Lipschitz domain and let be a Carathèodory integrand such that F(x,·) is polyconvex for L²-a.e. x . Moreover assume thatF is bounded from below and satisfies the condition as det for L²-a.e. x . The paper describes the effect of domain topologyon the existence and multiplicity of strong local minimizersof the functional wherethe map u lies in the Sobolev space W_id^1,p (, R²) with p 2and satisfies the pointwise condition u(x) >0 for L²-a.e.x . The question is settled by establishing that F[·]admits a set of strong local minimizers on that can be indexed by the group P_n Zⁿ, the directsum of Artin's pure braid group on n strings and n copies ofthe infinite cyclic group. The dependence on the domain topologyis through the number of holes n in and the different mechanismsthat give rise to such local minimizers are fully exploitedby this particular representation.     

Recurrence and asymptotics for orthonormal rational functions on an interval

^{Let µ be a positive bounded Borel measure on a subsetI of the real line and = {₁, ..., _n} a sequence of arbitrary ‘complex’poles outside I. Suppose {₁, ..., _n} is the sequence of rationalfunctions with poles in orthonormal on I with respect to µ. First, we are concernedwith reducing the number of different coefficients in the three-termrecurrence relation satisfied by these orthonormal rationalfunctions. Next, we consider the case in which I = [–1, 1] and µ satisfies the Erdos–Turán conditionµ' > 0 a.e. on I (where µ' is the Radon–Nikodymderivative of the measure µ with respect to the Lebesguemeasure) to discuss the convergence of _n+1(x)/_n(x) as n tendsto infinity and to derive asymptotic formulas for the recurrencecoefficients in the three-term recurrence relation. Finally,we give a strong convergence result for _n(x) under the morerestrictive condition that µ satisfies the Szeg condition(1 – x2)–1/2 log µ'(x) L1([– 1, 1]).}     

Localization of Closed (or Periodic) Solutions of a Differential System with Concave Nonlinearities

Consider a scalar differential equation , where I is an open interval containing [0,T]. Assumethat f(t, x) is continuous with a continuous derivative , and weakly concave (or weakly convex)in x for all t I, though strictly concave (or strictly convex)for some t [0, T]. It is well known that in this case therecan be either no, one or two closed solutions; that is, solutions(t) for which (0) = (T) If there are two closed solutions, thenthe greater has a negative characteristic exponent and the smallerhas a positive one. It is easily seen that this is equivalentto a statement on localization of closed solutions. It is shownhow this statement can be generalized to systems of differentialequations . The requirements are that the coordinate functions ) be continuous with continuous derivatives with respect to x₁,x₂, ...,x_n, that the f_j are weakly concave (or weakly convex)in , and that a certain condition pertaining to strict concavity (or strict convexity) is fulfilled.2000 Mathematics Subject Classification 34C25, 34C12.     

Immanant Inequalities, Induced Characters, and Rank Two Partitions

If = {₁, ₂, ..., _s}, where _{1 2} ... _s > 0, is a partitionof n then denotes the associated irreducible character of S_n,the symmetric group on {1, 2, ..., n}, and, if cCS_n, the groupalgebra generated by C and S_n, then d_c(·) denotes thegeneralized matrix function associated with c. If c₁, c₂ CS_nthen we write c₁ c₂ in case (A) (A) for each n x n positivesemi-definite Hermitian matrix A. If cCS_n and c(e) 0, wheree denotes the identity in S_n, then or denotes (c(e))^–1 c. The main result, an estimate for the norms of tensors of a certainanti-symmetry type, implies that if = {₁, ₂, ..., _s, 1^t} isa partition of n such that s > 1 and _s = 2, and ' denotes{₁, ₂, ..., _s-1, 1^t} then (, _{2}) where denotes characterinduction from S_n–2 x S₂ to S_n. This in turn implies thatif = {₁, ₂, ..., _s, 1^t} with s > 1, _s = 2, and ßdenotes {₁ + 2, ₂, ..., _s-1, 1^t} then _ß which,in conjunction with other known results, provides many new inequalitiesamong immanants. In particular it implies that the permanentfunction dominates all normalized immanants whose associatedpartitions are of rank 2, a result which has proved elusivefor some years. We also consider the non-relationship problem for immanants– that is the problem of identifying pairs, (,ß)such that _ß and _ß are both false.     

The n-Fold Central Haagerup Tensor Product

This paper defines the n-fold central Haagerup tensor product of a von Neumann algebra R, and shows that the map given by ([a₁a₂...a_n])(x₁, x₂, ...,x_n–1)=a₁x₁a₂x₂...x_n–1a_n is an isometry.