variable step size predictor corrector schemes for second kind volterra integral equations in this paper a family of implicit multistep methods for the solution of volterra integral equations is derived these methods together with an explicit euler predictor permit the use of a variable step size when solving integral equations means of controlling the error and stability by varying the step size and the order of the method are described extensive numerical results are presented

time finite element discretization of hamilton s law of varying action hamilton s law of varying action is used as a variational source for the derivation of finite element discretization procedure in the time domain three different versions of the proposed algorithms are presented and verified for accuracy and stability the first one is the high precision finite time element analogous to the standard finite elements with cutoff frequency the second version is the step by step one time element from which the unconditionally stable with slightly altered accuracy third algorithm is derived the new operator connected with the proposed algorithms bears attractive properties of much greater accuracy than other existing stable methods and easy computer implementation thus the work herein shows that the reservations expressed against the use of finite elements in time domain seem unjustified

multiple periodic solutions of hamiltonian systems with prescribed energy abstract consider the periodic solutions of autonomous hamiltonian systems x j h x on the given compact energy hypersurface h 1 1 if is convex or star shaped there have been many remarkable contributions for existence and multiplicity of periodic solutions it is a hard problem to discuss the multiplicity on general hypersurfaces of contact type in this paper we prove a multiplicity result for periodic solutions on a special class of hypersurfaces of contact type more general than star shaped ones

a multi purpose system for the numerical integration of ode s in this paper we describe the development of an experimental multi purpose integration system strongly modular the system allows the use of wide classes of multistep methods for solving ode s singulsr perturbation problems as well as implicit differential algebraic equations we focus our attention on a new family of collocation methods for nonstiff problems while these methods are strongly related to the stiffly stable bdf methods they are particularly well suited for the implementation of a type insensitive integration algorithm compound discretizations in which certain variables are integrated with the first method and others with the second can also be used a stability analysis is given and some implementation issues are discussed

closed characteristics on partially symmetric compact convex hypersurfaces in r2n abstract in this paper we study the multiplicity of closed characteristics on partially symmetric convex compact hypersurfaces in r 2n the main ingredient of the proof is a new p index function and its iteration theory for symplectic matrix paths and some symplectic orthogonal matrix p this theory yields estimates on this multiplicity via index estimates for a new type of iterations of any symplectic matrix path 0 sp 2n defined by t p t p for any t 0

the dynamics on three dimensional strictly convex energy surfaces we show that a hamiltonian flow on a three dimensional strictly convex energy surface s c r4 possesses a global surface of section of disc type it follows in particular that the number of its periodic orbits is either 2 or oc by a recent result of j franks on area preserving homeomorphisms of an open annulus in the plane the construction of this surface of section is based on partial differential equations of cauchy riemann type for maps from punctured riemann surfaces into r x s3 equipped with special almost complex structures

on one leg multistep methods the properties and possibilities of one leg methods are presented in a manner that admits the generalization to smoothly varying step size a new definition is given of the local truncation error and of the order of consistency some data are given for the most accurate one leg methods i e for p k and p k 1 where p is the order of consistency and k is the step number of the method

type insensitive ode codes based on implicit a stable formulas a special concept of stiffness is appropriate for implict a stable formulas it is possible to recognize this kind of stiffness economically and reliably using information readily available during the integration of an ode using this development a variety of effective ode solvers could be made insensitive to the type of problem i e the code would automatically recognize and alter automatically its algorithm at any step depending on whether the problem is stiff there

run time estimation of the spectral radius of jacobians it is useful for ordinary differential equation ode solvers to include an estimator of the spectral radius of the jacobian matrix of the system of ode s since this determines the numerical stability of the method hence a knowledge of spectral radius enables the run time selection of a more efficient integrator some techniques for estimating spectral radius are described and compared they include methods suitable for use with any ode solvers but which require additional computation other methods are described which are suitable with runge kutta or rosenbrock methods and which require little extra computation

