import numpy as np i=0 j=0 ki=[] old_ln=1 xy=[] m=0 ll=[] for k in range(1,100): ak=2**np.sqrt(k) ln=int(np.log10(ak))+1 n=int(ak/(10**(ln-1))) print(k,ak) if n == 1: i+=1 #print(x) if ln!=old_ln: ki.append([old_ln,m]) xy.append([old_ln,i]) ll.append(i/m) i=0 m=0 old_ln=ln else: m+=1 ki.append([old_ln,n]) #print(ki) #print(xy) #print(ll)