function I=mtaylorind(m,P)
  
I=cell(1,P-1);
rows = ones(1,m);

for p=2:P
  rows = fliplr(cumsum(fliplr(rows)));
  M=rows(1);
  for i=1:m
     I{p-1}(i)=M-rows(i);
  end
  I{p-1}(m+1)=M;
end

   



