[hat_rho_ML_t, hat_nu_ML]=copulafit('t',U);
hat_theta_g_ML=copulafit('gumbel',U);
hat_theta_c_ML=copulafit('NAME_1',U);

LL_t=sum(log(copulapdf('gumbel',U,hat_rho_ML_t,hat_nu_ML)));
LL_gumbel=sum(log(copulapdf('gumbel',U,hat_theta_g_ML)));
LL_NAME_1=sum(log(copulapdf('NAME_1',U,hat_theta_c_ML)));

BIC_t=log(n)*5-2*LL_t;
BIC_gumbel=log(n)-2*LL_gumbel;
BIC_NAME_1=log(n)-2*LL_NAME_1;