Numerical methods for solving initial value problems on the Infinity Computer

New algorithms for the numerical solution of Ordinary Differential Equations (ODEs) with initial condition are proposed. They are designed for work on a new kind of a supercomputer – the Infinity Computer, – that is able to deal numerically with finite, infinite and infinitesimal numbers. Due to this fact, the Infinity Computer allows one to calculate the exact derivatives of functions using infinitesimal values of the stepsize. As a consequence, the new methods described in this paper are able to work with the exact values of the derivatives, instead of their approximations

Tails for the Einstein-Yang-Mills system

We study numerically the late-time behaviour of the coupled Einstein Yang-Mills system. We restrict ourselves to spherical symmetry and employ Bondi-like coordinates with radial compactification. Numerical results exhibit tails with exponents close to -4 at timelike infinity $i^+$ and -2 at future null infinity \Scri.Comment: 12 pages, 5 figure