gfun - Maple Programming Help

Home : Support : Online Help : Mathematics : Discrete Mathematics : Combinatorics : gfun : gfun/algeqtodiffeq

gfun

 algeqtodiffeq
 compute a differential equation satisfied by an algebraic function

 Calling Sequence algeqtodiffeq(p, y(z), ini)

Parameters

 p - polynomial in y and z (or a polynomial equation) y - name; holonomic function name z - name; variable of the holonomic function y ini - (optional) set; specify computation of initial conditions for the resulting differential equation

Description

 • The algeqtodiffeq(p, y(z)) command computes a linear differential equation with polynomial coefficients verified by the function y(z).  The polynomial p defines an algebraic function, $\mathrm{RootOf}\left(p,y\right)$ in Maple terms.  The resulting equation is of order at most $\mathrm{degree}\left(p,y\right)-1$.
 • The resulting linear differential equation contains initial conditions in zero ($y\left(0\right)$, $\mathrm{D}\left(y\right)\left(0\right)$, and so on), and can thus be passed directly to dsolve. In general, $y\left(0\right)$ is a RootOf a polynomial, $\mathrm{D}\left(y\right)\left(0\right)$ a rational expression in $y\left(0\right)$, ${\mathrm{D}}^{\left(2\right)}\left(y\right)\left(0\right)$ a rational expression in $y\left(0\right)$, $\mathrm{D}\left(y\right)\left(0\right)$, and so on.
 • If initial conditions are specified using ini, the algeqtodiffeq function attempts to compute initial conditions for the resulting differential equation.
 If a particular solution of eq is selected by specifying initial terms of its power series expansion at the origin using the same syntax as that of initial conditions for dsolve, the algeqtodiffeq function returns the corresponding initial conditions together with the differential equation.

Examples

 > $\mathrm{with}\left(\mathrm{gfun}\right):$
 > $\mathrm{algeqtodiffeq}\left(y=1+z{y}^{2},y\left(z\right)\right)$
 ${1}{+}\left({-}{1}{+}{2}{}{z}\right){}{y}{}\left({z}\right){+}\left({4}{}{{z}}^{{2}}{-}{z}\right){}\left(\frac{{ⅆ}}{{ⅆ}{z}}{}{y}{}\left({z}\right)\right)$ (1)
 > $\mathrm{algeqtodiffeq}\left(56{a}^{3}+7{a}^{3}{y}^{3}-14yz,y\left(z\right),\left\{y\left(0\right)=-2\right\}\right)$
 $\left\{{-}{z}{}{y}{}\left({z}\right){+}{3}{}{{z}}^{{2}}{}\left(\frac{{ⅆ}}{{ⅆ}{z}}{}{y}{}\left({z}\right)\right){+}\left({-}{108}{}{{a}}^{{9}}{+}{2}{}{{z}}^{{3}}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{z}}^{{2}}}{}{y}{}\left({z}\right)\right){,}{y}{}\left({0}\right){=}{-}{2}{,}{\mathrm{D}}{}\left({y}\right){}\left({0}\right){=}{-}\frac{{1}}{{3}{}{{a}}^{{3}}}\right\}$ (2)

