The function call pdsolve(DE, series) returns an integral (exponents are nonnegative integers), formal (the convergence issue is not addressed), power series solution of the differential ideal defined by DE, truncated at order = n.

The initial values for the unknowns, and/or their derivatives, and/or an expansion point, can be provided in DE by enclosing all the arguments into a set or a list. The set of initial values must have only one expansion point and must be linear in the unknowns and derivatives evaluated at the expansion point. Typically, these initial values are specified as equations of the form , where the  are dependent variables and/or their derivatives evaluated at the expansion point - say  or  when the unknown is  - and  represent their values, possibly being expressions involving independent variables (in which case keep the independent variable as is, in the left-hand-side ).

If initial values are not specified, the expansion point is the origin, all the independent variables equal to , and all the functions and their derivatives different from zero. It is also possible to specify only an expansion point different than zero directly in DE - say including equations of the form .

The values attributed to the unknowns and their derivatives at the expansion point, also cannot be chosen freely. They must satisfy the system , returned by the DifferentialAlgebra command PowerSeriesSolution when using the option conditions. These conditions, which are checked by default, can be disregarded using the option disregard = ...; see the Examples section.

If the attributed values are solutions of the system , then, the returned power series are truncated solutions of the differential equation system DE.

If these values do not satisfy , then, the returned power series are not truncated solutions of the differential ideal.

If they satisfy , but do not satisfy , then, the returned power series may, or may not, be truncated solutions of the of DE. If they are, they do not need to be unique.

The opts arguments may contain one or more of the following options.

disregard = <nothing (default), inequations, equations, all>. This option permits to restrict the test for the initialvalues (see above).

equations: the equations are disregarded; pdsolve only checks that the inequations are satisfied by initialvalues.

inequations: the inequations are disregarded; pdsolve only checks that the equations are satisfied by initialvalues.

all: equations and inequations are disregarded.

order = n. This option specifies the truncation order for the series expansion; if not given, the value of n used is that of Order.