On Piecewise-Constant Approximation of Continuous Functions of n Variables in Integral Metrics

We consider the approximation by piecewise-constant functions for classes of functions of many variables defined by moduli of continuity of the form (₁, ..., _n ) = ₁(₁) + ... + _n ( _n ), where _i ( _i ) are ordinary moduli of continuity that depend on one variable. In the case where _i ( _i ) are convex upward, we obtain exact error estimates in the following cases: (i) in the integral metric L ₂ for (₁, ..., _n ) = ₁(₁) + ... + _n ( _n ); (ii) in the integral metric L _p (p 1) for (₁, ..., _n ) = c ₁₁ + ... + c _{n n} ; (iii) in the integral metric L _{(2, ..., 2, 2r)} (r = 2, 3, ...) for (₁, ..., _n ) = ₁(₁) + ... + _{n – 1}( _{n – 1}) + c _{n n} .     

Lower functions for asymmetric Lévy processes

Summary Let {X _t } be a ¹ process with stationary independent increments and its Lévy measurev be given byv{yy>x}=x ^–L ₁ (x), v{yy<–x}=x ^–L ₂ (x) whereL ₁,L ₂ are slowly varying at 0 and and 0<1. We construct two types of a nondecreasing functionh(t) depending on 0<<1 or =1 such that lim inf a.s. ast 0 andt for some positive finite constantC.This research is partialy supported by a grant from Korea University     

General technique for solving nonlinear,two-point boundary-value problems via the method of particular solutions

In this paper, a general technique for solving nonlinear, two-point boundary-value problems is presented; it is assumed that the differential system has ordern and is subject top initial conditions andq final conditions, wherep+q=n. First, the differential equations and the boundary conditions are linearized about a nominal functionx(t) satisfying thep initial conditions. Next, the linearized system is imbedded into a more general system by means of a scaling factor , 01, applied to each forcing term. Then, themethod of particular solutions is employed in order to obtain the perturbation x(t)=A(t) leading from the nominal functionx(t) to the varied function (t); this method differs from the adjoint method and the complementary function method in that it employs only one differential system, namely, the nonhomogeneous, linearized system.The scaling factor (or stepsize) is determined by a one-dimensional search starting from =1 so as to ensure the decrease of the performance indexP (the cumulative error in the differential equations and the boundary conditions). It is shown that the performance index has a descent property; therefore, if is sufficiently small, it is guaranteed that <P. Convergence to the desired solution is achieved when the inequalityP is met, where is a small, preselected number.In the present technique, the entire functionx(t) is updated according to (t)=x(t)+A(t). This updating procedure is called Scheme (a). For comparison purposes, an alternate procedure, called Scheme (b), is considered: the initial pointx(0) is updated according to (0)=x(0)+A(0), and the new nominal function (t) is obtained by forward integration of the nonlinear differential system. In this connection, five numerical examples are presented; they illustrate (i) the simplicity as well as the rapidity of convergence of the algorithm, (ii) the importance of stepsize control, and (iii) the desirability of updating the functionx(t) according to Scheme (a) rather than Scheme (b).This research, supported by the National Science Foundation, Grant No. GP-18522, is based on Ref. 1. The authors are indebted to Mr. A. V. Levy for computational assistance.     

Solution of the nonlinear functional equations representing the roll gap relationships in a cold mill

In the development of a roll force model for cold rolling, techniques were developed for solving the system equations which are of general interest. This paper gives a brief introduction to the physical model but concentrates on the solution of the model equations and the simulation. An unusual feature of the model was that the calculated profiles had to satisfy a number of boundary conditions at different points throughout the roll arc. A new method was developed for calculating these profiles and for determining the gradient functions which satisfied the boundary constraints.Nomenclature p() pressure at roll angle - h() gauge - a() roll radius - y() yield stress - g _i () gradient function on iterationi - e() gauge error - (, ) transition function - H(–) Heaviside unit step function at = - (–) unit impulse function at = - H(, ₁, ₂) defined asH( – ₁) –H( – ₂) - angular position from the roll center line - _T angular limits of roll arc represented - _n angular position of the neutral angle - _i angular position ofith strip elastic-plastic boundary - _pi pressure change at the boundaryi - _i , _i , _i constants defined in Appendix A - k ₁,k ₂ elastic region constants - k total number of strip boundaries (elastic-plastic and entry and exit points) - R undeformed work roll radius - R _s roll separation—distance between roll centers - h ₀₁ unstrained gauge in an elastic region - h _in gauge of the strip at the entry to the roll gap - J gauge error cost function - <x, y> inner product ofx andy - x norm ofx - L ₂[0, _T ] the space of Lebesgue square-integrable functions defined on the interval [0, _T ] - JUVY denotes (Dx)() =dx()/d The author would like to acknowledge the help given by Dr. G. F. Bryant, Director, and Mr. M. A. Fuller, Senior Research Engineer, the Industrial Automation Group, Imperial College of Science and Technology, London. He is also grateful to M. J. G. Henderson of the University of Birmingham for his advice and encouragement during the project. He would like to thank the Directors of GEC Electrical Projects Limited for allowing him to undertake the work and also Mr. J. McTaggart and Mr. C. McKenzie (GEC), Professor H. A. Prime of the University of Birmingham, and Dr. G. F. Bryant for arranging the project.     

