Eu tenho a integral tripla da figura. Eu usei a documentação do scipy
para resolver ela.
Agora eu quero trocar os valores dos limites superiores de x, y e z, que nesse problema são os valores 2, 3 e 2 por arrays.
Tipo onde for 2 (x) eu quero um numpy array [2,3,4,5]
, onde for o valor 3 (y) eu quero o np.array [5,6,7,8]
e onde for 2 (z) eu quero o np.array [0,1,2,3]
.
Como eu faço essa substituição?
Abaixo segue o meu código:
from scipy import integrate
f = lambda z, y, x: (12*x)*(y**2)*(z**3)
a = integrate.tplquad(f, -1, 2, lambda x: 0, lambda x: 3,
lambda x, y: 0, lambda x, y: 2)
b = np.array([2,3,4,5])
c = np.array([5,6,7,8])
d = np.array([0,1,2,3])
g = integrate.tplquad(f, -1, b, lambda x: 0, lambda x: c,
lambda x, y: 0, lambda x, y: d)