Our new method, named variational non-Gaussian exact diagonalization (NGSED), is published at Physical Review Research! It solves strongly correlated electronic systems with strong electron-phonon coupling, revealing a superconducting state in the 2D Hubbard-Holstein model.