I am having trouble using broyden1 from scipy to solve a differential equation - it never converges.
The specific problem I am solving is exercise 6.4 from Miranda an
