Source: https://remi.flamary.com/cours/methodnumpy/04_integdif_2x2.pdf