Tweaking the epsilon value for gradient estimation lets Ellison (2007).
Now terminates with reasonable output on a substantial dataset.