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\ne 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 $\mathrm{ui}=\mathrm{si}\left(\mathrm{x1}\mathrm{...}\mathrm{xm}\right)$, where the $\mathrm{ui}$ are the differential indeterminates and the $\mathrm{si}$ are series in the derivation variables.

The series involve parameters corresponding to initial conditions to be given.

The parameters appear as  $\mathrm{\_C}$u, 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.

$\mathrm{with}\left(\mathrm{diffalg}\right)\:$

$R≔\mathrm{differential\_ring}\left(\mathrm{derivations}=\left[x\right]\,\mathrm{ranking}=\left[u\,v\right]\right)$

$R≔\mathrm{ODE\_ring}$

$\mathrm{p1}≔{u\left[x\right]}^2-v\left[\right]\;$$\mathrm{p2}≔{v\left[x\right]}^2-4v\left[\right]$

$\mathrm{p1}≔u_x^2-v\left[\right]$

$\mathrm{p2}≔v_x^2-4v\left[\right]$

$J≔\mathrm{Rosenfeld\_Groebner}\left(\left[\mathrm{p1}\,\mathrm{p2}\right]\,\left[v\left[x\right]\right]\,R\right)$

$J≔\left[\mathrm{characterizable}\right]$

$\mathrm{sol}≔\mathrm{power\_series\_solution}\left(\left[x=0\right]\,50\,J\left[1\right]\,'\mathrm{syst}'\,'\mathrm{params}'\right)$

$\mathrm{sol}≔\left[u\left(x\right)=\mathrm{\_Cu}+x\mathrm{\_Cu\_x}+\frac{x^2\mathrm{\_Cv\_x}}{4\mathrm{\_Cu\_x}}\,v\left(x\right)=x\mathrm{\_Cv\_x}+x^2+\mathrm{\_Cv}\right]$

$\mathrm{syst}\;$$\mathrm{params}$

$\left[{\mathrm{\_Cu\_x}}^2-\mathrm{\_Cv}=0\,{\mathrm{\_Cv\_x}}^2-4\mathrm{\_Cv}=0\,\mathrm{\_Cu\_x}\ne 0\,\mathrm{\_Cv\_x}\ne 0\right]$

$\left\{\mathrm{\_Cv}\right\}$

$\mathrm{initial\_cond}≔\mathrm{\_Cv}=c^2,\mathrm{\_Cu\_x}=c,\mathrm{\_Cv\_x}=-2c,\mathrm{\_Cu}=\mathrm{u0}$

$\mathrm{initial\_cond}≔\mathrm{\_Cv}=c^2,\mathrm{\_Cu\_x}=c,\mathrm{\_Cv\_x}=-2c,\mathrm{\_Cu}=\mathrm{u0}$

$\mathrm{subs}\left(\mathrm{initial\_cond}\,\mathrm{syst}\right)$

$\left[0=0\,0=0\,c\ne 0\,-2c\ne 0\right]$

$\mathrm{sol}≔\mathrm{simplify}\left(\mathrm{subs}\left(\mathrm{initial\_cond}\,\mathrm{sol}\right)\right)$

$\mathrm{sol}≔\left[u\left(x\right)=\mathrm{u0}+xc-\frac{1}{2}{x}^2\,v\left(x\right)={\left(c-x\right)}^2\right]$

$R≔\mathrm{differential\_ring}\left(\mathrm{ranking}=\left[u\right]\,\mathrm{derivations}=\left[x\,y\right]\right)\:$

$\mathrm{p1}≔{u\left[x\right]}^2-4u\left[\right]\;$$\mathrm{p2}≔u\left[y\right]-u\left[\right]$

$\mathrm{p1}≔u_x^2-4u\left[\right]$

$\mathrm{p2}≔u_y-u\left[\right]$

$u\left(x\,y\right)=u\left(0\,0\right)+xu\left[x\right]\left(0\,0\right)+yu\left[y\right]\left(0\,0\right)+\frac{x^{2}u\left[x,x\right]\left(0\,0\right)}{2}+xyu\left[x,y\right]\left(0\,0\right)$

$u\left(x\,y\right)=u\left(0\,0\right)+x{u}_x\left(0\,0\right)+y{u}_y\left(0\,0\right)+\frac{x^2{u}_{x,x}\left(0\,0\right)}{2}+xy{u}_{x,y}\left(0\,0\right)$

$\mathrm{solve}\left(\mathrm{subs}\left(u\left[\right]=c^2\,u\left[x\right]=u\left[x\right]\left(0\,0\right)\,u\left[y\right]=u\left[y\right]\left(0\,0\right)\,u\left[x,x\right]=u\left[x,x\right]\left(0\,0\right)\,\left\{\mathrm{p1}\,\mathrm{p2}\,\mathrm{differentiate}\left(\mathrm{p1}\,x\,R\right)\right\}\right)\,\left\{u\left[x\right]\left(0\,0\right)\,u\left[y\right]\left(0\,0\right)\,u\left[x,x\right]\left(0\,0\right)\right\}\right)$

$\left\{u_x\left(0\,0\right)=2c\,u_y\left(0\,0\right)=c^2\,u_{x,x}\left(0\,0\right)=2\right\},\left\{u_x\left(0\,0\right)=-2c\,u_y\left(0\,0\right)=c^2\,u_{x,x}\left(0\,0\right)=2\right\}$

$\mathrm{solve}\left(\mathrm{subs}\left(u\left[x,y\right]=u\left[x,y\right]\left(0\,0\right)\,u\left[x\right]=2c\,u\left[y\right]=c^2\,\mathrm{differentiate}\left(\mathrm{p1}\,y\,R\right)\right)\,\left\{u\left[x,y\right]\left(0\,0\right)\right\}\right)$

$\left\{u_{x,y}\left(0\,0\right)=c\right\}$

$\mathrm{solve}\left(\mathrm{subs}\left(u\left[x,y\right]=u\left[x,y\right]\left(0\,0\right)\,u\left[x\right]=2c\,u\left[y\right]=c^2\,\mathrm{differentiate}\left(\mathrm{p2}\,x\,R\right)\right)\,\left\{u\left[x,y\right]\left(0\,0\right)\right\}\right)$

$\left\{u_{x,y}\left(0\,0\right)=2c\right\}$

$J≔\mathrm{Rosenfeld\_Groebner}\left(\left[\mathrm{p1}\,\mathrm{p2}\right]\,R\right)$

$J≔\left[\mathrm{characterizable}\right]$

$\mathrm{rewrite\_rules}\left(J\right)$

$\left[\left[u\left[\right]=0\right]\right]$

$\mathrm{power\_series\_solution}\left(\left[x\,y\right]\,1000\,J\left[1\right]\right)$

$\left[u\left(x\,y\right)=0\right]$