I want to make this kind of row vector:
I have p=(p^(0),p^(1),p^(2),...) %^(i) is not rank, its just index. p^(i)
p^(i)=(p_(4i)^(i),p_(4i+1)^(i),p_(4i+2)^(i),p_(4i+3)^(i)), where
p_j^(i)=(p_(j,0,1)^(i),p_(j,0,2)^(i),...,p_(j,0,N)^(i),
p_(j,1,1)^(i),p_(j,1,2)^(i),...,p_(j,1,N)^(i),...,
p_(j,M,1)^(i),p_(j,M,2)^(i),...,p_(j,M,N)^(i)) i,j>=0 , 4i<=j<=4i+3
for example if i have N=3 M=4, i get p^(0)=(p_(0)^(0),p_(1)^(0),p_(2)^(0),p_(3)^(0))
inside p_(0)^(0)=(p_(0,0,1)^(0),p_(0,0,2)^(0),p_(0,0,3)^(0),
p_(0,1,1)^(0),p_(0,1,2)^(0),p_(0,1,3)^(0),
p_(0,2,1)^(0),p_(0,2,2)^(0),p_(0,2,3)^(0),
p_(0,3,1)^(0),p_(0,3,2)^(0),p_(0,3,3)^(0)),
p_(0,4,1)^(0),p_(0,4,2)^(0),p_(0,4,3)^(0)).
I have this code,
p=sym('p',[1 N]) %but it can only get 1 row vector of N column.
