matrix := [
[0,2,0,0],
[2,0,0,0],
[0,0,2,1],
[0,0,0,1]
]*Z(3);