function treePoint=treeMat(r,c)	
global row col
if (r==1)&&(c==1)
    treePoint=[];
elseif (r>row/2)||(c>col/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 (tp(4,1)<=row/2)&&(tp(4,2)<=col/2)
        t1=treeMat(chPoint(1,1),chPoint(1,2));
        tm1=[tm1;t1];
        t2=treeMat(chPoint(2,1),chPoint(2,2));
        tm2=[tm2;t2];
        t3=treeMat(chPoint(3,1),chPoint(3,2));
        tm3=[tm3;t3];
        t4=treeMat(chPoint(4,1),chPoint(4,2));
        tm4=[tm4;t4];
    end
    treePoint=[tp;tm1;tm2;tm3;tm4];
end