a bug in charpoly
Publié : lun. mai 28, 2018 6:11 am
There is apparently some bug in the command charpoly. If I enterI get "Not invertible Error". However, the commandgives the expected result.
Code : Tout sélectionner
A:=[[0,1,0,1,1,1,0,0,0,0,0,0,1,0,0,1],[1,0,1,0,0,1,1,0,0,0,0,0,1,1,0,0],[0,1,0,1,0,0,1,1,0,0,0,0,0,1,1,0],[1,0,1,0,1,0,0,1,0,0,0,0,0,0,1,1],[1,0,0,1,0,1,0,1,1,1,0,0,0,0,0,0],[1,1,0,0,1,0,1,0,0,1,1,0,0,0,0,0],[0,1,1,0,0,1,0,1,0,0,1,1,0,0,0,0],[0,0,1,1,1,0,1,0,1,0,0,1,0,0,0,0],[0,0,0,0,1,0,0,1,0,1,0,1,1,1,0,0],[0,0,0,0,1,1,0,0,1,0,1,0,0,1,1,0],[0,0,0,0,0,1,1,0,0,1,0,1,0,0,1,1],[0,0,0,0,0,0,1,1,1,0,1,0,1,0,0,1],[1,1,0,0,0,0,0,0,1,0,0,1,0,1,0,1],[0,1,1,0,0,0,0,0,1,1,0,0,1,0,1,0],[0,0,1,1,0,0,0,0,0,1,1,0,0,1,0,1],[1,0,0,1,0,0,0,0,0,0,1,1,1,0,1,0]];
charpoly(A,x)Code : Tout sélectionner
det(x*idn(16)-A)