on critical point theory for indefinite functionals in the presence of symmetries we consider functionals which are not bounded from above or from below even modulo compact perturbations and which exhibit certain symmetries with respect to the action of a compact lie group we develop a method which permits us to prove the existence of multiple critical points for such functionals the proofs are carried out directly in an infinite dimensional hilbert space and they are based on minimax arguments the applications given here are to hamiltonian systems of ordinary differential equations where the existence of multiple time periodic solutions is established for several classes of hamiltonians symmetry properties of these hamiltonians such as time translation invariancy or evenness are exploited introduction many variational problems which arise from physics or mathematical physics are indefinite in the sense that the functionals involved are not bounded from below or from above however some of these functionals defined in an appropriate function space are semidefinite in the sense that they are bounded from below or from above if perturbed with a weakly continuous functional this paper deals with functionals which are not semidefinite usually problems involving indefinite functionals are more difficult to handle and only very recently a method has been developed which permits us to treat such problems directly in an infinite dimensional function space 7 in 7 some theorems have been proved which establish the existence of at least one nontrivial critical point for such functionals in many physical situations there are problems which have symmetries with respect to the action of some lie group in this case we expect the existence of many critical points this has been established for semidefinite functionals cf e g 8 and 3 for even functionals 9 for a zp action with p prime numbers 15 16 10 and 6 for an s action in this paper we develop a method which allows us to estimate the number of critical points of indefinite functionals in the presence of symmetries more precisely received by the editors july 16 1980 and in revised form june 26 1981 1980 mathematics subject classification primary 58e05 70h05 secondary 34c25

an algorithm for starting multistep methods the present paper deals with the problem of starting multistep methods we take into consideration an adams bashforth moulton pece pair with the predictor of order q and the corrector of order q 1 to start this method q 1 starting values are necessary in addition to y0 a well known result from theory says that the order of convergence of the whole integration is q 1 if all those starting values are accurate of that same order present production codes start with a predictor of order 1 and a corrector of order 2 at the first step and then proceed step by step each time raising the order by 1 until all the necessary starting values have been obtained but in this manner all the starting errors keep of order 3 and so the whole integration converges no faster than that order this drawback is normally compensated for by taking very small step sizes in the starting phase the general algorithm we propose furnishes at a reasonably low cost the necessary number q 1 of starting values each of the appropriate order q 1 whatever q might be it is independent of the particular multistep formula considered and is mainly designed to be used for high values of q q 10 where the alternative strategies are too expensive or do not exist at all the numerical results reported show the validity of our approach 2003 elsevier science ltd all rights reserved

closed characteristics on compact convex hypersurfaces in r 2n for any given compact c 2 hypersurface in r 2n bounding a strictly convex set with nonempty interior in this paper an invariant n is defined and satisfies n n 2 1 where a denotes the greatest integer which is not greater than a r the following results are proved in this paper there always exist at least n geometrically distinct closed characteristics on if all the geometrically distinct closed characteristics on are nondegenerate then n n if the total number of geometrically distinct closed characteristics on is finite there exists at least an elliptic one among them and there exist at least n 1 of them possessing irrational mean indices if this total number is at most 2 n 2 there exist at least two elliptic ones among them

existence of multiple periodic orbits of hamiltonian systems on positive type hypersurfaces in r2n abstract this paper deals with the fixed energy problem of hamiltonian systems in r 2n the main result is an existence theorem of multiple periodic orbits for hamiltonian vector fields on a class of contact type but not necessarily star shaped hypersurfaces this is the first multiplicity result for hypersurfaces that are more general than star shaped ones the main result is achieved by using critical point theory for a special hamiltonian function and central gravity of the proof relies on a delicate construction of such hamiltonian function

multiplicity of closed characteristics on symmetric convex hypersurfaces in mathbb r 2n let sigma be a compact c 2 hypersurface in mathbb r 2n bounding a convex set with non empty interior in this paper it is proved that there always exist at least n geometrically distinct closed characteristics on sigma if sigma is symmetric with respect to the origin

solutions of minimal period for a class of nonconvex hamiltonian systems and applications to the fixed energy problem on considere le systeme hamiltonien jz h z ou j x y y x x y r n r n h est non convexe a un bon comportement superquadratique et satisfait une condition simple reliant ses coefficients de superquadraticite et l exposant dans ce cas on demontre que pour t 0 l existence de solutions pour le systeme hamiltonien de periode minimale t

periodic solutions on hypersurfaces and a result by c viterbo on considere un champ vectoriel hamiltonien x j h x x h x sur x r 2n h etant une fonction lisse dont le gradient h est defini par rapport a la metrique euclidienne on cherche des solutions periodiques sur une surface d energie donnee s x r 2n h x const supposee reguliere

on the number of periodic orbits of hamiltonian systems on positive type hypersurfaces in r2n abstract in this paper we prove some existence theorems of multiple periodic orbits of hamiltonian systems on a class of hypersurfaces in r 2n that are more general than star shaped ones our results extend some earlier works on this subject

convex hamiltonian energy surfaces and their periodic trajectories in this paper we introduce symplectic invariants for convex hamiltonian energy surfaces and their periodic trajectories and show that these quentities satisfy several nontrivial relations in particular we show that they can be used to prove multiplicity results for the number of periodic trajectories

