A=zeros(7)
B=A;

A(2:4,2:4)=ones(size(A(2:4,2:4)));
B(2+1:4+1,2+1:4+1)=ones(size(A(2:4,2:4)));

est_trans2(A,B)

A(2:4,2:4)=ones(size(A(2:4,2:4)));
B(2+1:4+1,2:4)=ones(size(A(2:4,2:4)));