Numerical integration of a coupled Korteweg-de Vries system