Source: https://perso.telecom-paristech.fr/bloch/ANIM/ImagesTP/tp_morphoSI241-Matlab.pdf