代码之家  ›  专栏  ›  技术社区  ›  LetzerWille

与scipy integrate.tplquad的结果不同

  •  1
  • LetzerWille  · 技术社区  · 6 年前

    wolfram alpha 用integrate.tplquad我得到-3。
    这是Matematica代码。

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

    我看不出integrate.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 回复  |  直到 6 年前
        1
  •  3
  •   Warren Weckesser    6 年前

    再次检查 tplquad docstring ; 您将看到函数的签名 f f(z, y, x)

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