  We propose a numerical method for resummation of perturbative series, which
is based on the stochastic perturbative solution of Schwinger-Dyson equations.
The method stochastically estimates the coefficients of perturbative series,
and incorporates Borel resummation in a natural way. Similarly to the "worm"
algorithm, the method samples open Feynman diagrams, but with an arbitrary
number of external legs. As a test of our numerical algorithm, we study the
scale dependence of the renormalized coupling constant in a theory of
one-component scalar field with quartic interaction. We confirm the triviality
of this theory in four and five space-time dimensions, and the instability of
the trivial fixed point in three dimensions.
