p=[];
for A=[0 1]
for B=[0 1]
for C=[0 1];
p=[p;[A B C]];
end
end
end
[r,c] = size(p);
A = [1,2,4,7;0,3,5,8;0,0,6,9;0,0,0,10];
i=1;
j=1;
B = zeros(r,4);
for i = 1:r
B(i,1)=1;
for j=2:4
if j == 2
B(i,j)=A(1+p(i,j-1),j);
elseif j ==3
B(i,j)=A(1+p(i,j-1)+p(i,j-2),j);
else
B(i,j)=A(1+p(i,j-1)+p(i,j-2)+p(i,j-3),j);
end
end
end