D = ChainComplex({0: identity_matrix(ZZ, 4), 4: identity_matrix(ZZ, 30)})
H = D.homology()

H

H.get(2,0)

v = {1:5, 2:3}

v.get(1,0)

v.get(7,0)