fstr = '1<2;1<3;2<4;2<5;3<6;3<7'
T=readtree(fstr)
A=randtreeA(T,m)
A0=zeros(m); A0(1,1)=1; A0(1,2) = 1;
A0 = reshape(A0,[1 1 m m]);
A(2,5,:,:) = A0;
A(3,7,:,:) = A0;
[n,n]=size(T)
P = dtreefillprob1(m,n,T,A,C)
[hh7,Ht,H] = gethhn(P,m,n)

fstr = '1<2;1<3;2<4;3<5;'
T=readtree(fstr);
A0 = A;
A = [];
A(1,2,:,:) = A0(1,2,:,:);
A(1,3,:,:) = A0(1,3,:,:);
A(2,4,:,:) = A0(2,4,:,:);
A(3,5,:,:) = A0(3,7,:,:);
[n,n]=size(T)
P = dtreefillprob1(m,n,T,A,C)
[hh5,Ht,H] = gethhn(P,m,n)

%%%%%%%%%%%%%%%%%%%%%%%%%%%

fstr = '1<2;1<3;2<4'
T=readtree(fstr);
A=randtreeA(T,m)
[n,n]=size(T)
P = dtreefillprob1(m,n,T,A,C)
[hh4,Ht,H] = gethhn(P,m,n)


fstr = '1<2;1<3'
T=readtree(fstr);
A0 = A;
A=[];
[n,n]=size(T)
A(1,2,:,:) = A0(1,2,:,:);
A(1,3,:,:) = A0(1,3,:,:);
P = dtreefillprob1(m,n,T,A,C)
[hh3,Ht,H] = gethhn(P,m,n)