Source: https://stoltz.perso.math.cnrs.fr/Enseignements/TPs-Matlab/TP1.pdf