function GLRLM0=get_GLRLM(vec)
% N=size(v,2);
% for i=N:-1:1
%     for j=1:1:(N-i+1)
%         vtt=v(j:1:(j+i-1));
%         mean_t=mean(vtt);
%         t=abs(vtt-mean_t);
%         y=find(t>0.5);
%         if size(y,2)==0
%             if i==N
%                 GLRLM0=[mean_t;N];
%             elseif j~=1 && (i+j-1)~=N
%                 GLRLM0=[[mean_t;N] get_GLRLM(v(1:(j-1))) get_GLRLM(v((i+j-1):end))];
%             elseif j~=1 && (i+j-1)==N
%                 GLRLM0=[[mean_t;N] get_GLRLM(v(1:(j-1)))];
%             elseif j==1 && (i+j-1)~=N
%                 GLRLM0=[[mean_t;N] get_GLRLM(v((i+j-1):end))];
%             end
%             return
%         end
%     end
% end


vv{1}=vec;
NN=size(vec,2);
GLRLM0=[];
while (NN~=0)
    vvt=[];
    n_new=0;
    n=size(vv,2);
    for k=1:1:n
        vt=vv{k};
        N=size(vt,2);
        tt=0;
        i=N;
        while tt==0 && i>=1
            j=1;
            while tt==0 && j<=(N-i+1)
                vtt=vt(j:1:(j+i-1));
                mean_t=mean(vtt);
                t=abs(vtt-mean_t);
                y=find(t>0.5);
                if size(y,2)==0
                    NN=NN-i;
                    tt=1;
                    if i==N
                        GLRLM0=[GLRLM0;mean_t,i];
                    elseif j~=1 && (i+j-1)~=N
                        GLRLM0=[GLRLM0;mean_t,i];
                        vvt{n_new+1}=vt(1:(j-1));
                        vvt{n_new+2}=vt((i+j):end);
                        n_new=n_new+2;
                    elseif j~=1 && (i+j-1)==N
                        GLRLM0=[GLRLM0;mean_t,i];
                        vvt{n_new+1}=vt(1:(j-1));
                        n_new=n_new+1;
                    elseif j==1 && (i+j-1)~=N
                        GLRLM0=[GLRLM0;mean_t,i];
                        vvt{n_new+1}=vt((i+j):end);
                        n_new=n_new+1;
                    end
                end
                j=j+1;
            end
            i=i-1;
        end
    end
    vv=vvt;
end
    