Quantum Computing, Quantum Physics, Programming, Physics, Science.
multiply the starting # vector by it to get our result simul = scipy.linalg.expm(1j * H * (1 / (2 * np.pi))) vec = simul @ start vec = vec / vec[0] # global phase cvec = vec / np.linalg.norm(vec) # normalize Matrix(np.round(cvec, 5)) Output from numerical calculation By simply comparing the two