Vectorization of a nested for loop, where inner counter depends on outer one.

Hi, I'm trying to vectorize the following code:

 for i=3:6
  for j=i:6
  howmuch = howmuch + factorial(j) / factorial(j - i);
  end
 end

I tried the meshgrid approach, however couldn't find a way to connect "j to i". How can I do this?

Thanks!

Answer 1

[Stephen23]

>> [X,Y] = meshgrid(3:6);
>> idx = Y>=X;
>> sum(factorial(Y(idx)) ./ factorial(Y(idx)-X(idx)))
ans =
        2274

[Bidsitlov]

Thank you!