function [V]=matV(tp,M,C)

V=zeros(tp,size(C,2));

for i=1:tp,
  V(i,:)=C*M^(i-1);
end;