 LIB "realizationMatroids.lib";
ring r = 0,(x(1..4)),dp;
ideal I = x(1)+x(2)+x(3)+x(4);
list C = list(intvec(2,2,0,0),intvec(0,0,2,2));
//C represents the tropical fan curve which consists of the cones
//cone([(1,1,0,0)]) and cone([(1,1,0,0)]), both with weight 2
realizationDimPoly(I,C);
==> 0 x(1)^2+2*x(1)*x(2)+x(2)^2
C = list(intvec(0,0,0,4),intvec(0,1,3,0),intvec(1,0,1,0),intvec(0,2,0,0),intvec(3,1,0,0));
//C represents the tropical fan curve which consists of the cones
//cone([(0,0,0,1)]) with weight 4,
//cone([(0,1,3,0)]), cone([(1,0,1,0)]) both with weight 1,
//cone([(0,1,0,0)]) with weight 2, and
//cone([(3,1,0,0)]) with weight 1
realizationDimPoly(I,C);
==> 7 x(1)*x(2)^3+x(1)^3*x(3)+x(2)^3*x(3)+x(1)*x(3)^3
