Kunfu Peng

Running `num_harmonics()` in `fundamental_set.py`, we get all results correct, e.g., `num_harmonics(3 + 1, 34 - 1) = 1156, num_harmonics(4 + 1, 20 - 1) = 2870, ...`, except one, which...