# Otter's constant belongs to the unlabelled rooted, free, and planted classes.
# A quick check uses the recurrence for A000081 and removes the leading
# n^(-3/2) ratio correction.
from sage.arith.misc import divisors
from sage.rings.rational_field import QQ
N = 200
a = [0] * (N + 1)
a[1] = 1
for n in range(2, N + 1):
total = 0
for k in range(1, n):
sigma = sum(d * a[d] for d in divisors(k))
total += sigma * a[n - k]
a[n] = total // (n - 1)
float(a[N] / a[N - 1] / (1 - QQ(3) / (2 * N)))The tree-species asymptotic constants used here are tabulated by Harary, Robinson and Schwenk [1]. The two integer entries are exact, and the labelled entries are the exact constant $e$ stored as a real. The other decimal entries are transcribed from the OEIS constant pages cited in the comments.
They were independently checked against the defining sequences or radii: A000081 for Otter's constant, A004111 and A000220 for the identity-tree constant, A001678, A001679 and A000014 for the series-reduced constant, A001190 [11] for the weakly binary constant, and A000598 with A261340 [13] for the rooted ternary constant. These checks support the displayed digits but are not interval proofs of every transcribed decimal.