Hi,

I want to find the roots of a system of non-linear equations like this one:

```
function f!(F::Vector{Float64},
Δ::Vector{Float64},
L::Int64,
N::Int64,
g::Float64,
Γ::Float64,
Z::AbstractMatrix{Float64})
F[1:L] .= Z*Δ
for i=1:L
F[i] += g^2*N*(L-N)*Γ - 2Δ[i] + g*sum(Z[i, :])*Δ[i] - Δ[i]^2
end
F[L+1] = sum(Δ) + 2N
return F
end
```

My problem is that I wanted to use `NLsolve`

to solve this problem, but the function `f!`

has an argument `Z`

that is an array and, therefore, `NLsolve`

can’t handle it. There seems to be an issue already filed about this but there has not been any response yet in some months.

Do you know what solver I could use to find the roots of my function?