matrix := [
[21,0,0],
[1,10,0],
[1,0,10]
]*Z(31);