function crit = dist_double_power_criterion(k1,s1,f1,s2,f2)

c1 = log(s1)*log(1 - f2.^k1);

c2 = log(s2)*log(1 - f1.^k1);

crit = (c1 - c2)^2;