The Lambert W-functions and some of their integrals: a case study of high-precision computation

The real-valued Lambert W-functions considered here are w ₀(y) and w _− 1(y), solutions of we ^w  = y, − 1/e < y < 0, with values respectively in ( − 1,0) and ( − ∞ , − 1). A study is made of the numerical evaluation to high precision of these functions and of the integrals ò₁^￥ [-w₀(-xe^-x)]^a x^-bx\int_1^\infty [-w_0(-xe^{-x})]^\alpha x^{-\beta}\d x, α > 0, β ∈ ℝ, and ò₀¹ [-w_-1(-x e^-x)]^a x^-bx\int_0^1 [-w_{-1}(-x e^{-x})]^\alpha x^{-\beta}\d x, α > − 1, β < 1. For the latter we use known integral representations and their evaluation by nonstandard Gaussian quadrature, if α ≠ β, and explicit formulae involving the trigamma function, if α = β.     

A note on a Liouville type result of Gilbarg and Weinberger for the stationary Navier–Stokes equations in 2<Emphasis Type="Italic">D</Emphasis>

We show that a velocity field u satisfying the stationary Navier–Stokes equations on the entire plane must be constant under the growth condition lim sup |x|^−α |u(x)| < ∞ as |x| → ∞ for some α ∈ [0, 1/7).^† Bibliography: 10 titles.     

Fractional Poisson equations and ergodic theorems for fractional coboundaries

For a given contractionT in a Banach spaceX and 0<α<1, we define the contractionT _α=Σ _j=1 ^∞ a _j T ^j , where {a _j } are the coefficients in the power series expansion (1-t)^α=1-Σ _j=1 ^∞ a _j t ^j in the open unit disk, which satisfya _j >0 anda _j >0 and Σ _j=1 ^∞ a _j =1. The operator calculus justifies the notation(I−T) ^α :=I−T _α (e.g., (I−T _1/2)²=I−T). A vectory∈X is called an, α-fractional coboundary for T if there is anx∈X such that(I−T) ^α x=y, i.e.,y is a coboundary forT _α . The fractional Poisson equation forT is the Poisson equation forT _α . We show that if(I−T)X is not closed, then(I−T) ^α X strictly contains(I−T)X (but has the same closure). ForT mean ergodic, we obtain a series solution (converging in norm) to the fractional Poisson equation. We prove thaty∈X is an α-fractional coboundary if and only if Σ _k=1 ^∞ T ^k y/k ^1-α converges in norm, and conclude that lim _n ‖(1/n ^1-α)Σ _k=1 ⁿ T ^k y‖=0 for suchy. For a Dunford-Schwartz operatorT onL ₁ of a probability space, we consider also a.e. convergence. We prove that iff∈(I−T) ^α L ₁ for some 0<α<1, then the one-sided Hilbert transform Σ _k=1 ^∞ T ^k f/k converges a.e. For 1<p<∞, we prove that iff∈(I−T) ^α L _p with α>1−1/p=1/q, then Σ _k=1 ^∞ T ^k f/k ^1/p converges a.e., and thus (1/n ^1/p ) Σ _k=1 ⁿ T ^k f converges a.e. to zero. Whenf∈(I−T) ^1/q L _p (the case α=1/q), we prove that (1/n ^1/p (logn)^1/q )Σ _k=1 ⁿ T ^k f converges a.e. to zero.     

Piecewise-linear,discontinuous Galerkin method for a fractional diffusion equation

We use a piecewise-linear, discontinuous Galerkin method for the time discretization of a fractional diffusion equation involving a parameter in the range − 1 < α < 0. Our analysis shows that, for a time interval (0,T) and a spatial domain Ω, the error in L_￥((0,T);L₂(W))L_\infty\bigr((0,T);L_2(\Omega)\bigr) is of order k ^2 + α , where k denotes the maximum time step. Since derivatives of the solution may be singular at t = 0, our result requires the use of non-uniform time steps. In the limiting case α = 0 we recover the known O(k ²) convergence for the classical diffusion (heat) equation. We also consider a fully-discrete scheme that employs standard (continuous) piecewise-linear finite elements in space, and show that the additional error is of order h ²log(1/k). Numerical experiments indicate that our O(k ^2 + α ) error bound is pessimistic. In practice, we observe O(k ²) convergence even for α close to − 1.     

Harnack Inequalities for some Lévy Processes

In this paper we prove Harnack inequality for nonnegative functions which are harmonic with respect to random walks in ℝ ^d . We give several examples when the scale invariant Harnack inequality does not hold. For any α ∈ (0,2) we also prove the Harnack inequality for nonnegative harmonic functions with respect to a symmetric Lévy process in ℝ ^d with a Lévy density given by $c|x|^{-d-\alpha}1_{\{|x|\leq 1\}}+j(|x|)1_{\{|x|>1\}}$c|x|^{-d-\alpha}1_{\{|x|\leq 1\}}+j(|x|)1_{\{|x|>1\}}, where 0 ≤ j(r) ≤ cr ^{− d − α} , ∀ r > 1, for some constant c. Finally, we establish the Harnack inequality for nonnegative harmonic functions with respect to a subordinate Brownian motion with subordinator with Laplace exponent ϕ(λ) = λ ^α/2ℓ(λ), λ > 0, where ℓ is a slowly varying function at infinity and α ∈ (0,2).     

