擅长:python、mysql、java
<p>你可以把一个解析积分。假设我没有因为起床太晚而错过一些愚蠢的事情,下面的方法有帮助吗?在</p>
<pre><code>import sympy as sy
sys.displayhook = sy.pprint
A, y, z, tm, t, tni = sy.symbols('A, y, z, tm, t, tni')
A = 2*z*sy.exp(-2*z*y + z**2)
expr = sy.integrate(A, (z,0,t)) # patience - this takes a while
expr
# check:
(sy.diff(expr,t).simplify() - A.replace(z,t)).simplify()
# thus, the result:
expr.replace(y,tm/(2*tni)).replace(t,t/tm)
</code></pre>
<p>最后一行给出了解析形式的A函数的积分,尽管它需要计算虚误差函数(可以用它来计算scipy.special.erfi()). 在</p>