fact: Definition aep_bound epsilon := (aep_sigma2 P / epsilon ^ 3)%R.
imports: mathcomp(all_ssreflect ssralg ssrnum matrix) Reals mathcomp(Rstruct) ssrR Reals_ext realType_ext ssr_ext bigop_ext ssralg_ext logb fdist proba entropy