p1 = 0.75;
p0 = 1 - p1;
pp = [p1; p0];

q1 = 0.35;
q0 = 1 - q1;
qq = [q1; q0];

ce = hbits(pp, pp);
h2 = h2bits(p1);
testclose(ce, h2);

ce = hbits(pp, qq);
h2 = h2bits(p1, q1);
testclose(ce, h2);

test_exit