Monotonicity and boundedness in implicit Runge-Kutta methods

Summary This paper concerns the analysis of implicit Runge-Kutta methods for approximating the solutions to stiff initial value problems. The analysis includes the case of (nonlinear) systems of differential equations that are essentially more general than the classical test equationU=U (with a complex constant). The properties of monotonicity and boundedness of a method refer to specific moderate rates of growth of the approximations during the numerical calculations. This paper provides necessary conditions for these properties by using the important concept of algebraic stability (introduced by Burrage, Butcher and by Crouzeix). These properties will also be related to the concept of contractivity (B-stability) and to a weakened version of contractivity.     

Stability analysis of numerical methods for delay differential equations

Summary This paper deals with the stability analysis of step-by-step methods for the numerical solution of delay differential equations. We focus on the behaviour of such methods when they are applied to the linear testproblemU(t)=U(t)+U(t–) with >0 and , complex. A general theorem is presented which can be used to obtain complete characterizations of the stability regions of these methods.     

Numerical implementation of coherence for the example of the “Swiss Cross”

Summary We consider the two-dimensional Helmholtz equation u+u=0 inD with the boundary conditionsu=0 on D. D is the Swiss Cross — a region consisting of five unit squares. A method based on the concept of Coherence is utilized to determine an approximation for the first eigenvalue= ₁ more accurate than calculated by classical difference methods. The numerical result is used to illustrate isoperimetric upper and lower bounds for ₁, and to test some conjectures on its relations with torsional rigidity.Dedicated to the memory of Professor Lathar Collatz     

Subordination des résolvantes

