LIB "interval.lib";
ring R = 0,(x,y,z),lp;
ideal I = xyz3+z2y2+x,x4+y3+2z+3,xyz+1/2;
interval J = bounds2(1/2, 3/2);
box B = list(J,J,J);
ivmat A = evalJacobianAtBox(I, B); A;
==> [17/16, 97/16],[5/16, 189/16],[7/16, 351/16]
==> [1/2, 27/2],[3/4, 27/4],[2, 2]
==> [1/4, 9/4],[1/4, 9/4],[1/4, 9/4]
==> 
A*B;
==> [29/32, 1911/32],[13/8, 267/8],[3/8, 81/8]
==> 
unitMatrix2(3) * B;
==> [1/2, 3/2],[1/2, 3/2],[1/2, 3/2]
==> 
diagMatrix(3, bounds2(0, 1)) * B;
==> [0, 3/2],[0, 3/2],[0, 3/2]
==>