% theoretical run-times for x98

  MM_max = 10;

  K_max = 10;

  th_perm = zeros(K_max-1,1);

for K = 2:K_max

  K

  nstates( MM_max, K);

  th_perm(K-1) = M*prod([1:K]);

end
