3
gsl_rng *r = apop_rng_alloc(2312311);
4
int empirical_size = 5e3;
5
apop_model *expo = apop_model_set_parameters(apop_exponential, 1.7);
6
assert (apop_kl_divergence(expo, expo) < 1e-4);
7
apop_data *empirical = apop_data_alloc(empirical_size, 1);
8
for (int i=0; i<empirical_size; i++)
9
apop_draw(apop_data_ptr(empirical, i, 0), r, expo);
10
apop_model *pmf = apop_estimate(empirical, apop_pmf);
11
assert(apop_kl_divergence(pmf,expo) < 1e-4);
12
apop_data_free(empirical);