#calculating firing rate proportional to the membrane potential freq = interp(pot[l][m], [-1.069,2.781], [1,20]) #print(l,m,pot[l][m], freq) # print freq Could anyone please elaborate more about the spike_train.py and recep_field.py code. I could not get it properly .