docs: def _plot_MLmodel(ax, sampler, modelidx, e_range, e_npoints, e_unit, sed):
      """compute and plot ML model"""
  
      ML, MLp, MLerr, ML_model = _calc_ML(
          sampler, modelidx, e_range=e_range, e_npoints=e_npoints
      )
      f_unit, sedf = sed_conversion(ML_model[0], ML_model[1].unit, sed)
  
      ax.loglog(
          ML_model[0].to(e_unit).value,
          (ML_model[1] * sedf).to(f_unit).value,
          color="k",
          lw=2,
          alpha=0.8,
      )