function [treePoint,chPoint]=treeMat(r,c)	
global rMat cMat
if (r==1)&&(c==1)
    treePoint=[];
elseif (r>rMat/2)||(c>cMat/2)
    treePoint=[];
else
    chPoint=[2*r-1,2*c-1;2*r-1,2*c;2*r,2*c-1;2*r,2*c];
    tm1=[];tm2=[];tm3=[];tm4=[];
    if (chPoint(4,1)<=rMat/2)&&(chPoint(4,2)<=cMat/2)
        [t1,c1]=treeMat(chPoint(1,1),chPoint(1,2));
        tm1=[tm1;t1];
        [t2,c2]=treeMat(chPoint(2,1),chPoint(2,2));
        tm2=[tm2;t2];
        [t3,c3]=treeMat(chPoint(3,1),chPoint(3,2));
        tm3=[tm3;t3];
        [t4,c4]=treeMat(chPoint(4,1),chPoint(4,2));
        tm4=[tm4;t4];
    end
    treePoint=[chPoint;tm1;tm2;tm3;tm4];
end