Superstable implicit Runge-Kutta methods for second order initial value problems

Using the well known properties of thes-stage implicit Runge-Kutta methods for first order differential equations, single step methods of arbitrary order can be obtained for the direct integration of the general second order initial value problemsy=f(x, y, y),y(x _o)=y _o,y(x _o)=y _o. These methods when applied to the test equationy+2y+ ² y=0, ,0, +>0, are superstable with the exception of a finite number of isolated values ofh. These methods can be successfully used for solving singular perturbation problems for which f/y and/or f/y are negative and large. Numerical results demonstrate the efficiency of these methods.     

On the Eigenvalues and Eigenfunctions of the Sturm--Liouville Operator with a Singular Potential

In this paper we consider the Sturm--Liouville operators generated by the differential expression -y+q(x)y and by Dirichlet boundary conditions on the closed interval [0,]. Here q(x) is a distribution of first order,, i.e., q(x)dx L ₂[0,]. Asymptotic formulas for the eigenvalues and eigenfunctions of such operators which depend on the smoothness degree of q(x) are obtained.     

Detailed Error Analysis for a Fractional Adams Method

We investigate a method for the numerical solution of the nonlinear fractional differential equation D _* y(t)=f(t,y(t)), equipped with initial conditions y ^(k)(0)=y ₀ ^(k), k=0,1,...,–1. Here may be an arbitrary positive real number, and the differential operator is the Caputo derivative. The numerical method can be seen as a generalization of the classical one-step Adams–Bashforth–Moulton scheme for first-order equations. We give a detailed error analysis for this algorithm. This includes, in particular, error bounds under various types of assumptions on the equation. Asymptotic expansions for the error are also mentioned briefly. The latter may be used in connection with Richardson's extrapolation principle to obtain modified versions of the algorithm that exhibit faster convergence behaviour.     

Blow-up rates for parabolic systems

Let ⁿ be a bounded domain andB _R be a ball in ⁿ of radiusR. We consider two parabolic systems: ut=u +f(), i= +g(u) in × (0,T) withu=v=0 on × (0,T) andu _t =u, v _t =v inB _r × (0,T) withe/v=f (v), e/v=g(u) onB _R × (0,T). Whenf(v) andg(u) are power law or exponential functions, we establish estimates on the blow-up rates for nonnegative solutions of the systems.     

A phase transition for the coupled branching process

Summary We consider a particular Markov process _t ^u on ^S ,S= ⁿ . The random variable _t ^u (x) is interpreted as the number of particles atx at timet. The initial distribution of this process is a translation invariant measure withf(x)d<. The evolution is as follows: At rateb(x) a particle is born atx but moves instantaneously toy chosen with probabilityq(x, y). All particles at a site die at ratepd withp[0, 1],d, ⁺ and individual particles die independently from each other at rate (1–p)d. Every particle moves independently of everything else according to a continuous time random walk.We are mainly interested in the caseb=d andn3. The process exhibits a phase transition with respect to the parameterp: Forp<p ^* all weak limit points of ( _t ^µ ) ast still have particle density (x)d. Forp>p ^*, _t ^µ ) converges ast to the measure concentrated on the configuration identically 0. We calculatep ^* as well asp ⁽ⁿ⁾ , the points with the property that the extremal invariant measures have forp>p ⁽ⁿ⁾ infiniten-th moment of (x) and forp<p ⁽ⁿ⁾ finiten-th moment. We show the case 1>p ^*>p⁽²⁾>p⁽³⁾...p ⁽ⁿ⁾ >0, p⁽ⁿ⁾0 occurs for suitable values of the other parameters. Forp<p ⁽²⁾ we prove the system has a one parameter set of extremal invariant measures and we determine their domain of attraction. Part I contains statements of all results but only the proofs of the results about the process for values ofp withp<p ⁽²⁾ and the behaviour of then-th moments andp ⁽ⁿ⁾ .     

