Multivariate Rootfinding: problem with NLsolve (dual numbers not applicable for a function I am using)

[nilshg]

Also consider using one of my all-time favourite packages GitHub - JuliaApproximation/FastGaussQuadrature.jl: Julia package for Gaussian quadrature