Résumé Soit (V ₎₀ une résolvante définie sur un espace mesurable telle que le noyau initial est borné; on trouve une condition nécéssaire et suffisante pour qu'un noyau borné U possède une résolvante (U ₎₀ telle que U V pour tout 0. On donne plusieurs applications de ce résultat.     

The dependence of critical parameter bounds on the monotonicity of a Newton sequence

Summary A new method is proposed for the inclusion of the critical parameter ^* of some convex operator equationu=Tu (appearing e.g. in thermal explosion theory). It is based on the fact that for a fixed Newton's method starting with a suitable subsolution is not monotonically if and only if >^*. Several numerical examples arising from nonlinear boundary value problems illustrate the efficiency of the method.     

Ultrafilter translations

We develop a method for extending results about ultrafilters into a more general setting. In this paper we shall be mainly concerned with applications to cardinality logics. For example, assumingV=L, Gödel's Axiom of Constructibility, we prove that if > then the logic with the quantifier there exist many is (,)-compact if and only if either is weakly compact or is singular of cofinality<. As a corollary, for every infinite cardinals and , there exists a (,)-compact non-(,)-compact logic if and only if either < orcf<cf or < is weakly compact.Counterexamples are given showing that the above statements may fail, ifV=L is not assumed.However, without special assumptions, analogous results are obtained for the stronger notion of [,]-compactness.     

A sandwich theorem and Hyers—Ulam stability of affine functions

It is shown that two real functionsf andg, defined on a real intervalI, satisfy the inequalitiesf(x + (1 – )y) g(x) + (1 – )g(y) andg(x + (1 – )y) f(x) + (1 – )f(y) for allx, y I and [0, 1], iff there exists an affine functionh: I such thatf h g. As a consequence we obtain a stability result of Hyers—Ulam type for affine functions.     

Reverse Functional Inequalities and Their Applications to Nonlinear Elliptic Boundary Value Problems

We establish some reverse inequalities. We give applications to nonlinear elliptic boundary value problems containing a parameter which have two branches of solutions u (0) and U (>0) of which the first is continuous at the origin and the second increases indefinitely as 0.     

A new interpolation procedure for adapting Runge-Kutta methods to delay differential equations

This paper deals with adapting Runge-Kutta methods to differential equations with a lagging argument. A new interpolation procedure is introduced which leads to numerical processes that satisfy an important asymptotic stability condition related to the class of testproblemsU(t)=U(t)+U(t–) with , C, Re()<–||, and >0. Ifc _i denotes theith abscissa of a given Runge-Kutta method, then in thenth stept _n–1t _n :=t _n–1+h of the numerical process our interpolation procedure computes an approximation toU(t _n–1+c _i h–) from approximations that have already been generated by the process at pointst _j–1+c _i h(j=1,2,3,...). For two of these new processes and a standard process we shall consider the convergence behaviour in an actual application to a given, stiff problem.     

Shift invariant spaces inL 2

If is a complex, separable Hilbert space, letL ² () denote theL ²-space of functions defined on the unit circle and having values in . The bilateral shift onL ²() is the operator (U f)()=f(). A Hilbert spaceH iscontractively contained in the Hilbert spaceK ifHK and the inclusion mapH–K is a contraction. We describe the structure of those Hilbert spaces, contractively contained inL ²(), that are carried into themselves contractively byU . We also do this for the subcase of those spaces which are carried into themselves unitarily byU .     

Resolution of the Symmetric Nonnegative Inverse Eigenvalue Problem for Matrices Subordinate to a Bipartite Graph

There is a symmetric nonnegative matrix A, subordinate to a given bipartite graph G on n vertices, with eigenvalues _{12 n} if and only if, ₁ + _n 0, ₂ + _n-10,..., _m + _{n - m + 1}0, _{m + 1}0,..., _{n - m} 0, in which m is the matching numberof G. Other observations are also made about the symmetric nonnegative inverse eigenvalue problem with respect to a graph     

Diagonals of positive semigroups

LetE be a Dedekind complete complex Banach lattice and letD denote the diagonal projection from the spaceL _r (E) onto the centerZ(E) ofE. Let {T(t)} _t0 be a positive strongly continuous semigroup of linear operators with generatorA. The first main result is that if the spectral bounds(A) equals to zero, then the functionD(T(t)) is a center valuedp-function. The second main result is that if for >0 the diagonalD(R(, A)) of the resolvent operatorR(, A) is strictly positive, then (D(R(, A))) ^–1 is a center valued Bernstein function. As an application of these results it follows that the order limit lim₀D(R(,A)) exists inZ(E) and equals the order limit lim _m D((R(, A)) ^m ) for any >0.     

Remarks on the perturbation of analytic matrix functions III

As in [N], [LN] the Newton diagram is used in order to get information about the first terms of the Puiseux expansions of the eigenvalues () of the perturbed matrix pencilT(, )=A()+B(, ) in the neighbourhood of an unperturbed eigenvalue () ofA(). In fact sufficient conditions are given which assure that the orders of these first terms correspond to the partial multiplicities of the eigenvalue ₀ ofA().     

Spectral parameter asymptotics of the Weyl solutions of Sturm-Liouville equations

In the paper one investigates the dependence of Weyl's solution ,)=c(,)+n()s(,) of the Sturm-Liouville equation y+q()y=²y on the spectral parameter . Under the condition that the potential q is bounded from below and q()exp(c₀+c[in1 ¦¦), it is proved for {ie217-01} for any positive values and A. If q()>1 and {ie217-02} for all >0, then in the semiplane >0 the Weyl solution (, ) is obtained from the Weyl solution (,x) is obtained from the Weyl solution e^ix with zero potential, with the aid of a generalization of B. Ya Levin's transformation operators.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova Akademii Nauk SSSR, Vol. 170, pp. 184–206, 1989.I express my sincere gratitude to L. A. Pastur and I. V. Ostrovskii for valuable advice and discussions.     

Preassigning the boundary of diametrically-complete sets

In a Minkowski space with unit ballS, a setD() is diametrically complete ifD()={. In this paper a method is given for generating such setsD() which contain an arbitrary setX. This technique is then employed to find necessary and sufficient conditions for a setA to lie on the boundary of some suchD() and bounds are given on the acceptable values of .     

Principal part of resolvent and factorization of an increasing nonanalytic operator-function

Let a selfadjoint operator-valued functionL() be given on the interval [a,b] such thatL(a)0,L(b)0,L()0 (ab), andL() has a certain smoothness (for instance, it satisfies Hölder's condition). It turns out that the spectral theory of the operator-valued functionL() can be reduced to the spectral theory of one operatorZ, the spectrum of which lies on (a, b) and which is similar to a selfadjoint operator. In particular, the factorization takes place:L()=M()(I–Z), where the operator-valued functionM() is invertible on [a, b]. Earlier similar results were known only for analytic operator-valued functions. The authors had to use new methods for the proof of the described theorem. The key moment is the decomposition ofL ^–1() into the sume of its principal and regular parts.     

Topological influence of the presence of fully-λ-bases in locally convex spaces

This article reveals the topological impact of fully--bases in locally convex spaces where carries either the traditional normal topology or the fairly generalized-topology of Ruckle. It has been established that the generalized nuclearity of plays a significant role in influencing the topology of the space. Further, the equivalence of normal topology and the topology arising out of the fully--base ( being equipped with normal topology or-topology) has been investigated.We acknowledge with thanks the suggestions of the referee.     

The spectral bound of Schrödinger operators

Let V: R ^N [0, ] be a measurable function, and >0 be a parameter. We consider the behaviour of the spectral bound of the operator 1/2–V as a function of . In particular, we give a formula for the limiting value as , in terms of the integrals of V over subsets of R ^N on which the Laplacian with Dirichlet boundary conditions has prescribed values. We also consider the question whether this limiting value is attained for finite .     

Numerical solution of eigentuple-eigenvector problems in Hilbert spaces by a gradient method

Summary A gradient technique previously developed for computing the eigenvalues and eigenvectors of the general eigenproblemAx=Bx is generalized to the eigentuple-eigenvector problem . Among the applications of the latter are (1) the determination of complex (,x) forAx=Bx using only real arithmetic, (2) a 2-parameter Sturm-Liouville equation and (3) -matrices. The use of complex arithmetic in the gradient method is also discussed. Computational results are presented.This research was partially supported by NSF Grants MPS74-13332 and MCS76-09172     

On the stability properties of Brown's multistep multiderivative methods

Summary Brown introducedk-step methods usingl derivatives. We investigate for whichk andl the methods are stable or unstable. It is seen that to anyl the method becomes unstable fork large enough. All methods withk2(l+1) are stable. Fork=1,2,..., 18 there exists a _k such that the methods are stable for anyl _k and unstable for anyl < _k . The _k are given.