Series solutions of ODEs?



Can Maxima produce the solution of an ODE in form of a power series? 
For example, the Riccati equation:

dy/dx=x^2+y^2, y(0)=1

has a complicated closed-form solution using Bessel functions, but
sometimes a few terms of the series is all you need.

How do I get a series solution in Maxima?

As a related question, how do I use the diffequations/ode1_riccati.mac
solver to include an initial condition?

thanks,
Alasdair