from math import comb
def plane_gromov_witten(limit):
values = {1: 1}
for d in range(2, limit + 1):
total = 0
for a in range(1, d):
b = d - a
total += (values[a] * values[b] * a * a * b
* (b * comb(3 * d - 4, 3 * a - 2)
- a * comb(3 * d - 4, 3 * a - 1)))
values[d] = total
return values
plane_gromov_witten(15)[15]The generator evaluates Formula (1) in exact integer arithmetic. The values were compared with OEIS A013587 [5] for every entry, and the entries with $d\leq 12$ were compared with the table in [2].