I’m currently working on a code which requires the calculation of a
Fourier series more than a million times. Fourier series are
calculated using the same function below.