i1 : (f0,f1,f2) = genericLaurentPolynomials(1,2,2)
2 2 2 2
o1 = (a x + a x + a , b x + b x x + b x + b x + b x + b , c x + c x x + c x + c x + c x + c )
2 1 1 2 0 5 1 4 1 2 2 2 3 1 1 2 0 5 1 4 1 2 2 2 3 1 1 2 0
o1 : Sequence
|
i2 : time denseResultant(f0,f1,f2); -- using Poisson formula
-- used 0.297241 seconds
|
i3 : time denseResultant(f0,f1,f2,Algorithm=>"Macaulay"); -- using Macaulay formula
-- used 0.788567 seconds
|
i4 : time (denseResultant(1,2,2)) (f0,f1,f2); -- using sparseResultant
-- used 0.594973 seconds
|
i5 : assert(o2 == o3 and o3 == o4) |
The object denseResultant is a function closure.