A note on polynomial matrix functions over a finite field

Let F_n denote the ring of n×n matrices over the finite field F=GF(q) and let A(x)=A_Nx^N+ ?+ A₁x+A₀?F_n[x]. A function $\text{?:F}{}_{\mathrm{n}}\text{→F}{}_{\mathrm{n}}$ is called a right polynomial function iff there exists an A(x)?F_n[x] such that $\text{?(B)=A}{}_{\mathrm{N}}\text{B}{}^{\mathrm{N}}\text{+?+A}{}_1\text{B+ A}{}_0$ for every B?F_n. This paper obtains unique representations for and determines the number of right polynomial functions.     

F-mixing and weak disjointness

A topological system (X,f) is $\text{F}$-transitive if for each pair of opene subsets U and V of X, $\text{N}{}_{\mathrm{f}}\text{(U,V)={n∈}\text{Z}{}_{\mathrm{+}}\text{:}\text{f}{}^{\mathrm{n}}\text{U∩V≠∅}∈}\text{F}$, where $\text{F}$ is a collection of subsets of $\text{Z}{}_{\mathrm{+}}$ which is hereditary upward. (X,f) is $\text{F}$-mixing if (X×X,f×f) is $\text{F}$-transitive. In this paper $\text{F}$-mixing systems are characterized in terms of the chaoticity of the systems. Moreover, weak disjointness is studied via family. We will give conditions such that a dual theorem of the Weiss–Akin–Glasner theorem holds. Examples with this dual theorem fails for some “good” families are obtained.     

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.     

Generalizations of two inequalities involving hermitian forms

