Source: http://www.lmm.jussieu.fr/~hoepffner/enseignement/matlab2013_tp2correc.pdf