Hello everyone! I am calculating the longitudinal resistance (SdHo oscillation) using a square lattice.The order of leads is as follows: #Left Lead sys.attach_lead(leadA)# lead no. 0 #right lead sys.attach_lead(leadA.reversed())# lead no. 1 # First upper lead sys.attach_lead(leadB)# lead no. 2 # First lower lead sys.attach_lead(leadB.reversed())# lead no. 3 # Second upper lead sys.attach_lead(leadC) # lead no. 4 #second lower lead sys.attach_lead(leadC.reversed()) # lead no. 5 And using the Folowing piece of code to calculate the Longitudinal resistance: energy =2 def plot_conductance(sys, energy, Bs): data = [] for B in Bs: smatrix = kwant.smatrix(sys, energy, args=[B]) Trans = np.array([[smatrix.transmission(i, j) for j in xrange(6)] for i in xrange(6)]) Trans -= np.diag(Trans.sum(axis=0)) cm = Trans[:-1, :-1] resistance = np.linalg.solve(cm, [-1,1,0, 0, 0])[3] data.append(resistance) My question is that it is correct also what is the unit of current being used here, and units of rest of the quantities?Also How I am supposed to write code gate voltage. I am trying the same way accept now my arguement is not B but G (gate voltage). But its wrong I know.Please help me about these questions. Thank you in advance.