# Tag Info

I believe newton_krylov may be what you are looking for. The documentation is confusing but I think you can do what you want. There doesn't seem to be any way to give an analytic Jacobian for use in a Krylov solver in this software, but if you give it a function that takes in $x$ and returns a sparse matrix in the method argument then it looks like that ...