ASKSAGE: Sage Q&A Forum - Individual question feedhttp://ask.sagemath.org/questions/Q&A Forum for SageenCopyright Sage, 2010. Some rights reserved under creative commons license.Fri, 09 Jun 2017 16:33:25 -0500indicial equationhttp://ask.sagemath.org/question/37861/indicial-equation/Let N be an integer. Let a, b, A_1,...,A_N be constant real numbers. Let g(n) be a real function of integer variable n.
g(n) satisfies the recurrence equation:
(n + 1)*(n + b)*g(n+1) = (n + a)*g(n) + sum(i=1 to N)(A_i * g(n-i))
with g(n)=0 for n<0 . g(0) is obtained through boundary conditions, so it can be considered a given constant.
I need to know the general functional form of the term g(n), as a function of the given constants of the problem: a,b,N,A_1,...,A_N and g(0).
Could you help me? Could I programm in SAGE the code to provide the searched general term g(n)?
Thanks for your attention.
Javier GarciaThu, 08 Jun 2017 12:34:46 -0500http://ask.sagemath.org/question/37861/indicial-equation/Comment by fjgg1549 for <p>Let N be an integer. Let a, b, A_1,...,A_N be constant real numbers. Let g(n) be a real function of integer variable n.</p>
<p>g(n) satisfies the recurrence equation:</p>
<p>(n + 1)<em>(n + b)</em>g(n+1) = (n + a)*g(n) + sum(i=1 to N)(A_i * g(n-i))</p>
<p>with g(n)=0 for n<0 . g(0) is obtained through boundary conditions, so it can be considered a given constant.</p>
<p>I need to know the general functional form of the term g(n), as a function of the given constants of the problem: a,b,N,A_1,...,A_N and g(0).</p>
<p>Could you help me? Could I programm in SAGE the code to provide the searched general term g(n)?</p>
<p>Thanks for your attention.</p>
<p>Javier Garcia</p>
http://ask.sagemath.org/question/37861/indicial-equation/?comment=37899#post-id-37899Thank you dan_fulea. Typically N would be <20, and the rest of the constants would not be integers. It is a pity that there is not a readily available general term. I was hoping to find out some general properties about this new function: some derivatives and simple integrals. But 100 terms in the series is not so bad: it should provide a graph good enough. I am not sure if I should dare to ask you for help. In any event, thanks for your attention.
Javier GarciaFri, 09 Jun 2017 16:33:25 -0500http://ask.sagemath.org/question/37861/indicial-equation/?comment=37899#post-id-37899Comment by dan_fulea for <p>Let N be an integer. Let a, b, A_1,...,A_N be constant real numbers. Let g(n) be a real function of integer variable n.</p>
<p>g(n) satisfies the recurrence equation:</p>
<p>(n + 1)<em>(n + b)</em>g(n+1) = (n + a)*g(n) + sum(i=1 to N)(A_i * g(n-i))</p>
<p>with g(n)=0 for n<0 . g(0) is obtained through boundary conditions, so it can be considered a given constant.</p>
<p>I need to know the general functional form of the term g(n), as a function of the given constants of the problem: a,b,N,A_1,...,A_N and g(0).</p>
<p>Could you help me? Could I programm in SAGE the code to provide the searched general term g(n)?</p>
<p>Thanks for your attention.</p>
<p>Javier Garcia</p>
http://ask.sagemath.org/question/37861/indicial-equation/?comment=37890#post-id-37890In case there is a special situation of interest, please give us the special constants.
(One or more cases of interest.)
(I hope $N$ is at most two, and $a,b$ are rational.)
(As a comment here, or as an addition to the initial post.)
There is no problem to compute the first say some $100$ terms in the sequence / coefficients in the series, although the denominators will quickly grow in maginitude, and the result may be of no use. But at least writing the code will take not so much time. If this is the question, then there is an answer. A 20 minutes job, maybe. For a special situation then an answer (for guessing the other coefficients) can maybe be given, or at least guessing the shape may have a starting point. (But for the general problem...)Fri, 09 Jun 2017 12:55:20 -0500http://ask.sagemath.org/question/37861/indicial-equation/?comment=37890#post-id-37890Comment by dan_fulea for <p>Let N be an integer. Let a, b, A_1,...,A_N be constant real numbers. Let g(n) be a real function of integer variable n.</p>
<p>g(n) satisfies the recurrence equation:</p>
<p>(n + 1)<em>(n + b)</em>g(n+1) = (n + a)*g(n) + sum(i=1 to N)(A_i * g(n-i))</p>
<p>with g(n)=0 for n<0 . g(0) is obtained through boundary conditions, so it can be considered a given constant.</p>
<p>I need to know the general functional form of the term g(n), as a function of the given constants of the problem: a,b,N,A_1,...,A_N and g(0).</p>
<p>Could you help me? Could I programm in SAGE the code to provide the searched general term g(n)?</p>
<p>Thanks for your attention.</p>
<p>Javier Garcia</p>
http://ask.sagemath.org/question/37861/indicial-equation/?comment=37872#post-id-37872(homework?)
First, one may look at the particular case when the $A$-constants are all $0$. And $g(0)=1$. Then $g(1) = \frac a{1\cdot b}$, then $g(2)= \frac {a(a+1)}{2!\cdot b(b+1)}$, and so on. Looks as in a general hypergeometric dream. So it is natural to associate the *formal* sum:
$$ F(z) = \sum_{n\ge 0} g(n)z^n $$
and to rephrase the given recursion, as an equation in $F$, that sage may eventually solve. For instance, $zF'(z) = \sum_{n\ge 1} n\,g(n)\, z^n = \sum_{n\ge 0} (n+1)\,g(n+1)\, z^{n+1} $, and $(z^bF'(z))'=\dots$
The recursion needs in fact the preceding $(N+1)$ coefficients for the one to be computed. Well, defining all negative ones as zero is ok, but then the recursion holds for $n\ge 0$. For $b\in-\mathbb N$ we have problems.) Sage can compute first coeffs, but generalThu, 08 Jun 2017 14:27:29 -0500http://ask.sagemath.org/question/37861/indicial-equation/?comment=37872#post-id-37872Answer by fjgg1549 for <p>Let N be an integer. Let a, b, A_1,...,A_N be constant real numbers. Let g(n) be a real function of integer variable n.</p>
<p>g(n) satisfies the recurrence equation:</p>
<p>(n + 1)<em>(n + b)</em>g(n+1) = (n + a)*g(n) + sum(i=1 to N)(A_i * g(n-i))</p>
<p>with g(n)=0 for n<0 . g(0) is obtained through boundary conditions, so it can be considered a given constant.</p>
<p>I need to know the general functional form of the term g(n), as a function of the given constants of the problem: a,b,N,A_1,...,A_N and g(0).</p>
<p>Could you help me? Could I programm in SAGE the code to provide the searched general term g(n)?</p>
<p>Thanks for your attention.</p>
<p>Javier Garcia</p>
http://ask.sagemath.org/question/37861/indicial-equation/?answer=37879#post-id-37879Thank you dan_fulea. Yes, you are right. My problem is related to the hypergeometric function, but it has a greater complexity. I have manually found the first 4 terms of the series, g(1) to g(4), and they look rather complicated. Do you think there is any hope that SageMath could solve this problem?
Thanks, Javier Garcia.
Thu, 08 Jun 2017 16:59:16 -0500http://ask.sagemath.org/question/37861/indicial-equation/?answer=37879#post-id-37879Comment by fjgg1549 for <p>Thank you dan_fulea. Yes, you are right. My problem is related to the hypergeometric function, but it has a greater complexity. I have manually found the first 4 terms of the series, g(1) to g(4), and they look rather complicated. Do you think there is any hope that SageMath could solve this problem? </p>
<p>Thanks, Javier Garcia.</p>
http://ask.sagemath.org/question/37861/indicial-equation/?comment=37882#post-id-37882Thank you dan_fulea. Yes, you are right. My problem is related to the hypergeometric function, but it has a greater complexity. I have manually found the first 4 terms of the series, g(1) to g(4), and they look rather complicated. Do you think there is any hope that SageMath could solve this problem?Thu, 08 Jun 2017 17:00:16 -0500http://ask.sagemath.org/question/37861/indicial-equation/?comment=37882#post-id-37882