A new form of trigonometric orthogonality and Gaussian-type quadrature

Gauss's (2n+1)-point trigonometric interpolation formula, based upon f(x_i), i = 1(1)2n+1, gives a trigonometric sum of the n^th order, S_2n+1(x = a₀ + ∑_jⁿ = 1^(a_jcos jx + b_jsin jx), which may be integrated to provide formulas for either direct quadrature or stepwise integration of differential equations having periodic (or near-periodic) solutions. An “orthogonal” trigonometric sum S_2r+1(x) is one that satisfies

$\text{∫}{}_{\mathrm{a}}{}^{\mathrm{b}}\text{S}{}_{\mathrm{2r+1}}\text{(x)S}{}_{\mathrm{2r\prime +1}}\text{(x)dx=0, r′<r}$

and two other arbitrarily imposable conditions needed to make S_2r1(x) unique. Two proofs are given of a fundamental factor theorem for any S_2n+1(x) (somewhat different from that for polynomials) from which we derive 2r-point Gaussian-type quadrature formulas, r = [n/2] + 1, which are exact for any S_4r?1(x). We have

$\text{∫}{}_{\mathrm{a}}{}^{\mathrm{b}}\text{S}{}_{\mathrm{4r?1(x)}}\text{dx=∑}{}_{\mathrm{j=1}}{}^{\mathrm{2r}}\text{A}{}_{\mathrm{j}}\text{S}{}_{\mathrm{4r?1}}\text{(x}{}_{\mathrm{j}}\text{)}$

where the nodes x_j, j = 1(1)2r, are the zeros of the orthogonal S_2r+1(x). It is proven that A_j > 0 and that 2r-1 of the nodes must lie within the interval [a,b], and the remaining node (which may or may not be in [a,b]) must be real. Unlike Legendre polynomials, any [a′,b′] other than a translation of [a,b], requires different and unrelated sets of nodes and weights. Gaussian-type quadrature formulas are applicable to the numerical integration of the Gauss (2n+1)-point interpolation formulas, with extra efficiency when the latter are expressed in barycentric form. S_2r+1(x), x_jandA_j, j = 1(1)2r, were calculated for [a,b] = [0, π/4], 2r = 2 and 4, to single-precision accuracy.     

Phase reconstruction from magnitude of band-limited multidimensional signals

In this paper, the problem of phase reconstruction from magnitude of multidimensional band-limited functions is considered. It is shown that any irreducible band-limited function f(z₁…,z_n), z_i ? $\text{C}$, i=1, …, n, is uniquely determined from the magnitude of f(x₁…,x_n): | f(x₁…,x_n)|, x_i ? $\text{R}$, i=1,…, n, except for (1) linear shifts: ^{i(α1z1+…+α_n2n+β}), β, α_i?$\text{R}$, i=1,…, n; and (2) conjugation: $\text{f}{}^1\text{(z}{}_1{}^1\text{,…,z}{}_{\mathrm{n}}{}^1\text{)}$.     

Temporally inhomogeneous scattering theory

A theory of scattering for the time dependent evolution equations $\text{du}\text{dt}\text{= iH}{}_{\mathrm{j}}\text{(t)u, j = 0, 1}$ (¹) is developed. The wave operators are defined in terms of the evolution operators U_j(t, s), which govern (¹). The scattering operator remains unitary. Sufficient conditions for existence and completeness of the wave operators are obtained; these are the main results. General properties, such as the chain rule and various intertwining relations, are also established. Applications include potential scattering (H₀(t) = ?Δ, Δ denoting the Laplacian, and H₁(t) = ?Δ + q(t, ·)) and scattering for second-order differential operators with coefficients constant in the spatial variable ($\text{H}{}_{\mathrm{j}}\text{(t) = ∑}{}_{\mathrm{m,\; k\; =\; 1}}{}^{\mathrm{n}}\text{a}{}_{\mathrm{mk}}{}^{\mathrm{(j)}}\text{(t)(}\text{?}{}^2\text{?x}{}_{\mathrm{m}}\text{?x}{}_{\mathrm{k}}\text{) + b}{}_{\mathrm{j}}\text{(t)}\text{for}\text{j = 0, 1}$).     

The Dirichlet problem for the Kohn Laplacian on the Heisenberg group,II

