A Python implementation of an algorithm to obtain a phase diagram of two-dimensional interacting lattice bosons based on cluster meanfield theory in physics.