Source: https://maths.ac-noumea.nc/IMG/pdf/ScilabMupad2-2-1.pdf