Controllability in a filled domain for the multidimensional wave equation with a singular boundary control

In accordance with the demands of the so-called local approach to inverse problems, the set of “waves” u^f (·, T) is studied, where u^f (x,t) is the solution of the initial boundary-value problem u_tt−Δu=0 in Ω×(0,T), u|_t<0=0, u|_∂Ω×(0,T)=f, and the (singular) control f runs over the class L²((0,T); H^−m (∂Ω)) (m>0). The following result is established. Let Ω^T={x ∈ Ω : dist(x, ∂Ω)<T)} be a subdomain of Ω ⊂ ℝⁿ (diam Ω<∞) filled with waves by a final instant of time t=T, let T_*=inf{T : Ω^T=Ω} be the time of filling the whole domain Ω. We introduce the notation D_m=Dom((−Δ)^m/2), where (−Δ) is the Laplace operator, Dom(−Δ)=H²(Ω)∩H ₀ ¹ (Ω);D_−m=(D_m)′;D_−m(Ω^T)={y∈D_−m:supp y ⋐ Ω^T. If T<T., then the reachable set R _m ^T ={u^t(·, T): f ∈ L₂((0,T), H^−m (∂Ω))} (∀m>0), which is dense in D_−m(Ω^T), does not contain the class C ₀ ^∞ (Ω^T). Examples of a ∈ C ₀ ^∞ , a ∈ R _m ^T , are presented. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 210, 1994, pp. 7–21. Translated by T. N. Surkova.     

Failure rates of regenerative systems with heavy tails

Let ξ_t be a regenerative process and assume that, at each state x, the process can fail with intensity α(x). If the inter-regeneration times have a finite exponential moment orinf _xα(x)>0, then α(ξ_t) tends to some limiting positive intensity as t→∞ (under mild additional restrictions). This fact is widely used in engineering because the limiting intensity can be employed in various calculations, say, in reliability theory. The paper contains a variety of examples showing that α(ξ_t) provided that inter-regeneration time has no exponential moment andinf _x α(x)=0. The speed of convergence depends, in general, on both the tail of inter-regeneration time and α(x). Proceedings of the Seminar on Stability Problems for Stochastic Models, Hajdúszoboszló, Hungary, 1997, Part III.     

Global asymptotics for Laguerre polynomials with large negative parameter—a Riemann-Hilbert approach

In this paper, we study the asymptotic behavior of the Laguerre polynomials as n→∞. Here α _n is a sequence of negative numbers and −α _n /n tends to a limit A>1 as n→∞. An asymptotic expansion is obtained, which is uniformly valid in the upper half plane ℂ₊={z:Im z≥0}. A corresponding expansion is also given for the lower half plane ℂ₋={z:Im z≤0}. The two expansions hold, in particular, in regions containing the curve Γ in the complex plane, on which these polynomials are orthogonal. Our method is based on the Riemann-Hilbert approach introduced by Deift and Zhou. The work of R. Wong is partially supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 102504).     

The finite-dimensionalP λ spaces with small λ

It is proved that there is a positive function Φ(∈) defined for sufficiently small ∈>0 such that lim_∈→0 Φ(∈)=0 and for every integerk and everyk-dimensionalP _1+∈ spaceE, d(E, l _∞ ^k )<1+Φ(∈). Author was partially supported by N.S.F. Grant MCS 79-03042. An erratum to this article is available at .     

On meromorphic univalent functions omitting a disc

The class Σ_b is defined to consist of meromorphic univalent functionsH omitting a disc with the radiusb:H(z)=z+ Σ ₀ ^∞ A _n z ⁻ⁿ ,z>1,H(b)>b ∈ (0, 1). By aid of FitzGerald inequalities the inverse coefficients of odd Σ_b-functions are maximized. The result extends the corresponding estimation, due to Netanyahu and Schober, fromb=0 to the whole interval (0, 1). The author wishes to express her gratitude to Professor O. Tammi for valuable discussions connected with the problem. This work was supported by a grant from the Finnish Ministry of Education.     

Asymptotic behavior of solutions for a class of retarded difference equations

Consider the retarded difference equationx _n −x _n−1 =F(−f(x _n )+g(x _n−k )), wherek is a positive integer,F,f,g:R→R are continuous,F andf are increasing onR, anduF(u)>0 for allu≠0. We show that whenf(y)≥g(y) (resp. f(y)≤g(y)) fory∈R, every solution of (*) tends to either a constant or −∞ (resp. ∞) asn→∞. Furthermore, iff(y)≡g(y) fory∈R, then every solution of (*) tends to a constant asn→∞. Project supported by NNSF (19601016) of China and NSF (97-37-42) of Hunan     

