Correction
On tape :
A:=[[1,a,0,0],[0,1,a,0],[0,0,1,a],[0,0,0,1]]
On définit la matrice de passage P en mettant en colonne les coordonnées
des nouveaux vecteurs de base :
P:=[[1,1,1,1],[0,1,1,1],[0,0,1,1],[0,0,0,1]]
On définit la matrice de u dans la nouvelle base :
B:=normal(changebase(A,P))
On obtient :
[[1,a,0,0],[0,1,a,0],[0,0,1,a],[0,0,0,1]]
On a :
B = inv(P)*A*P
Vérifions :
On tape :
E1:=[1,0,0,0]
E2:=[1,1,0,0]
E3:=[1,1,1,0]
E4:=[1,1,1,1]
On tape :
A*E1
On obtient :
[1,0,0,0]
On tape :
A*E2
On obtient :
[1,1,0,0]
On tape :
A*E3
On obtient :
[1,1,1,0]
On tape :
A*E4
On obtient :
[1,1,1,1]