File:Catagolue distinct still life scatterplot b3s23 C1.png

Script:

import os,sys
import numpy as np
import matplotlib.pyplot as plt
import statsmodels.api as sm

# read data file
datafile=open(os.path.join(sys.path[0],'xs_distinctplot.dat'),"rb")
x,y=np.loadtxt(datafile,delimiter="\t",skiprows=1,unpack=True)

fig,ax=plt.subplots()

# total number of distinct still lifes
td=(2,1,5,4,9,10,25,46,121,240,619,1353,3286,7773,19044,45759,112243,273188,672172,1646147,4051732,9971377,24619307,60823008,150613157,373188952,926068847,2299616637,5716948683,14223867298,35422864104)
plt.plot(x[:29], td[:29], label='OEIS A019473')

# generate a scatter plot
ax.scatter(x, y, label='Distinct')
ax.set_yscale('symlog', linthresh=1)
ax.set_xlabel('Population', fontsize=15)
ax.set_ylabel("Distinct", fontsize=15)
ax.set_title('Distinct still lifes by population', fontsize=20)
plt.ylim(ymin=0)
plt.xlim(xmin=0)
ax.grid(True)
fig.tight_layout()

# add legend. This needs to come after everything's been plotted.
plt.legend(loc='upper right')

# save plot to file
fig.set_size_inches(8, 8)
plt.savefig('xs_distinctplot.png', dpi=80)

Input data:

Population	Count
4	2
5	1
6	5
7	4
8	9
9	10
10	25
11	46
12	121
13	240
14	619
15	1352
16	3210
17	6833
18	13203
19	22041
20	32872
21	43481
22	52111
23	55331
24	54494
25	47691
26	40329
27	30567
28	22840
29	15194
30	10268
31	6014
32	3727
33	1983
34	1227
35	616
36	403
37	137
38	139
39	44
40	58
41	14
42	21
43	4
44	13
45	3
46	9
47	1
48	4
50	1
56	2