Some Oscillation Theorems for Second Order Differential Equations

In this paper we establish some oscillation or nonoscillation criteria for the second order half-linear differential equation

Average B-width and infinite dimensional G-width of some smooth function classes on the line

In this paper, we get the exact values of average σ-B width and infinite dimensional σ-G width of Sobolev class B^r _p(R) in the metric L_p(R) (1≤p≤∞) and obtain the exact (σ∈N) and strong asymptotic (σ>1) results of infinite dimensional σ-G widths of Sobolev-Wiener class W^r _pq (R) in the metric L_q(R) and its dual case W^r _p(R) in the metric L_qp(R) (1≤q≤p≤∞). Supported by the National Natural Science Foundation of China (No. 19671012) and by the Doctoral Programme Foundation of Institution of Higher Education of National Education Committee of China.     

Some characterizations of weighted Herz-type Hardy spaces and their applications

In this paper, the authors first establish some new real-variable characterizations of Herz-type Hardy spaces and , where ω₁,ω₃ ∈ A₁-weight, 1<q>∞,n(1−1/q)≤α<∞ and 0<p<∞. Then, using these new characterizations, they investigate the convergence of a bounded set in these spaces, and study the boundedness of some potential operators on these spaces. Supported by the NNSF of China     

<Emphasis Type="Italic">L</Emphasis><Superscript>∞</Superscript> Estimates of solution for <Emphasis Type="Italic">m</Emphasis>-Laplacian parabolic equation with a nonlocal term

In this paper, we consider the global existence, uniqueness and L ^∞ estimates of weak solutions to quasilinear parabolic equation of m-Laplacian type u _t − div(|∇u| ^m−2∇u) = u|u| ^β−1 ∫_Ω |u| ^α dx in Ω × (0,∞) with zero Dirichlet boundary condition in tdΩ. Further, we obtain the L ^∞ estimate of the solution u(t) and ∇u(t) for t > 0 with the initial data u ₀ ∈ L ^q (Ω) (q > 1), and the case α + β < m − 1.     

On the degree of continuity and smoothness of sine and cosine Fourier transforms of Lebesgue integrable functions

We consider complex-valued functions f∈L ¹(ℝ₊), where ℝ₊:=[0,∞), and prove sufficient conditions under which the sine Fourier transform [^(f)]_s\hat{f}_{s} and the cosine Fourier transform [^(f)]_c\hat{f}_{c} belong to one of the Lipschitz classes Lip (α) and lip (α) for some 0<α≦1, or to one of the Zygmund classes Zyg (α) and zyg (α) for some 0<α≦2. These sufficient conditions are best possible in the sense that they are also necessary if f(x)≧0 almost everywhere.     

On the law of the iterated logarithm for subsequences for a stable subordinator

Let {X(t): 0≤t<∞) be a stable subordinator with α∈(0,1). For increasing sequences t_k we give normalizing constants a_k such thatliminf _k→∞ a _k ⁻¹ X(t_k) is a.s constant. We also derive a.s. upper bounds. Proceedings of the XVI Seminar on Stability Problems for Stochastic Models Part I, Eger, Hungary, 1994.     

On the order of growth of solutions of algebraic differential equations

Assume thatf is an integer transcendental solution of the differential equationP _n (z, f, f′)=P _n−1(z, f, f′, ... f ^(p)), whereP _n andP _n−1 are polynomials in all variables, the degree ofP _n with respect tof andf′ is equal ton, and the degree ofP _n−1 with respect tof, f′, ... f ^(p) is at mostn−1. We prove that the order ρ of growth off satisfies the relation 1/2≤ρ<∞. We also prove that if ρ=1/2, then, for a certain real ν, in the domain {z: ν<argz<ν+2π}/E _*, whereE _* is a certain set of disks with finite sum of radii, the estimate lnf(z)=z ^1/2 (β+o(1)), β∈C, holds forz=re ^iϕ,r≥r(ϕ)≥0. Furthermore, on the ray {z: argz=ν}, the following relation is true: ln‖f(re ^iν)‖=o(r ^1/2),r→+∞,r>0, , where Δ is a certain set on the semiaxisr>0 with mes Δ<∞. “L'vivs'ka Politekhnika” University, Lvov. Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 51, No. 1, pp. 69–77, January, 1999.     

The Besov subspace consisting of most non-smooth functions

Under the assumptions that Δ(f, h)(t) = |f(t + h) − f(t)|, X is a symmetric space of functions in [0, 1], α ∈ (0, 1) and p ∈ [1, ∞) are any fixed number, by the triple (X, α, p) a Besov type space Λ _X,p ^α is constructed, where the norm is given by the equality

Spectral invariance for pseudodifferential operators on weighted function spaces

The spectrum of each symmetric ψ DO of the symbol class S⁰ _{1, γ}, 0≤γ<1, acting on B³ _p,q(w(x)) and F³ _p,q(w(x)), is independent of the choice ofs ∈ℝ, 0<p≤∞ (p<∞ in the F-case), 0<q≤∞ and the weight w(x)∈W.