We can use algeqtodiffeq with diffeqtorec to determine fast Taylor expansions.

 > $p≔y=1+zy+z{y}^{5}$
 ${p}{:=}{y}{=}{{y}}^{{5}}{}{z}{+}{y}{}{z}{+}{1}$ (3)
 > $\mathrm{deq}≔\mathrm{algeqtodiffeq}\left(p,y\left(z\right)\right)$
 ${\mathrm{deq}}{:=}\left\{\left({-}{147840}{}{{z}}^{{3}}{+}{169920}{}{{z}}^{{2}}{-}{22320}{}{z}{+}{240}\right){}{y}{}\left({z}\right){+}\left({393216}{}{{z}}^{{8}}{-}{294912}{}{{z}}^{{7}}{-}{466944}{}{{z}}^{{6}}{+}{92160}{}{{z}}^{{5}}{+}{447360}{}{{z}}^{{4}}{+}{65664}{}{{z}}^{{3}}{-}{265488}{}{{z}}^{{2}}{+}{29064}{}{z}{-}{120}\right){}\left(\frac{{ⅆ}}{{ⅆ}{z}}{}{y}{}\left({z}\right)\right){+}\left({589824}{}{{z}}^{{9}}{-}{1130496}{}{{z}}^{{8}}{-}{30720}{}{{z}}^{{7}}{+}{878592}{}{{z}}^{{6}}{-}{132480}{}{{z}}^{{5}}{-}{444384}{}{{z}}^{{4}}{+}{203256}{}{{z}}^{{3}}{+}{67800}{}{{z}}^{{2}}{-}{1392}{}{z}\right){}\left(\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{z}}^{{2}}}{}{y}{}\left({z}\right)\right){+}\left({196608}{}{{z}}^{{10}}{-}{606208}{}{{z}}^{{9}}{+}{489472}{}{{z}}^{{8}}{+}{165888}{}{{z}}^{{7}}{+}{83200}{}{{z}}^{{6}}{+}{626112}{}{{z}}^{{5}}{+}{271848}{}{{z}}^{{4}}{+}{24488}{}{{z}}^{{3}}{-}{1408}{}{{z}}^{{2}}\right){}\left(\frac{{{ⅆ}}^{{3}}}{{ⅆ}{{z}}^{{3}}}{}{y}{}\left({z}\right)\right){+}\left({16384}{}{{z}}^{{11}}{-}{69632}{}{{z}}^{{10}}{+}{105472}{}{{z}}^{{9}}{-}{56064}{}{{z}}^{{8}}{+}{188480}{}{{z}}^{{7}}{+}{166896}{}{{z}}^{{6}}{+}{38012}{}{{z}}^{{5}}{+}{1333}{}{{z}}^{{4}}{-}{256}{}{{z}}^{{3}}\right){}\left(\frac{{{ⅆ}}^{{4}}}{{ⅆ}{{z}}^{{4}}}{}{y}{}\left({z}\right)\right){,}{y}{}\left({0}\right){=}{1}\right\}$ (4)
 > $\mathrm{rec}≔\mathrm{diffeqtorec}\left(\mathrm{deq},y\left(z\right),u\left(n\right)\right)$
 ${\mathrm{rec}}{:=}\left\{\left({16384}{}{{n}}^{{4}}{+}{98304}{}{{n}}^{{3}}{+}{180224}{}{{n}}^{{2}}{+}{98304}{}{n}\right){}{u}{}\left({n}\right){+}\left({-}{69632}{}{{n}}^{{4}}{-}{466944}{}{{n}}^{{3}}{-}{1060864}{}{{n}}^{{2}}{-}{958464}{}{n}{-}{294912}\right){}{u}{}\left({n}{+}{1}\right){+}\left({105472}{}{{n}}^{{4}}{+}{700416}{}{{n}}^{{3}}{+}{1332224}{}{{n}}^{{2}}{+}{208896}{}{n}{-}{995328}\right){}{u}{}\left({n}{+}{2}\right){+}\left({-}{56064}{}{{n}}^{{4}}{-}{170496}{}{{n}}^{{3}}{+}{1257216}{}{{n}}^{{2}}{+}{5973504}{}{n}{+}{6543360}\right){}{u}{}\left({n}{+}{3}\right){+}\left({188480}{}{{n}}^{{4}}{+}{1968000}{}{{n}}^{{3}}{+}{7213120}{}{{n}}^{{2}}{+}{11107200}{}{n}{+}{6572160}\right){}{u}{}\left({n}{+}{4}\right){+}\left({166896}{}{{n}}^{{4}}{+}{2962656}{}{{n}}^{{3}}{+}{18918576}{}{{n}}^{{2}}{+}{51195456}{}{n}{+}{49204800}\right){}{u}{}\left({n}{+}{5}\right){+}\left({38012}{}{{n}}^{{4}}{+}{956064}{}{{n}}^{{3}}{+}{8804404}{}{{n}}^{{2}}{+}{35087184}{}{n}{+}{50788512}\right){}{u}{}\left({n}{+}{6}\right){+}\left({1333}{}{{n}}^{{4}}{+}{53814}{}{{n}}^{{3}}{+}{747191}{}{{n}}^{{2}}{+}{4381134}{}{n}{+}{9313488}\right){}{u}{}\left({n}{+}{7}\right){+}\left({-}{256}{}{{n}}^{{4}}{-}{8064}{}{{n}}^{{3}}{-}{95216}{}{{n}}^{{2}}{-}{499464}{}{n}{-}{982080}\right){}{u}{}\left({n}{+}{8}\right){,}{u}{}\left({0}\right){=}{1}{,}{u}{}\left({1}\right){=}{2}{,}{u}{}\left({2}\right){=}{12}{,}{u}{}\left({3}\right){=}{112}{,}{u}{}\left({4}\right){=}{1232}{,}{u}{}\left({5}\right){=}{14832}{,}{u}{}\left({6}\right){=}{189184}{,}{u}{}\left({7}\right){=}{2512064}\right\}$ (5)
 > p_generator:=rectoproc(rec,u(n),list):
 > $\mathrm{p_generator}\left(30\right)$
 $\left[{1}{,}{2}{,}{12}{,}{112}{,}{1232}{,}{14832}{,}{189184}{,}{2512064}{,}{34358784}{,}{480745984}{,}{6848734464}{,}{99003237376}{,}{1448575666176}{,}{21411827808256}{,}{319255531155456}{,}{4796005997940736}{,}{72520546008219648}{,}{1102912584949792768}{,}{16859182461720526848}{,}{258886644574700699648}{,}{3991711460817459806208}{,}{61775021926688418365440}{,}{959229931916911121530880}{,}{14940323391360408046796800}{,}{233352506098550016631111680}{,}{3654109325605190169830359040}{,}{57356169042767373344673103872}{,}{902258550678887785413876908032}{,}{14222150544386213214581667397632}{,}{224605955315162319622246867402752}{,}{3553395680818726488887774467325952}\right]$ (6)