function [n,ne,nr,crl]=CummulatedRestrainedList3(nnj,rl)
kk1=0;
for jx=1:nnj
    if rl(jx)==0
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
n=kk1;%非約束位移數
for jx=1:nnj
    if rl(jx)==1
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
ne=kk1-n;%彈性約束位移數
for jx=1:nnj
    if rl(jx)==2
        kk1=kk1+1;
        crl(jx)=kk1;
    end
end
nr=kk1-n-ne;%彈性約束位移數