#pragma once
void DEM_moments(double *T_arr, double *lnT_arr, double *DEM_arr, int N, double *n_avg, double *T_avg);
void DDM_moments(double *T_arr, double *lnT_arr, double *DDM_arr, int N, double *n_avg, double *T_avg);