driver file to reproduce numerical results