Making matrix whose submatrices are from for loop

I'm trying to construct matrix whose submatrices are obtained in for loop. for loop is

for j = 2:Jmax
    P(j-1)=(sqrt(-2*p(j-1)+1)/sqrt(-2*p(j-1)))*(I-(p(j)+p(j-1))*inv(A+p(j)*I))
    V(j)=P(j-1)*V(j-1);  
    W=[W V(j)]

Is this ok? Also where to put V(1), can it be outside for loop?

[James Tursa]

What are the dimensions of the variables? In particular, is P(j-1) expected to be a scalar or a matrix?

[Ljix]

P(j-1) is a matrix, and V(j) is a vector.

Answer 1

[Walter Roberson]

for j = 2:Jmax
  P{j-1}=(sqrt(-2*p(j-1)+1)/sqrt(-2*p(j-1)))*(I-(p(j)+p(j-1))*inv(A+p(j)*I))
  V{j}=P{j-1}*V{j-1};  
  W=[W V{j}]

It is not clear whether p (lower-case) is intended to be a matrix or a vector. If it is intended to be a matrix then

for j = 2:Jmax
  P{j-1}=(sqrt(-2*p{j-1}+1)/sqrt(-2*p{j-1}))*(I-(p{j}+p{j-1})*inv(A+p{j}*I))
  V{j}=P{j-1}*V{j-1};  
  W=[W V{j}]

Based upon the p(j)*I, my suspicion would be that p(j) is intended to be a scalar, but that is not certain.