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