A small Python code that implements a way to do a very specific integral transform — the Hankel transform.
This transform is particularly common in astrophysics and cosmology (for instance, it is required to transform a correlation function into a power spectrum and vice versa). It is also a bit of a tricky integral to do, because it is highly oscillatory. This code implements a neat idea (not mine!) that is able to accurately perform this integral with good computational performance (in most cases).