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.     

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     

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.     

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.     

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}.     

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.     

Local error estimates for radial basis function interpolation of scattered data

Introducing a suitable variational formulation for the localerror of scattered data interpolation by radial basis functions(r), the error can be bounded by a term depending on the Fouriertransform of the interpolated function f and a certain ‘Krigingfunction’, which allows a formulation as an integral involvingthe Fourier transform of . The explicit construction of locallywell-behaving admissible coefficient vectors makes the Krigingfunction bounded by some power of the local density h of datapoints. This leads to error estimates for interpolation of functionsf whose Fourier transform f is ‘dominated’ by thenonnegative Fourier transform of (x) = (||x||) in the sense . Approximation orders are arbitrarily high for interpolationwith Hardy multiquadrics, inverse multiquadrics and Gaussiankernels. This was also proven in recent papers by Madych andNelson, using a reproducing kernel Hilbert space approach andrequiring the same hypothesis as above on f, which limits thepractical applicability of the results. This work uses a differentand simpler analytic technique and allows to handle the casesof interpolation with (r) = r^s for s R, s > 1, s 2N, and(r) = r^s log r for s 2N, which are shown to have accuracy O(h^s/2)     

Maximum Precision in Elliott-Donaldson Formulae

This paper considers the approximation of by quadratures which are constructed inthe manner described by Elliott & Donaldson. Results areestablished for choosing the nodes so that the quadrature is exact on as wide a range of positive and negativepowers of z as possible.     

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.     

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.     

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.     

Stability Analysis of Runge-Kutta Methods for Volterra Integral Equations of the Second Kind

Stability analysis of Volterra-Runge-Kutta methods based onthe basic test equation of the form where is a complex parameter, and on the convolution test equation where and are real parameters, is presented. General stabilityconditions are derived and applied to construct numerical methodswith good stability properties. In particular, a family of second-orderV^o-stable Volterra-Runge-Kutta methods is obtained. No V^o-stablemethods of order greater than one have been presented previouslyin the literature.     

Some Applications of the Boundary-Locus Method and the Method of D-Partitions

We summarize aspects of the method of D-partitions and its analogue,the boundary-locus technique. We apply these techniques to adiscussion of the stability of the delay equation where y(t)=(t), t [–, 0]. Some comparisons with resultsvalid for other functional equations are introduced. Of interestare the following features: (i) the results corresponding toour results but applicable to delay-differential equations arenot found in the literature, (ii) some unconventional parameterspaces are introduced, and (iii) the possibility of extendingresults which hold for a basic test equation to more generalequations is established.     

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.     

Periodic Solutions of Lienard Equations with Asymmetric Nonlinearities at Resonance

The existence of 2-periodic solutions of the second-order differentialequation where a, b satisfy and p(t)=p(t+2),t R, is examined. Assume that limits lim_x±F(x)=F(±)(F(x)=) and lim_x±g(x)=g(±)exist and are finite. It is proved that the equation has atleast one 2-periodic solution provided that the zeros of thefunction ₁ are simple and the zeros of the functions ₁, ₂ aredifferent and the signs of ₂ at the zeros of ₁ in [0,2/n) donot change or change more than two times, where ₁ and ₂ aredefined as follows: Moreover, it is also proved that the given equation has at leastone 2-periodic solution provided that the following conditionshold: with 1 p < q 2.     

Numerical Evaluation of a Class of Double Integrals of Oscillatory Functions

A general method is described for the numerical evaluation ofintegrals of the form where, in a typical case, (x, y) is a lengthy polynomial ofdegree 10 in (x, y) and A is the region common to three overlappingcircles, although the method is in no way restricted to suchcases. Illustrative numerical examples are given.     

Convergence of Univariate Quasi-Interpolation Using Multiquadrics

Quasi-interpolants to a function f: RR on an infinite regularmesh of spacing h can be defined by where :RR is a function with fast decay for large argument. In the approach employing the radial-basis-function : RR, thefunction is a finite linear combination of basis functions(|•–jh|) (jZ). Choosing Hardy's multiquadrics (r)=(r²+c²)^?,we show that sufficiently fast-decaying exist that render quasi-interpolationexact for linear polynomials f. Then, approximating f C²(R),we obtain uniform convergence of s to f as (h, c)0, and convergenceof s' to f' as (h, c²/h)0. However, when c stays bounded awayfrom 0 as h0, there are f C(R) for which s does not convergeto f as h0. We also show that, for all which vanish at infinity but arenot integrable over R, there are no finite linear combinations of the given basis functions allowing the construction of admissiblequasi-interpolants. This includes the case of the inverse multiquadncs(r)=(r²+c²)^–?.     

Convergence of a finite-element approximation of surfactant spreading on a thin film in the presence of van der Waals forces

We consider a fully practical finite-element approximationof the following system of nonlinear degenerate parabolic equations: (u)/(t) + . (u² [(v)]) - (1)/(3) .(u³ w)= 0, w = - c u - u^-+ a u^-3 , (v)/(t) + . (u v [(v)]) - v - .(u² v w) = 0. The above models a surfactant-driven thin-film flow in the presenceof both attractive, a>0, and repulsive, >0 with >3,van der Waals forces; where u is the height of the film, v isthe concentration of the insoluble surfactant monolayer and(v):=1-v is the typical surface tension. Here 0 and c>0 arethe inverses of the surface Peclet number and the modified capillarynumber. In addition to showing stability bounds for our approximation,we prove convergence, and hence existence of a solution to thisnonlinear degenerate parabolic system, (i) in one space dimensionwhen >0; and, moreover, (ii) in two space dimensions if inaddition 7. Furthermore, iterative schemes for solving the resultingnonlinear discrete system are discussed. Finally, some numericalexperiments are presented.     

Oscillation Criteria for First-Order Delay Equations

This paper is concerned with the oscillatory behaviour of first-orderdelay differential equations of the form (1) where is non-decreasing, (t)< t for t t₀ and . Let the numbers k andL be defined by It is proved here that when L < 1 and 0 < k 1/e all solutionsof equation (1) oscillate in several cases in which the condition holds, where ₁ is the smaller root of the equation = e^k. 2000Mathematics Subject Classification 34K11 (primary); 34C10 (secondary).     

Spectral Multipliers for Hardy Spaces Associated with Schrodinger Operators with Polynomial Potentials

Let T_t be the semigroup of linear operators generated by a Schrödingeroperator – A = – V, where V is a non-negative polynomial,and let be the spectral resolution of A. We say that f is an element of if the maximal function Mf(x) = sup_t>0|T_tf(x)| belongs toL_p. We prove a criterion of Mihlin type on functions F whichimplies boundedness of the operators on , 0 < p 1. 1991 MathematicsSubject Classification 42B30, 35J10.