Numerical solution of gravitational dynamics à la Chesler & Yaffe