擅长:python、mysql、java
<p><a href="https://math.stackexchange.com/questions/2063507/solving-this-integral-involving-ei-function">https://math.stackexchange.com/questions/2063507/solving-this-integral-involving-ei-function</a>中的答案提供了一个公式,您要对一系列的<em>t</em>值进行计算。在构建这样的脚本时,我建议您采用自下而上的方法。以下是我目前为止所写的。在</p>
<pre><code>from numpy import exp
import matplotlib.pyplot as plt
K = 1000
r = 2.5
N_0 = 10
N_t = 50
def Ei(x):
return x
def Ei_inv(x):
return 1/x
def N(t):
return K*Ei_inv(exp(r)*t+Ei(r*N_0/K))/r
t = [_ for _ in range(40)]
N = [N(_) for _ in t]
plt.plot(t,N)
plt.show()
</code></pre>
<ul>
<li>我知道一些常量是什么;因此我将它们放入代码中。在</li>
<li>我不知道什么是Ei或它的逆函数,也不知道在哪里可以找到这些函数;作为临时措施,我编写了“虚拟”函数来代替它们。在</li>
<li>现在,我可以编写一个函数,实际上我想要计算的函数是另一个具有单个参数的Python函数,即伪函数。在这一点上,我尝试使用<code>print (N(50))</code>来验证代码的各个部分在一起的行为是否有礼貌。在</li>
<li>我添加了两行,一行定义绘图的<em>x</em>或时间值,第二行是<em>y</em>或<em>N</em>值。在</li>
</ul>
<p>轮到你了,用计算Ei和它的逆函数所需要的代替伪函数。在</p>