Fourier spectral method for KdV: two solitons pass through each other

Afterwards you can solve a periodic PDE with numpy.fft and solve_ivp, treat its stiff term exactly with an integrating factor, and dealias it by the 2/3 rule.