Rewrote `plotypus.preprocessing.Fourier.trigonometric_coefficient_matrix`.
Now, instead of using a list comprehension passed to `numpy.array`, it
uses `numpy`'s bulk indexing to modify many elements of an array
initialized to the size `Nx2n+1`at once, making it almost 100 times faster.