Demonstrates how to find the ground state of the 1D Hubbard model with infinite boundary conditions using Bethe ansatz. This python implementation is based on the original paper by Lieb and Wu (2003). The bottom left plot resembles that of Knizia and Chan (2012), Fig. 1(a).
