Python使用集成绘制函数

import numpy as np

def g(t):

g_value= integrate.quad (lambda t: t, 0, t)
return (g_value)

def f(t):

f_value = t**3 - g(t)
return (f_value)

t1 = np.arange (-5, 5, 0.1)
ft1 = []
for tt in t1:

ft1.append(f(tt)[0])

plt.figure(1)
plt.subplot(211)
plt.plot(t1, ft1)
plt.show()
矢量化版本可以像这样完成
import numpy as np
def g(t):

g_value= integrate.quad (lambda t: t ,0,t)
return (g_value)

g_vectorized = np.vectorize(g)

def f(t):

f_value = t**3 - g_vectorized(t)
return (f_value)

t1 = np.arange (-5, 5, 0.1)

plt.figure(1)
plt.subplot(211)
plt.plot(t1, f(t1)[0])
plt.show()