Let λ₁ and λ_N be, respectively, the greatest and smallest eigenvalues of an N×N hermitian matrix H=(h_ij), and x=(x₁,x₂,…,x_N) with (x,x)=1. Then, it is known that (1) λ₁?(x,Hx)?λ_N and (2) if, in addition, H is positive definite, $\text{(λ}{}_1\text{+λ}{}_{\mathrm{N}}\text{)}{}^2\text{4λ}{}_1\text{λ}{}_{\mathrm{N}}\text{?(x,Hx)(x,H}{}^{\mathrm{?1}}\text{x)?1}$. Assuming that y=(y₁,y₂,…, y_N) and |y_i|?1, i=1,2,…,N, it is shown in this paper that these inequalities remain true if H and H^?1 are, respectively, replaced by the Hadamard products $\text{M(y)1H}$ and $\text{M(y)1H}{}^{\mathrm{?1}}$, where M(y) is a matrix defined by $\text{M(y)=(δ}{}_{\mathrm{ij}}\text{+(1?δ}{}_{\mathrm{ij}}\text{)}\text{y}{}_{\mathrm{i}}\text{y}{}_{\mathrm{j}}$. Subsequently, these results are extended to improve the spectral bounds of $\text{M(y)1H}$.     

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.     

A global genericity theorem for bifurcations in variational problems

Let C be a Banach space, H a Hilbert space, and let F(C,H) be the space of C^∞ functions f: C × H → R having Fredholm second derivative with respect to x at each (c, x) ?C × H for which $\text{D?}{}_{\mathrm{c}}\text{(x) = 0}$; here we write $\text{?}{}_{\mathrm{c}}\text{(x)}$ for $\text{?(c, x)}$. Say ? is of standard type if at all critical points of ?_c it is locally equivalent (as an unfolding) to a quadratic form Q plus an elementary catastrophe on the kernel of Q. It is proved that if f?F (A × B, H) satisfies a certain ‘general position’ condition, and dim B ? 5, then for most a?A the function f_o?F(B,H) is of standard type. Using this it is shown that those f?F(B,H) of standard type form an open dense set in F(B,H) with the Whitney topology. Thus both results are Hilbert-space versions of Thom's theorem for catastrophes in $\text{R}$ⁿ.     

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.     

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.     

Positive dependence properties of elliptically symmetric distributions

Let X₁, …, X_p have p.d.f. g(x₁² + … + x_p²). It is shown that (a) X₁, …, X_p are positively lower orthant dependent or positively upper orthant dependent if, and only if, X₁,…, X_p are i.i.d. N(0, σ²); and (b) the p.d.f. of |X₁|,…, |X_p| is TP₂ in pairs if, and only if, In g(u) is convex. Let X₁, X₂ have p.d.f. $\text{f(x}{}_1\text{, x}{}_2\text{) =}\text{|Σ|}{}^{\mathrm{<mtext>?</mtext><mtext>1</mtext><mtext>2</mtext>}}\text{g((x}{}_1\text{, x}{}_2\text{)}\text{Σ}{}^{\mathrm{?1}}\text{(x}{}_1\text{, x}{}_2\text{)′)}$. Necessary and sufficient conditions are given for f(x₁, x₂) to be TP₂ for fixed correlation ?. It is shown that if f is TP₂ for all ? >0. then (X₁, X₂)′ ～ N(0, Σ). Related positive dependence results and applications are also considered.     

On some properties of the Cesàro limit of a stochastic matrix

Let$\text{?(x}{}_1\text{,…,x}{}_{\mathrm{p}}\text{)}$ be a polynomial in the variables x1,…,xp with nonnegative real coefficients which sum to one, let A1,…,Ap be stochastic matrices, and let $\text{?}\text{?}\text{(A}{}_1\text{,…,A}{}_{\mathrm{p}}\text{)}$ be the stochastic matrix which is obtained from ? by substituting the Kronecker product of An11,…,Anppfor each term Xn11·?·Xnpp. In this paper, we present necessary and sufficient conditions for the Cesàro limit of the sequence of the powers of $\text{?}\text{?}\text{(A}{}_1\text{,…,A}{}_{\mathrm{p}}\text{)}$ to be equal to the Kronecker product of the Cesàro limits associated with each of A1,…,Ap. These conditions show that the equality of these two matrices depends only on the number of ergodic sets under$\text{?}\text{?}\text{(A}{}_1\text{,…,A}{}_{\mathrm{p}}\text{)}$ and?or the cyclic structure of the ergodic sets under A1,…,Ap, respectively. As a special case of these results, we obtain necessary and sufficient conditions for the interchangeability of the Kronecker product and the Cesàro limit operator.     

A canonical decomposition of the probability measure of sets of isotropic random points in Rn

The probability measure of X = (x₀,…, x_r), where x₀,…, x_r are independent isotropic random points in $\text{R}$ⁿ (1 ≤ r ≤ n ? 1) with absolutely continuous distributions is, for a certain class of distributions of X, expressed as a product measure involving as factors the joint probability measure of (ω, ?), the probability measure of p, and the probability measure of $\text{Y}{}^1\text{= (y}{}_0{}^1\text{,…, y}{}_{\mathrm{r}}{}^1\text{)}$. Here ω is the r-subspace parallel to the r-flat η determined by X, ? is a unit vector in ω^⊥ with ‘initial’ point at the origin [ω^⊥ is the (n ? r)-subspace orthocomplementary to ω], p is the norm of the vector z from the origin to the orthogonal projection of the origin on η, and $\text{y}{}_{\mathrm{i}}{}^1\text{=}\text{(x}{}_{\mathrm{i}}\text{? z)}\text{α(p}{}^2\text{)}$, where α is a scale factor determined by p. The probability measure for ω is the unique probability measure on the Grassmann manifold of r-subspaces in $\text{R}$ⁿ invariant under the group of rotations in $\text{R}$ⁿ, while the conditional probability measure of ? given ω is uniform on the boundary of the unit (n ? r)-ball in ω^⊥ with centre at the origin. The decomposition allows the evaluation of the moments, for a suitable class of distributions of X, of the r-volume of the simplicial convex hull of {x₀,…, x_r} for 1 ≤ r ≤ n.     

The asymptotic distribution of lattice points in hyperbolic space

Suppose x and y are two points in the upper half-plane H⁺, and suppose Γ is a discontinuous group of conformal automorphisms of H⁺ having compact fundamental domain S. Denote by N_T(x, y) the number of points of the form γy (γ?Γ) in the closed disc of hyperbolic radius T centered about x, and set $\text{Q}{}_{\mathrm{T}}\text{(x, y) = N}{}_{\mathrm{T}}\text{(x, y) ?}\text{V(T)}\text{A}\text{,}\text{where}\text{V(T)}$ is the hyperbolic area of the disc, and A is the hyperbolic area of S. The asymptotic behavior of the quantity ?_LxL(Q_T(x,y))² is estimated in terms of small eigenvalues of the Laplacian on functions automorphic under Γ.     

On Berry-Esséen rates,a law of the iterated logarithm and an invariance principle for the proportion of the sample below the sample mean

Let F_n(x) be the empirical distribution function based on n independent random variables X₁,…,X_n from a common distribution function F(x), and let $\text{X}\text{= Σ}{}_{\mathrm{i=1}}{}^{\mathrm{n}}\text{X}{}_{\mathrm{i}}\text{n}$ be the sample mean. We derive the rate of convergence of $\text{F}{}_{\mathrm{n}}\text{(}\text{X}\text{)}$ to normality (for the regular as well as nonregular cases), a law of iterated logarithm, and an invariance principle for $\text{F}{}_{\mathrm{n}}\text{(}\text{X}\text{)}$.     

A completeness theorem for a nonlinear multiparameter eigenvalue problem

In this paper we study the linked nonlinear multiparameter system

$\text{y}{}_{\mathrm{r}}{}^{\mathrm{n}}\text{(X}{}_{\mathrm{r}}\text{) + M}{}_{\mathrm{r}}\text{Y}{}_{\mathrm{r}}\text{+}\text{∑}\text{s=1}\text{k}\text{λ}{}_{\mathrm{s}}\text{(a}{}_{\mathrm{rs}}\text{(X}{}_{\mathrm{r}}\text{) + P}{}_{\mathrm{rs}}\text{) Y}{}_{\mathrm{r}}\text{(X}{}_{\mathrm{r}}\text{) = 0, r = l,…, k}$

, where x_r? [a_r, b_r], y_r is subject to Sturm-Liouville boundary conditions, and the continuous functions a_rs satisfy ¦ $\text{A ¦ (x) =}\text{det}\text{a}{}_{\mathrm{rs}}\text{(x}{}_{\mathrm{r}}\text{) > 0}$. Conditions on the polynomial operators M_r, P_rs are produced which guarantee a sequence of eigenfunctions for this problem yⁿ(x) = Π_r=1^ky_rⁿ(x_r), n ? 1, which form a basis in $\text{L}{}^2\text{([a, b], ¦ A ¦)}$. Here [a, b] = [a₁, b₁ × … × [a_k, b_k].     

Zeta matrices of elliptic curves

Let $\text{O =}\text{lim}{}_{\mathrm{n}}\text{Z/p}{}^{\mathrm{n}}\text{Z}$, let $\text{A}\text{= O[g}{}_2\text{, g}{}_3\text{]}{}_{\mathrm{\Delta }}$, where g₂ and g₃ are coefficients of the elliptic curve: Y² = 4X³ ? g₂X ? g₃ over a finite field and Δ = g₂³ ? 27g₃² and let $\text{B}\text{=}\text{A}\text{[X, Y]}\text{(Y}{}^2\text{? 4X}{}^3\text{+ g}{}_2\text{X + g}{}_3\text{)}$. Then the p-adic cohomology theory will be applied to compute explicitly the zeta matrices of the elliptic curves, induced by the pth power map on the free $\text{A}{}^2\text{?}{}_{\mathrm{Z}}\text{Q}$-module $\text{H}{}^1\text{(X,}\text{A}{}^2\text{?}{}_{\mathrm{Z}}\text{Q)}$. Main results are; Theorem 1.1: X²dY and YdX are basis elements for $\text{H}{}^1\text{(}\text{X}\text{, Γ}{}_{\mathrm{<mtext>A</mtext>}}{}^1\text{(}\text{X}\text{)}{}^2\text{?}{}_{\mathrm{Z}}\text{Q)}$; Theorem 1.2: YdX, X²dY, Y^?1dX, Y^?2dX and XY^?2dX are basis elements for $\text{H}{}^1\text{(}\text{X}\text{? (Y = 0), Γ}{}_{\mathrm{<mtext>A</mtext>}}{}^1\text{(}\text{X}\text{)}{}^2\text{?}{}_{\mathrm{Z}}\text{Q)}$, where $\text{X}$ is a lifting of X, and all the necessary recursive formulas for this explicit computation are given.     

Decomposability of symmetric multilinear functions

Let V denote a finite dimensional vector space over a field K of characteristic 0, let T_n(V) denote the vector space whose elements are the K-valued n-linear functions on V, and let S_n(V) denote the subspace of T_n(V) whose members are the fully symmetric members of T_n(V). If $\text{L}$_n denotes the symmetric group on {1,2,…,n} then we define the projection $\text{P}{}_{\mathrm{<mtext>L</mtext>}}\text{: T}{}_{\mathrm{n}}\text{(V) → S}{}_{\mathrm{n}}\text{(V)}$ by the formula , where P_σ : T_n(V) → T_n(V) is defined so that P_σ(A)(y₁,y₂,…,y_n = A(y_σ(1),y_σ(2),…,y_σ(n)) for each A?T_n(V) and y_i?V, 1 ? i ? n. If $\text{x}{}_{\mathrm{i}}\text{? V}{}^1\text{, 1 ? i ? n}$, then x₁?x₂? … ?x_n denotes the member of T_n(V) such that $\text{(x}{}_1\text{?x}{}_2\text{· ? ? ?x}{}_{\mathrm{n}}\text{)(y}{}_1\text{,y}{}_2\text{,…,y}{}_{\mathrm{n}}\text{) = П}{}^{\mathrm{n}}{}_{\mathrm{i=1}}\text{x}{}_{\mathrm{i}}\text{(y}{}_{\mathrm{i}}\text{)}$ for each y₁ ,₂,…,y_n in V, and x₁·x₂… x_n denotes $\text{P}{}_{\mathrm{<mtext>L</mtext>}}\text{(x}{}_1\text{?x}{}_2\text{? … ?x}{}_{\mathrm{n}}\text{)}$. If B? S_n(V) and there exists $\text{x}{}_{\mathrm{i}}\text{? V}{}^1\text{, 1 ? i ? n}$, such that B = x₁·x₂…x_n, then B is said to be decomposable. We present two sets of necessary and sufficient conditions for a member B of S_n(V) to be decomposable. One of these sets is valid for an arbitrary field of characteristic zero, while the other requires that K = R or C.     

Some asymptotic formulas on generalized divisor functions,II

Let A be an infinite sequence of positive integers a₁ < a₂ <… and put $\text{f}{}_{\mathrm{A}}\text{(x) =}\text{Σ}{}_{\mathrm{a\in A,\; a\le x}}\text{(}\text{1}\text{a}\text{)}$, $\text{D}{}_{\mathrm{A}}\text{(x) =}\text{max}{}_{\mathrm{1\le n\le x}}\text{Σ}{}_{\mathrm{<mtext>a\in A,</mtext><mtext>a</mtext><mtext>n</mtext>}}\text{1}$. In Part I, it was proved that $\text{lim}{}_{\mathrm{x\to +\infty }}\text{sup}\text{D}{}_{\mathrm{A}}\text{(x)}\text{f}{}_{\mathrm{A}}\text{(x)}\text{= +∞}$. In this paper, this theorem is sharpened by estimating D_A(x) in terms of f_A(x). It is shown that lim_x→+∞sup D_A(x) exp(?c₁(logf_A(x))²) = +∞ and that this assertion is not true if c₁ is replaced by a large constant c₂.     

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

For (x,y,t)∈$\text{R}$ⁿ × $\text{R}$ⁿ × $\text{R}$, denote $\text{X}{}_{\mathrm{j}}\text{=}\text{?}\text{?x}{}_{\mathrm{j}}\text{+ 2y}{}_{\mathrm{j}}\text{?}\text{?t}\text{, y}{}_{\mathrm{j}}\text{=}\text{?}\text{?y}{}_{\mathrm{j}}\text{? 2x}{}_{\mathrm{j}}\text{?}\text{?t}$ and $\text{L}{}_{\mathrm{\alpha }}\text{=?}\text{1}\text{4}\text{∑}\text{j=1}\text{n}\text{X}{}_{\mathrm{j}}{}^2\text{+ Y}{}_{\mathrm{j}}{}^2\text{+}\text{?}\text{?t}$. When α = n ? 2q, $\text{L}$_a represents the action of the Kohn Laplacian □_b on q-forms on the Heisenberg group. For ?n < α < n, we construct a parametrix for the Dirichlet problem in smooth domains D near non-characteristic points of ?D. A point w of ?D is non-characteristic if one of X₁,…, X_n, Y₁,…, Y_n is transverse to ?D at w. This yields sharp local estimates in the Dirichlet problem in the appropriate non-isotropic Lipschitz classes. The main new tool is a “convolution calculus” of pseudo-differential operators that can be applied to the relevant layer potentials, for which the usual asymptotic composition formula is false. Characteristic points are treated in Part II.     

Some special vapnik-chervonenkis classes

For a class $\text{C}$ of subsets of a set X, let V($\text{C}$) be the smallest n such that no n-element set F?X has all its subsets of the form A ∩ F, A ∈ $\text{C}$. The condition V($\text{C}$) <+∞ has probabilistic implications. If any two-element subset A of X satisfies both A ∩ C = Ø and A ? D for some C, D∈$\text{C}$, then $\text{V(}\text{C}\text{)=2}$ if and only if $\text{C}$ is linearly ordered by inclusion. If $\text{C}$ is of the form $\text{C}\text{={∩}{}^{\mathrm{n}}{}_{\mathrm{i=1}}\text{C}{}_{\mathrm{i}}\text{:C}{}_{\mathrm{i}}\text{∈}\text{C}{}_{\mathrm{i}}$, i=1,2,…,n}, where each $\text{C}{}_{\mathrm{i}}$ is linearly ordered by inclusion, then $\text{V(}\text{C}\text{)?n+1}$. If H is an (n-1)-dimensional affine hyperplane in an n-dimensional vector space of real functions on X, and $\text{C}$ is the collection of all sets {x: f(x)>0} for f in H, then $\text{V(}\text{C}\text{)=n}$.     

Arnold's canonical matrices and the asymptotic simplification of ordinary differential equations

Let A(x,ε) be an n×n matrix function holomorphic for |x|?x₀, 0<ε?ε₀, and possessing, uniformly in x, an asymptotic expansion $\text{A(x,ε)?Σ}{}^{\mathrm{\infty }}{}_{\mathrm{r=0}}\text{A}{}_{\mathrm{r}}\text{(x) ε}{}^{\mathrm{r}}$, as ε→0⁺. An invertible, holomorphic matrix function P(x,ε) with an asymptotic expansion $\text{P(x,ε)?Σ}{}^{\mathrm{\infty }}{}_{\mathrm{r=0}}\text{P}{}_{\mathrm{r}}\text{(x)ε}{}^{\mathrm{r}}$, as ε→0⁺, is constructed, such that the transformation y = P(x,ε)z takes the differential equation $\text{ε}{}^{\mathrm{h}}\text{dy}\text{dx}\text{= A(x,ε)y,h}$ a positive integer, into $\text{ε}{}^{\mathrm{h}}\text{dz}\text{dx}\text{= B(x,ε)z}$, where B(x,ε) is asymptotically equal, to all orders, to a matrix in a canonical form for holomorphic matrices due to V.I. Arnold.