Let L = ∑_{j = 1}^mX_j² be sum of squares of vector fields in $\text{R}$ⁿ satisfying a Hörmander condition of order 2: span{X_j, [X_i, X_j]} is the full tangent space at each point. A point x??D of a smooth domain D is characteristic if X₁,…, X_m are all tangent to ?D at x. We prove sharp estimates in non-isotropic Lipschitz classes for the Dirichlet problem near (generic) isolated characteristic points in two special cases: (a) The Grushin operator $\text{?}{}^2\text{?x}{}^2\text{+ x}{}^2\text{?}{}^2\text{?t}{}^2$ in $\text{R}$². (b) The real part of the Kohn Laplacian on the Heisenberg group $\text{∑}{}_{\mathrm{j\; ?\; 1}}{}^{\mathrm{n}}\text{(}\text{?}\text{?x}{}_{\mathrm{j}}\text{+ 2y}{}_{\mathrm{j}}\text{?}\text{?t}\text{)}{}^2\text{+ (}\text{?}\text{?y}{}_{\mathrm{j}}\text{? 2x}{}_{\mathrm{j}}\text{?}\text{?t}\text{)}{}^2$ in $\text{R}$^{2n + 1}. In contrast to non-characteristic points, C^∞ regularity may fail at a characteristic point. The precise order of regularity depends on the shape of ?D at x.     

A method of feasible directions using function approximations,with applications to min max problems

This paper presents a demonstrably convergent method of feasible directions for solving the problem min{φ(ξ)| gⁱ(ξ)?0i=1,2,…,m}, which approximates, adaptively, both φ(x) and ▽φ(x). These approximations are necessitated by the fact that in certain problems, such as when $\text{φ(x) =}\text{max}\text{{f(x, y) ¦ y ? Ω}{}_{\mathrm{y}}\text{}}$, a precise evaluation of φ(x) and ▽φ(x) is extremely costly. The adaptive procedure progressively refines the precision of the approximations as an optimum is approached and as a result should be much more efficient than fixed precision algorithms.It is outlined how this new algorithm can be used for solving problems of the form $\text{min}{}_{\mathrm{y}}\text{? Ω}{}_{\mathrm{x}}\text{max}{}_{\mathrm{y}}\text{? Ω}{}_{\mathrm{y}}\text{f(x, y)}$ under the assumption that Ωm_ξ={x|gⁱ(x)?0, j=1,…,s} ∩$\text{R}$ⁿ, Ω_y={y|ζ_i(y)?0, i-1,…,t} ∩ $\text{R}$^m, with f, g^j, ζⁱ continuously differentiable, f(x, ·) concave, ζⁱ convex for $\text{i = 1,…, t,}\text{and}\text{Ω}{}_{\mathrm{x}}\text{, Ω}{}_{\mathrm{y}}$ compact.     

Elementary proof of the transformation formula for Lambert series involving generalized Dedekind sums

An elementary proof is given of the author's transformation formula for the Lambert series $\text{G}{}_{\mathrm{p}}\text{(x) =}\text{Σ}{}_{\mathrm{n?1}}{}^{\mathrm{<mtext>\infty </mtext>}}\text{n}{}^{\mathrm{?p}}\text{x}{}^{\mathrm{n}}\text{(1?x}{}^{\mathrm{n}}\text{)}$ relating G_p(e^2πiτ) to G_p(e^2πiAτ), where p > 1 is an odd integer and $\text{Aτ =}\text{(aτ + b)}\text{(cτ + d)}$ is a general modular substitution. The method extends Sczech's argument for treating Dedekind's function $\text{log}\text{η(τ) =}\text{πiτ}\text{12}\text{? G}{}_1\text{(e}{}^{\mathrm{2\pi i\tau }}\text{)}$, and uses Carlitz's formula expressing generalized Dedekind sums in terms of Eulerian functions.     

Arithmetic progressions that consist only of primes

Let N_m(x) be the number of arithmetic progressions that consist of m terms, all primes and not larger than x, and set $\text{F}{}_{\mathrm{m}}\text{(x) =}\text{C}{}_{\mathrm{m}}\text{x}{}^2\text{log}{}^{\mathrm{m}}\text{x}$ (C_m explicitly given). It is shown that Hardy and Littlewood's prime k-tuple conjecture implies that N_m(x) = F_m(x){1 + Σ_j=1^Na_jlog^?jx + O((log x)^?N?1)}, (here the bracket represents an asymptotic series with explicitly computable coefficients). This formula holds rather trivially for m = 1 and m = 2. It is proved here for m = 3, by the Vinogradov version of the Hardy-Ramanujan-Littlewood circle method.     

On the rate for uniform strong consistency of empirical distributions of independent nonidentically distributed multivariate random variables

