Hi,
I would like to integrate (symbolically or numerically) a rational function on some polyhedral domain. As an example, my function would be:
sage: f(a,b,c) = 1/((a+b)^2 * c)
and my domain would be the unbounded polyhedron:
sage: I = [
....: [0,1,0,0], # a > 0
....: [0,0,1,0], # b > 0
....: [0,0,0,1], # c > 0
....: [-1,1,1,1], # a + b + c > 1
....: [1,-1,-1,0], # a + b < 1
....: [0,0,-1,1] # b < c
....: ]
sage: P = Polyhedron(ieqs=I)
Is it possible to do it in Sage? With some optional package? With some other software?
(Note for suspicious minds: the integral is finite.)
