expand the non-singular zero of a characterizable differential ideal into integral power series
power_series_solution (point, order, J, 'syst', 'params')
list or set of names or equations
characterizable differential ideal
Important: The diffalg package has been deprecated. Use the superseding package DifferentialAlgebra instead.
The function power_series_solution computes a formal integral power series solution of the differential system equations J=0, inequations J≠0. Such a system is formally integrable. See the last example below.
The parameter point furnishes the point of expansion of the formal power series. It is a set or a list of equations x=v where x is one of the derivation variables and v is its value.
If point is a singular point of equations (J), then power_series_solution returns FAIL. Nevertheless, this does not mean that no formal power series solution exists at that point.
When point is not singular, the series is truncated at the order given by the parameter order. They could be expanded up to any order, though convergence is not guaranteed.
The result is presented as a list of equations ui=si⁡x1⁢...⁢xm, where the ui are the differential indeterminates and the si are series in the derivation variables.
The series involve parameters corresponding to initial conditions to be given.
The parameters appear as _Cu, where u is a differential indeterminate if it represents the value of the solution at point, or _Cu_x, where x is some derivation variable, if it represents the value of the value of the first derivative of u according to x at point.
The parameters must satisfy a triangular system of polynomial equations and inequations given by syst in terms of the parameters involved in the power series solution.
If present, the variable params receives the subset of the parameters involved in the power series solution that can almost be chosen arbitrarily if not for some inequations in syst.
If J is a radical differential ideal represented by a list of characterizable differential ideals, the function power_series_solution is mapped on its component.
The command with(diffalg,power_series_solution) allows the use of the abbreviated form of this command.
R ≔ differential_ring⁡derivations=x,ranking=u,v
p1 ≔ ux2−v;p2 ≔ vx2−4⁢v
J ≔ Rosenfeld_Groebner⁡p1,p2,vx,R
sol ≔ power_series_solution⁡x=0,50,J1,'syst','params'
initial_cond ≔ _Cv=c2,_Cu_x=c,_Cv_x=−2⁢c,_Cu=u0
sol ≔ simplify⁡subs⁡initial_cond,sol
Let us explain now why, in general, we have to start from a characterizable differential system instead of any differential system. Consider the differential system given by these two differential polynomials.
R ≔ differential_ring⁡ranking=u,derivations=x,y:
p1 ≔ ux2−4⁢u;p2 ≔ uy−u
We are looking for a solution starting as:
It seems that we can choose an initial condition u⁡0,0=c2 (c≠0) and that, by differentiating the equations, all the coefficients in the expansion can be expressed in terms of c.
The first terms do not lead to any problem:
To compute the next term we can either differentiate p1 or p2. The problem is that the results obtained are not compatible.
The system p1=0,p2=0 is not formally integrable as it stands. The only solution of the system is:
J ≔ Rosenfeld_Groebner⁡p1,p2,R
Download Help Document