другой результат с scipy integrate.tplquad - PullRequest
0 голосов
/ 09 октября 2018

С Matematica этот интеграл дает мне 3. wolfram alpha С integrate.tplquad я получаю -3.
Это код Matematica.

Integrate[1-x, {(0,2),(0,3-1.5*x),(0,6-3*x -2*y)}]   

Я не вижу, что я делаю неправильно с ntegrate.tplquad

f = lambda x,y,z: 1-x   

x1, x2 = 0,2

y1, y2 = lambda x : 0 , lambda x:3-1.5*x

z1, z2  = lambda x,y:0, lambda x,y: 6 -3*x -2*y

print(integrate.tplquad(f,x1,x2, y1, y2, z1, z2)[0])
-3.0

1 Ответ

0 голосов
/ 09 октября 2018

Дважды проверьте tplquad строка документа ;вы увидите, что подпись функции f равна f(z, y, x).Похоже, что ваша интеграция должна быть:

   f = lambda x, y, z: 1 - z
...