Я хочу построить график интегрирования функции с интервалами вокруг позиции x.
Я пробовал это с scipy.integrate.quad и scipy.integrate.cumtrapz, но ни один из них, похоже, не работал. Я думал, что это довольно распространенная задача, но не смог найти ни одного примера кода, который мог бы мне помочь. Приведенный ниже код — это попытка, которую я пытался выполнить, но она не вернула полезного результата.
Код: Выделить всё
import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate
x = np.linspace(0, 3, num=1000)
def f(v):
f=np.cos(np.pi*v)+1
return f
def y_int(v):
y_int = scipy.integrate.cumtrapz(f(v), x, initial=0)
return y_int
a=0.5
plt.plot(x, y_int(x+a)-y_int(x-a), 'r-')
plt.plot(x, f(x), 'b-')
plt.show()
Источник: https://stackoverflow.com/questions/561 ... ition-of-x
