Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Solving Ordinary Differential Equations

A new class of boundary conditions is described for quantum systems integrable by means of the quantum inverse scattering (R-matrix) method. The method proposed allows the author to treat open quantum chains with appropriate boundary terms in the Hamiltonian. The general considerations are applied to the XXZ and XYZ models, the nonlinear Schrodinger equation and Toda chain.

https://math.berkeley.edu/~reshetik/Sklyanin-boundary.pdf

Boundary conditions for integrable quantum systems

This paper gives a review of integration algorithms for finite dimensional mechanical systems that are based on discrete variational principles. The variational technique gives a unified treatment of many symplectic schemes, including those of higher order, as well as a natural treatment of the discrete Noether theorem. The approach also allows us to include forces, dissipation and constraints in a natural way. Amongst the many specific schemes treated as examples, the Verlet, SHAKE, RATTLE, Newmark, and the symplectic partitioned Runge–Kutta schemes are presented.

/pdf/discrete-mechanics-and-variational-integrators-4qm2t5dk6u.pdf

Discrete mechanics and variational integrators

https://digital.csic.es/bitstream/10261/15682/1/blanes.pdf

The Magnus expansion and some of its applications

Some new similarity reductions of the Boussinesq equation, which arises in several physical applications including shallow water waves and also is of considerable mathematical interest because it is a soliton equation solvable by inverse scattering, are presented. These new similarity reductions, including some new reductions to the first, second, and fourth Painleve equations, cannot be obtained using the standard Lie group method for finding group‐invariant solutions of partial differential equations; they are determined using a new and direct method that involves no group theoretical techniques.

New similarity reductions of the Boussinesq equation

Closed form expressions in terms of multi-sums of products have been given in \cite{Tranclosedform, KRQ} of integrals of sine-Gordon, modified Korteweg-de Vries and potential Korteweg-de Vries maps obtained as so-called $(p,-1)$-traveling wave reductions of the corresponding partial difference equations. We prove the involutivity of these integrals with respect to recently found symplectic structures for those maps. The proof is based on explicit formulae for the Poisson brackets between multi-sums of products.

Involutivity of integrals for sine-Gordon, modified KdV and potential KdV maps

https://www.massey.ac.nz/~rmclachl/apnum00.pdf

Numerical integrators that contract volume

In this paper we develop an integrating factor matrix method to derive conditions for the existence of first integrals. We use this novel method to obtain first integrals, along with the conditions for their existence, for two- and three-dimensional Lotka–Volterra systems with constant terms. The results are compared to previous results obtained by other methods.

An integrating factor matrix method to find first integrals

In this Letter we propose a systematic approach for detecting and calculating preserved measures and integrals of a rational map. The approach is based on the use of cofactors and Discrete Darboux Polynomials and relies on the use of symbolic algebra tools. Given sufficient computing power, all rational preserved integrals can be found. 
We show, in two examples, how to use this method to detect and determine preserved measures and integrals of the considered rational maps.

Using discrete Darboux polynomials to detect and determine preserved measures and integrals of rational maps.

We generalise the concept of duality to lattice equations. We derive a novel 3 dimensional lattice equation, which is dual to the lattice AKP equation. Reductions of this equation include Rutishauser's quotient-difference (QD) algorithm, the higher analogue of the discrete time Toda (HADT) equation and its corresponding quotient-quotient-difference (QQD) system, the discrete hungry Lotka-Volterra system, discrete hungry QD, as well as the hungry forms of HADT and QQD. We provide three conservation laws, we conjecture the equation admits N-soliton solutions and that reductions have the Laurent property and vanishing algebraic entropy.

G. R. W. Quispel

Papers

Involutivity of integrals for sine-Gordon, modified KdV and potential KdV maps

Numerical integrators that contract volume

An integrating factor matrix method to find first integrals

Using discrete Darboux polynomials to detect and determine preserved measures and integrals of rational maps.

Duality for discrete integrable systems II