docs: def plot_dop(bands, int_max, dop, hund_cu, name):
      """Plot of Quasiparticle weight for N degenerate bands
         under selected doping shows transition only at half-fill
         the rest are metallic states"""
      data = ssplt.calc_z(bands, dop, np.arange(0, int_max, 0.1), hund_cu, name)
      ssplt.plot_curves_z(data, name)