// remplir un chemin dans une matrice 
// M_[i,sigma(i)]=q(i)
// Il faut utiliser ce qui suit 

-->diag([1,2,3,4,5])
 ans  =
 
!   1.    0.    0.    0.    0. !
!   0.    2.    0.    0.    0. !
!   0.    0.    3.    0.    0. !
!   0.    0.    0.    4.    0. !
!   0.    0.    0.    0.    5. !
 
-->M=ans;
 
-->M(:,[3,5,1,2,4])=diag([1,2,3,4,5])
 M  =
 
!   0.    0.    1.    0.    0. !
!   0.    0.    0.    0.    2. !
!   3.    0.    0.    0.    0. !
!   0.    4.    0.    0.    0. !
!   0.    0.    0.    5.    0. !