Problem with NLSolve: cant find x st f(x) =0

Specific Domains
3

Does this work for you?

julia> theta  = nlsolve(x -> f_aux(max.(0.0, x)), x0).zero
1-element Vector{Float64}:
 0.06408989836281948