TTH = Matrix(4,4,[1/2,1/2,1/2,0,1/2,0,0,0,0,1/2,0,0,0,0,1/2,1])
HTH = Matrix(4,4,[1/2,0,1/2,0,1/2,1/2,0,0,0,1/2,0,0,0,0,1/2,1])
for i in range(10):
a = (TTH^i * vector([1,0,0,0]))[3]
b = (HTH^i * vector([1,0,0,0]))[3]
print i, float(a), float(b)
|
|
0 0.0 0.0
1 0.0 0.0
2 0.0 0.0
3 0.125 0.125
4 0.1875 0.25
5 0.25 0.34375
6 0.3125 0.421875
7 0.3671875 0.4921875
8 0.41796875 0.5546875
9 0.46484375 0.609375
0 0.0 0.0
1 0.0 0.0
2 0.0 0.0
3 0.125 0.125
4 0.1875 0.25
5 0.25 0.34375
6 0.3125 0.421875
7 0.3671875 0.4921875
8 0.41796875 0.5546875
9 0.46484375 0.609375
|