Using the Lehman.jl package to solve the large-N self-consistent equations for the Kondo (or Anderson) model, as seen in the papers https://doi.org/10.1103/PhysRevB.103.205147 or https://doi.org/10.1103/PhysRevResearch.4.L042011, or in general with the papers by Prof. Coleman's group at Rutgers (with Lebanon, Rech, Parcollet, Komijani, Drouin-Touchette...)