A Nonparametric Estimation Problem From Indirect Observations

This paper deals with a nonparametric estimation problem of an integral-type functional from indirect observations where the observation Y (t) is a sum of a known function of an unobservable process X (t) and a Gaussian white noise, and X (t) is a sum of an unknown function a(t) and a Gaussian process. The minimax lower bound on the quality of nonparametric estimation is derived and an asymptotically efficient estimator is proposed. The paper concludes with some examples including one about reduction to parameter estimation.     

Path properties of kernel generated two-time parameter Gaussian processes

Summary We study path properties of two-parameter Gaussian processes {X(t,v),t R} of the formX(t,v)= _{0 –} (t,v,x,y)dW(x,y), where the kernel function (t, v, x, y) is assumed to be square integrable in (x, y) onR×R ⁺, andW(x, y) is a standard two-parameter Wiener process.Work partially supported by an NSERC Canada Operating Grant at Carleton UniversityWork supported by an NSERC Canada International Scientific Exchange Award at Carleton University and by National Natural Science Foundation of China     

Intuitionistic weak arithmetic

We construct -framed Kripke models of i₁ and i₁ non of whose worlds satisfies xy(x=2yx=2y+1) and x,yzExp(x, y, z) respectively. This will enable us to show that i₁ does not prove ¬¬xy(x=2yx=2y+1) and i₁ does not prove ¬¬x, yzExp(x, y, z). Therefore, i₁¬¬lop and i₁¬¬i₁. We also prove that HAl₁ and present some remarks about i₂. Mathematics Subject Classification (2000):03F30, 03F55, 03H15.     

Phase-lag analysis of explicit Nyström methods fory″=f(x,y)

We examine phase-lag (frequency distortion) of the two-parameter familyM ₄(₁, ₃) of fourth order explicit Nyström methods of [1] by applying these to the test equation:y+ ² y=0, >0. While the methodM ₄(1/6, 5/6) possessing the largest interval of periodicity of size 3.46 has a phase-lag of (1/4320)H (H ⁴=h, h is the step-size), we show that there exist two fourth order methods ofM ₄(₁, ₃) for which the phase-lag is minimal and of size (1/40320)H ⁶; interestingly, both methods also possess a sizable interval of periodicity of length 2.75 each.     

Approximation and regular perturbation of optimal control problems via Hamilton-Jacobi theory

We present two convergence theorems for Hamilton-Jacobi equations and we apply them to the convergence of approximations and perturbations of optimal control problems and of two-players zero-sum differential games. One of our results is, for instance, the following. LetT andT _h be the minimal time functions to reach the origin of two control systemsy = f(y, a) andy = f _h (y, a), both locally controllable in the origin, and letK be any compact set of points controllable to the origin. If f _h –f Ch, then |T(x) – T _h (x)| C _K h , for all x K, where is the exponent of Hölder continuity ofT(x).     

A complete asymptotic expansion of the jacobi functions with error bounds

In this paper we give a complete asymptotic expansion of the Jacobi functions ^{(, )} (t) as + . The method we employed to get the complete expansion follows that of Olver in treating similar problems. By using a Gronwall-Bellman type inequality for an improper integral in which the integrand is an unbounded function and contains a parameter, we get an error bound of the asymptotic approximation which is different from that of Olver's.     

Hereditary optimal control problems: Numerical method based upon a padé approximation

In this paper, we consider a particular approximation scheme which can be used to solve hereditary optimal control problems. These problems are characterized by variables with a time-delayed argumentx(t – ). In our approximation scheme, we first replace the variable with an augmented statey(t) x(t - ). The two-sided Laplace transform ofy(t) is a product of the Laplace transform ofx(t) and an exponential factor. This factor is approximated by a first-order Padé approximation, and a differential relation fory(t) can be found. The transformed problem, without any time-delayed argument, can then be solved using a gradient algorithm in the usual way. Four problems are solved to illustrate the validity and usefulness of this technique.This research was supported in part by the National Aeronautics and Space Administration under NASA Grant NCC-2-106.     

A general minimax theorem

