Abstract. According to the Göttsche conjecture (now a theorem), the degree N d,δ of the Severi variety of plane curves of degree d with δ nodes is given by a polynomial in d, provided d is large enough. These "node polynomials" N δ (d) were determined by Vainsencher and Kleiman-Piene for δ ≤ 6 and δ ≤ 8, respectively. Building on ideas of Fomin and Mikhalkin, we develop an explicit algorithm for computing all node polynomials, and use it to compute N δ (d) for δ ≤ 14. Furthermore, we improve the threshold of polynomiality and verify Göttsche's conjecture on the optimal threshold up to δ ≤ 14. We also determine the first 9 coefficients of N δ (d), for general δ, settling and extending a 1994 conjecture of Di Francesco and Itzykson.