formulation and solution of inverse spaghetti problem application to beam deployment dynamics a finite element model of the dynamics of axially moving highly flexible beams is presented the model is based on a geometrically nonlinear beam formulation that allows for large overall motions to account for the varying length of a deploying beam a moving finite element reference grid is incorporated within the nonlinear beam formulation such that the number of finite element nodes remains fixed and the finite element length is allowed to vary hamilton s law is used to formulate the equations of motion and a transient integration solution procedure is derived from a space time finite element discretization of the hamiltonian variational principle computational results of the methodology are presented for a planar inverse spaghetti problem he dynamics of axially moving beamlike structures are becoming increasingly important for the analysis of spacecraft and large space structures that deploy flexible appendages such as antenna stabilizing booms solar arrays and long trusslike structures as well as for other applications such as magnetic tape drives printing machines traveling cables and band saws an extensive amount of research has focused on the modeling of axially moving continua to analyze the dynamic behavior of such structures one of the most common models of axially moving continua that has received the most attention is the traveling string problem in what he termed the spaghetti problem carrier1 analyzed the motion of a string being accelerated upward into a fixed orifice in other works since then analytical expressions representing both linear and nonlinear vibrations of axially moving strings subject to various support conditions have been derived 2 9 the dynamics of beamlike structures traveling between or over fixed supports have also been studied to address industrial applications such as high speed tape drives and band saws 10 13 to address spacecraft applications in which antenna or other long flexible structures are deployed from a host satellite the dynamics of a cantilever beam that is ejected from a single fixed guide are modeled using bernoulli euler beam theory in ref 14 in a quite different approach the extrusion of a beam from a rotating base is modeled using a series of elastically connected rigid links in which the number of rigid links are continuously increased decreased to account for deployment retrieval of the beam 15 general formulations of axially moving bernoulli euler beam models cantilevered to a rotating rigid host body have also been developed to study the effect of appendage deployment on the attitude dynamics of orbiting spacecraft 16 20 all of these given analyses are based on an assumption of small elastic deformations with the exception of the discrete spring mass model of ref 15 all of the other models that have been referenced thus far use linear combinations of space depende nt admissible functions of a time varying length weighted by time dependen t generalized coordi

on noncharacteristic boundary conditions for discrete hyperbolic initial boundary value problems abstract the stability properties of some discrete hyperbolic initial boundary value problems involving two dependent variables and noncharacteristic boundary conditions are examined approximations based on the lax wendroff leap frog and crank nicolson methods are considered the most accurate approximations are found to be the least stable the results for the model problems can be used in the stability analysis of larger problems where characteristics occur singly or in pairs in particular discretizations of the shallow water and gas dynamics equations which involve boundary conditions written in natural variables are considered some computations demonstrating the stability results are discussed an appendix discussing numerical stability verification techniques is included

bounds on nonlinear operators in finite dimensional banach spaces we consider lipschitz continuous nonlinear maps in finite dimensional banach and hilbert spaces boundedness and monotonicity of the operator are characterized quantitatively in terms of certain functionals these functionals are used to assess qualitative properties such as invertibility and also enable a generalization of some well known matrix results directly to nonlinear operators closely related to the numerical range of a matrix the gerschgorin domain is introduced for nonlinear operators this point set in the complex plane is always convex and contains the spectrum of the operator s jacobian matrices finally we focus on nonlinear operators in hilbert space and hint at some generalizations of the von neumann spectral theory

space time spectral element method for solution of second order hyperbolic equations abstract a two field mixed hybrid spectral element formulation in time and space is presented explicit stsee and implicit stsei algorithms are developed and used to solve a second order scalar hyperbolic equation the algorithms were analyzed and their characteristics in terms of truncation error and stability are presented this analysis was confirmed by numerical results the improvement in accuracy in comparison with semi discrete schemes is remarkable

existence of multiple periodic orbits on star shaped hamiltonian surfaces consider the hamiltonian system hs i 1 n here h c2 2n in this paper we investigate the existence of periodic orbits of hs on a given energy surface z 2n h z c c o is a constant the surface is required to verify certain geometric assumptions bounds a star shaped compact region and for some ellipsoid 2n o o depending in an explicit fashion on the lengths of the main axes of and one other geometrical parameter of such that if furthermore 2 2 1 then hs has at least n distinct geometric orbits on this result is shown to extend and unify several earlier works on this subject among them works by weinstein rabinowitz ekeland lasry and ekeland in proving this result we construct index theories for an s1 action from which we derive abstract critical point theorems for s1 invariant functionals we also derive an estimate for the minimal period of solutions to differential equatious