from PyNormaliz import *
C = Cone( inequalities = [
[1,0,0,0,0],
[0,1,0,0,0],
[0,0,1,0,0],
[0,0,0,1,0],
[-1,0,0,0,1],
[0,-1,0,0,1],
[0,0,-1,0,1],
[0,0,0,-1,1]],
equations =
[[-1,-1,-1,-1,2 ]],
polynomial =
"x[1]*x[2]*x[3]*x[4]*x[1]*x[2]*x[3]*x[4]*(x[1]-x[2])*(x[1]-x[2])*(x[1]-x[3])*(x[1]-x[3])*(x[1]-x[4])*(x[1]-x[4])*(x[2]-x[3])*(x[2]-x[3])*(x[2]-x[4])*(x[2]-x[4])*(x[3]-x[4])*(x[3]-x[4])" )