This paper is concerned with minimax theorems for two-person zero-sum games (X, Y, f) with payofff and as main result the minimax equality inf supf (x, y)=sup inff (x, y) is obtained under a new condition onf. This condition is based on the concept of averaging functions, i.e. real-valued functions defined on some subset of the plane with min {x, y}< (x, y)x, y} forx y and (x, x)=x. After establishing some simple facts on averaging functions, we prove a minimax theorem for payoffsf with the following property: Forf there exist averaging functions and such that for any x₁, x₂ X, > 0 there exists x₀ X withf (x₀, y) > f (x₁,y),f (x₂,y))– for ally Y, and for any y₁, y₂ Y, > 0 there exists y₀ Y withf (x, y₀) (f (x, y₁),f (x, y₂))+. This result contains as a special case the Fan-König result for concave-convex-like payoffs in a general version, when we take linear averaging with (x, y)=x+(1–)y, (x, y)=x+(1–)y, 0 <, < 1.Then a class of hide-and-seek games is introduced, and we derive conditions for applying the minimax result of this paper.

---

Zusammenfassung In dieser Arbeit werden Minimaxsätze für Zwei-Personen-Nullsummenspiele (X, Y,f) mit Auszahlungsfunktionf behandelt, und als Hauptresultat wird die Gültigkeit der Minimaxgleichung inf supf (x, y)=sup inff (x, y) unter einer neuen Bedingung an f nachgewiesen. Diese Bedingung basiert auf dem Konzept mittelnder Funktionen, d.h. reellwertiger Funktionen, welche auf einer Teilmenge der Ebene definiert sind und dort der Eigenschaft min {x, y} < < (x, y)x, y} fürx y, (x, x)=x, genügen. Nach der Herleitung einiger einfacher Aussagen über mittelnde Funktionen beweisen wir einen Minimaxsatz für Auszahlungsfunktionenf mit folgender Eigenschaft: Zuf existieren mittelnde Funktionen und, so daß zu beliebigen x₁, x₂ X, > 0 mindestens ein x₀ X existiert mitf (x₀,y) (f (x ₁,y),f (x₂,y)) – für alley Y und zu beliebigen y₁, y₂ Y, > 0 mindestens ein y₀ Y existiert mitf (x, y₀) (f (x, y₁),f (x, y ₂))+ für allex X. Dieses Resultat enthält als Spezialfall den Fan-König'schen Minimaxsatz für konkav-konvev-ähnliche Auszahlungsfunktionen in einer allgemeinen Version, wenn wir lineare Mittelung mit (x, y)=x+(1–)y, (x, y)= x+(1–)y, 0 <, < 1, betrachten.Es wird eine Klasse von Suchspielen eingeführt, welche mit dem vorstehenden Resultat behandelt werden können.

     

On the distribution of square-full and cube-full integers

LetN _r (x) be the number ofr-full integers x and let _r (x) be the error term in the asymptotic formula forN _r (x). Under Riemann's hypothesis, we prove the estimates ₂(x)x^1/7+, ₃(x)x^97/804+(>0), which improve those of Cao and Nowak. We also investigate the distribution ofr-full andl-free numbers in short intervals (r=2,3). Our results sharpen Krätzel's estimates.     

Stochastic processes with Gaussian interaction of components

Summary LetU(x), x ^d-|0}, be a nonnegative even function such that _x 0U(x)1. In this paper, we consider an infinite system of stochastic process _t (x); x ^d with the following mechanism: at each sitex, after mean 1 exponential waiting time, _t(x) is replaced by a Gaussian random variable with mean _{yx t} (y) U(y-x) and variance 1. It is understood here that all the interactions are independent of one another. The behavior of this system will be investigated and some ergodic theorems will be derived. The results strongly depend whether _{x 0} U(x)<1 or =1.     

Transformation of Wiener measure under anticipative flows

Summary LetT()=+F() be a transformation from the Wiener space to itself with the range ofF() assumed to be in the Cameron-Martin space. The absolute continuity and the density function associated withT is considered;T is assumed to be embedded in or defined through a parameterizationT _t =+F _t () andF _t is assumed to be differentiable int. The paper deals first with the case where the range of thet-derivative ofF _t () is also in the Cameron-Martin space and new representations for the Radon-Nikodym derivative and the Carleman-Fredholm determinant are derived. The case where thet-derivative ofF _t is not in the Cameron-Martin space is considered next and results on the absolute continuity and the density function, under conditions which are considerably weaker than previously known conditions, are presented.The work of the second author was supported by the fund for promotion of research at the Technion