擅长:python、mysql、java
<p>如果这个问题没有解析解(见下文),我会考虑采用线性代数方法<a href="http://stackoverflow.com/a/22699353/8146556">here</a>:</p>
<pre><code>from sympy import Function, rsolve
from sympy.abc import n
y = Function("y")
initial_conditions = {y(0): 1, y(1): 1, y(2): 1, y(3): 1}
f = -2*y(x-1)+11*y(x-2)+12*y(x-3)-36*y(x-4)+41**(x-4)+3-y(x)
print(rsolve(f, y(x)))
print(rsolve(f, y(x), initial_conditions))
</code></pre>
<p>退货:</p>
^{pr2}$