J = range(1,101)
A = [ZZ(gap.eval("NumberSmallGroups("+str(j)+")")) for j in J]
s = IndexedSequence(A,J)
P = s.plot_histogram()
show(P)