  We study the class of non-holonomic power series with integer coefficients
that reduce, modulo primes, or powers of primes, to algebraic functions. In
particular we try to determine whether the susceptibility of the square-lattice
Ising model belongs to this class, and more broadly whether the susceptibility
is a solution of a differentially algebraic equation. Initial results on
Tutte's non-linear ordinary differential equation (ODE) and other simple
quadratic non-linear ODEs suggest that a large set of differentially algebraic
power series solutions with integer coefficients might reduce to algebraic
functions modulo primes, or powers of primes. Here we give several examples of
series with integer coefficients and non-zero radius of convergence that reduce
to algebraic functions modulo (almost) every prime (or power of a prime). These
examples satisfy differentially algebraic equations with the encoding
polynomial occasionally possessing quite high degree (and thus difficult to
identify even with long series). Additionally, we have extended both the high-
and low-temperature Ising square-lattice susceptibility series to 5043
coefficients. We find that even this long series is insufficient to determine
whether it reduces to algebraic functions modulo $3$, $5$, etc. This negative
result is in contrast to the comparatively easy confirmation that the
corresponding series reduce to algebraic functions modulo powers of $2$.
