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