Look at Marcos' work in Automatic adaptive step size. See if marcos is willing to partner and investigate ways that it can be adopted in progpy