Local Solution of a Nonlinear ODE Using a Power Series Expansion

Consider the nonlinear ODE , subject to the initial conditions and (set by you).

The Demonstration finds approximate solutions of the ODE for near 0 in the form . You can vary the order of the expansion from 1 to 5.

The green curve is the numerical solution obtained using Mathematica's built-in function NDSolve and the approximate power series expansion is drawn in red. The approximate solution (a polynomial) is shown; it agrees very well with the numerical solution for near 0. The light blue pane defines the region where the absolute value of the difference between the two functions is less than 0.01.