File:Catagolue distinct still life scatterplot b3s23 SS.png

Script:

#!/usr/bin/python3

import numpy as np
import matplotlib.pyplot as plt
import statsmodels.api as sm

# read data file
datafile=open('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, 4051726]
plt.plot(x[0:21], td, label='OEIS A019473')

# generate a scatter plot
ax.scatter(x, y, label='Distinct')
ax.set_yscale('symlog', linthreshy=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=4)
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	43
12	72
13	99
14	165
15	162
16	172
17	105
18	110
19	81
20	66
21	42
22	46
23	19
24	8
25	4
26	2
27	1
28	2
30	2
31	1
32	1