Let X_j = (X_1j ,…, X_pj), j = 1,…, n be n independent random vectors. For x = (x₁ ,…, x_p) in R^p and for α in [0, 1], let F_j${}^1$(x) = αI(X_1j < x₁ ,…, X_pj < x_p) + (1 ? α) I(X_1j ≤ x₁ ,…, X_pj ≤ x_p), where I(A) is the indicator random variable of the event A. Let F_j(x) = E(F_j${}^1$(x)) and D_n = sup_{x, α} max_{1 ≤ N ≤ n} |Σ₀ⁿ(F_j${}^1$(x) ? F_j(x))|. It is shown that P[D_n ≥ L] < 4pL exp{?2(L²n^?1 ? 1)} for each positive integer n and for all L² ≥ n; and, as n → ∞, $\text{D}{}_{\mathrm{n}}\text{= 0((n}\text{log}\text{n)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}$ with probability one.     

Orthogonal polynomials on the negative multinomial distribution

Orthogonal polynomials on the multivariate negative binomial distribution,

$\text{(1 + Θ)}{}^{\mathrm{?\alpha ?x}}\text{(}\text{π}\text{j=0}\text{p}\text{Θ}{}_{\mathrm{j}}{}^{\mathrm{x}}\text{j}\text{x}{}_{\mathrm{j}}\text{!}\text{)}\text{Γ(α + x)}\text{Γ(α)}$

where α > 0, Θ₁ > 0, x = ΣΘ_i, x₀, x₁, …, x_p = 0,1, … are constructed and their properties studied.     

Infinite-variate wide-sense Markov processes and functional analysis for bounded operator-forming vectors

Let p, q be arbitrary parameter sets, and let $\text{H}$ be a Hilbert space. We say that x = (xⁱ)_i?q, xⁱ ? $\text{H}$, is a bounded operator-forming vector (?$\text{H}$_F^q) if the Gram matrix 〈x, x〉 = [(xⁱ, x^j)]_i?q,j?q is the matrix of a bounded (necessarily ≥ 0) operator on $\text{l}{}_{\mathrm{q}}{}^2$, the Hilbert space of square-summable complex-valued functions on q. Let A be p × q, i.e., let A be a linear operator from $\text{l}{}_{\mathrm{q}}{}^2$ to $\text{l}{}_{\mathrm{p}}{}^2$. Then exists a linear operator ǎ from (the Banach space) $\text{H}$_F^q to $\text{H}$_F^p on $\text{D}$(A) = {x:x ? $\text{H}$_F^q, $\text{A〈x, x〉}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ is p × q bounded on $\text{l}{}_{\mathrm{q}}{}^2$} such that y = ǎx satisfies y^j?σ(x) = {space spanned by the xⁱ}, 〈y, x〉 = A〈x, x〉 and $\text{〈y, y〉 = A〈x, x〉}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{(A〈x, x〉}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}{}^1$. This is a generalization of our earlier [J. Multivariate Anal.4 (1974), 166–209; 6 (1976), 538–571] results for the case of a spectral measure concentrated on one point. We apply these tools to investigate q-variate wide-sense Markov processes.     

Nonhomogeneous dirichlet problem for elliptic operators with axially discontinuous coefficients

Si studia, in un cilindro, il problema di Dirichlet per l'equazione ellittica del II ordine: $\text{L}{}_{\mathrm{\alpha }}\text{u = ?}$, dove $\text{L}{}_{\mathrm{\alpha }}\text{=}\text{αΔ + (1 ? 3α)∑}{}_{\mathrm{ij\; =\; 1}}{}^2\text{x}{}_{\mathrm{i}}\text{x}{}_{\mathrm{j}}\text{(x}{}_1{}^2\text{+ x}{}_2{}^2\text{)}{}^{\mathrm{?1}}\text{?}{}^2\text{?x}{}_{\mathrm{i}}\text{?x}{}_{\mathrm{j}}\text{, α ? (0,}\text{1}\text{3}\text{]}$è l'operatore a coefficienti discontinui sull'asse x₃ già introdotto da N. Ural'tseva per mostrare che l'equazione considerata può non avere soluzione nello spazio di Sobolev W^2,p(p > 2) per qualche f?L^p. In questo lavoro si danno limitazioni a priori e teoremi di esistenza e unicità in W^2,p quando p varia in un intervallo (p₁(α), p₂(α)), dipendente dalla costante di ellitticità α. Se p = p₂(α) le limitazioni a priori cadono: l'esempio è quello di Ural'tseva.     

Numerical approximation of product integrals

A technique for the numerical approximation of matrix-valued Riemann product integrals is developed. For a ? x < y ? b, I_m(x, y) denotes

$\text{∫}\text{χ}\text{y}\text{∫}\text{χ}\text{v}{}_2\text{?}\text{∫}\text{χ}\text{v}{}_2\text{∏}\text{i=1}\text{m}\text{F(ν}{}_{\mathrm{i}}\text{)dν}{}_1\text{dν}{}_2\text{?dν}{}_{\mathrm{m}}$

, and A_m(x, y) denotes an approximation of I_m(x, y) of the form

$\text{(y?x)}{}^{\mathrm{m}}\text{∑}\text{k=1}\text{n}\text{a}{}_{\mathrm{k}}\text{∏}\text{i=1}\text{m}\text{F(χ}{}_{\mathrm{ik}}\text{)}$

, where a_k and y_ik are fixed numbers for i = 1, 2,…, m and k = 1, 2,…, N and x_ik = x + (y ? x)y_ik. The following result is established. If p is a positive integer, F is a function from the real numbers to the set of w × w matrices with real elements and F⁽¹⁾ exists and is continuous on [a, b], then there exists a bounded interval function H such that, if n, r, and s are positive integers, $\text{(b ? a)}\text{n}\text{= h < 1, x}{}_{\mathrm{i}}\text{= a + hi}\text{for}\text{i = 0, 1,…, n}\text{and}\text{0 < r ? s ? n}$, then

$\text{∏}\text{χ}{}_{\mathrm{r?}}\text{χ}{}_{\mathrm{s}}\text{(I+F dχ)?}\text{∏}\text{i=r}\text{s}\text{I+}\text{∑}\text{j=1}\text{p}\text{I}{}_{\mathrm{j}}\text{(χ}{}_{\mathrm{i?1}}\text{,χ}{}_{\mathrm{i}}\text{)}$

$\text{=h}{}^{\mathrm{p}}\text{H(χ}{}_{\mathrm{r?1}}\text{,χ}{}_{\mathrm{s}}\text{)+O(h}{}^{\mathrm{p+1}}\text{)}$

Further, if F^(j) exists and is continuous on [a, b] for j = 1, 2,…, p + 1 and A is exact for polynomials of degree less than p + 1 ? j for j = 1, 2,…, p, then the preceding result remains valid when A_j is substituted for I_j.     

Abstract Cauchy problems that involve a product of two Euler-Poisson-Darboux operators

Let X be a Banach space, let B be the generator of a continuous group in X, and let A = B². Assume that D(A^r) is dense in X for r an arbitrarily large positive integer and that a and b are non-negative reals. Solution representations are developed for the abstract differential equation

$\text{(D}{}^2{}_{\mathrm{t}}\text{+}\text{b}\text{t}\text{D}{}_{\mathrm{t}}\text{? A) · (D}{}^2{}_{\mathrm{t}}\text{+}\text{a}\text{t}\text{D}{}_{\mathrm{t}}\text{? A) u(t) = 0, t > 0}$

corresponding to initial conditions of the form: (i) u(0+) = φ, u^(j)(0+) = 0, j = 1, 2, 3 and (ii) u²(0+) = φ, u^j(0+) = 0, j = 0, 1, 3 (with φ∈D(A^r)) for all choices of a and b.     

Positive operators on the n-dimensional ice cream cone

Let $\text{K}{}_{\mathrm{n}}\text{= {x ? R}{}^{\mathrm{n}}\text{: (x}{}_1{}^2\text{+ · +x}{}^2{}_{\mathrm{n?1}}\text{)}\text{1}\text{2}\text{? x}{}_{\mathrm{n}}\text{}}$ be the n-dimensional ice cream cone, and let Γ(K_n) be the cone of all matrices in $\text{R}$ⁿⁿ mapping K_n into itself. We determine the structure of Γ(K_n), and in particular characterize the extreme matrices in Γ(K_n).     

An asymptotic evaluation of the cycle index of a symmetric group

Let $\text{Z(S}{}_{\mathrm{n}}\text{;?(x))}$ denote the polynomial obtained from the cycle index of the symmetric group Z(S_n) by replacing each variable s_i by f(x¹). Let f(x) have a Taylor series with radius of convergence ? of the form f(x)=x^k + a_k+1x^k+1 + a_k+2x^k+2+? with every a₁?0. Finally, let 0<x<1 and let x??. We prove that

This limit is used to estimate the probability (for n and p both large) that a point chosen at random from a random p-point tree has degree n + 1. These limiting probabilities are independent of p and decrease geometrically in n, contrasting with the labeled limiting probabilities of $\text{1}\text{n!e}$.In order to prove the main theorem, an appealing generalization of the principle of inclusion and exclusion is presented.     

Temporal exponential decay for the Stark effect in atoms

Let $\text{H}{}_1\text{= ?∑}{}_{\mathrm{i\; =\; 1}}{}^{\mathrm{N}}\text{(Δ}{}_{\mathrm{i}}\text{+ V(}\text{x}{}_{\mathrm{i}}\text{)) + ∑}{}_{\mathrm{1\; ?\; i\; <j\; ?\; N}}\text{¦}\text{x}{}_{\mathrm{i}}\text{?}\text{x}{}_{\mathrm{j}}\text{¦}{}^{\mathrm{?1}}\text{, V(}\text{x}{}_{\mathrm{i}}\text{) = N ∝ ¦}\text{x}\text{?}\text{y}\text{¦}{}^{\mathrm{?1}}\text{?(}\text{y}\text{)d}\text{y}$, with ? a normalized Gaussian. Suppose $\text{E}$ ≠ 0 and that $\text{H = H}{}_1\text{+}\text{E}\text{· (∑}{}_{\mathrm{i\; =\; 1}}{}^{\mathrm{N}}\text{x}{}_{\mathrm{i}}\text{)}$ has no eigenfunctions in L²($\text{R}$^3N. If H₁ψ = μψ with μ < infσ_ess(H₁), then (ψ, e^?itHψ) decays exponentially at a rate governed by the positions of the resonances of H.     

Isospectral matrices and classical polynomials

The matrices of order n defined, in terms of the n arbitrary numbers x_j, by the formulae $\text{X=}\text{diag}\text{(x}{}_{\mathrm{j}}\text{)}\text{and}\text{Z}{}_{\mathrm{jk}}\text{=δ}{}_{\mathrm{jk}}\text{∑′}\text{l=1}\text{n}\text{(x}{}_{\mathrm{j}}\text{?x}{}_{\mathrm{l}}\text{)}{}^{\mathrm{?1}}\text{+(1?δ}{}_{\mathrm{jk}}\text{(x}{}_{\mathrm{j}}\text{?x}{}_{\mathrm{k}}\text{)}{}^{\mathrm{?1}}$, are representations of the multiplicative operator ξ and of the differential operator d/dξ in a space spanned by the polynomials in ξ of degree less than n. This elementary fact implies a number of remarkable formulae involving these matrices, including novel representations of the classical polynomials.     

Two-parameter Turan inequalities for ultraspherical and Laguerre polynomials

It is known that the classical orthogonal polynomials satisfy inequalities of the form U_n²(x) ? U_{n + 1}(x) U_{n ? 1}(x) > 0 when x lies in the spectral interval. These are called Turan inequalities. In this paper we will prove a generalized Turan inequality for ultraspherical and Laguerre polynomials. Specifically if P_n^λ(x) and L_n^α(x) are the ultraspherical and Laguerre polynomials and $\text{F}{}_{\mathrm{n}}{}^{\mathrm{\lambda }}\text{(x) =}\text{P}{}_{\mathrm{n}}{}^{\mathrm{\lambda }}\text{(x)}\text{P}{}_{\mathrm{n}}{}^{\mathrm{\lambda }}\text{(1)}\text{, G}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x) =}\text{L}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x)}\text{L}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(0)}\text{,}\text{then}\text{F}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x) F}{}_{\mathrm{n}}{}^{\mathrm{\beta }}\text{(x) ? F}{}_{\mathrm{n\; +\; 1}}{}^{\mathrm{\alpha }}\text{(x) F}{}_{\mathrm{n\; ?\; 1}}{}^{\mathrm{\beta }}\text{(x) > 0, ? 1 < x < 1, ?}\text{1}\text{2}\text{< α ? β ? α + 1}\text{and}\text{G}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x) G}{}_{\mathrm{n}}{}^{\mathrm{\beta }}\text{(x) ? G}{}_{\mathrm{n\; +\; 1}}{}^{\mathrm{\alpha }}\text{(x) G}{}_{\mathrm{n\; ?\; 1}}{}^{\mathrm{\beta }}\text{(x) > 0, x > 0, 0 < α ? β ? α + 1}$. We also prove the inequality $\text{(n + 1) F}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x) F}{}_{\mathrm{n}}{}^{\mathrm{\beta }}\text{(x) ? nF}{}_{\mathrm{n\; +\; 1}}{}^{\mathrm{\alpha }}\text{(x) F}{}_{\mathrm{n\; ?\; 1}}{}^{\mathrm{\beta }}\text{(x) > A}{}_{\mathrm{n}}\text{[F}{}_{\mathrm{n}}{}^{\mathrm{\alpha }}\text{(x)]}{}^2\text{, ?1 < x < 1, ?}\text{1}\text{2}\text{< α ? β < α + 1,}\text{where}\text{A}{}_{\mathrm{n}}$ is a positive constant depending on α and β.