|  |  D.15.30.7 ivmatGaussian Procedure from libraryrootisolation.lib(see  rootisolation_lib).
 
Example:Usage:
ivmatGaussian(A);A ivmat
Return:
0, if Anot invertible,(1, Ainv)ifAinvertible
whereAinvis the inverse matrix
Purpose:
Inverts an interval matrix using Gaussian elimination in the setting
of interval arithmetic. Pivoting is handled as a special case as
I/I != [1,1]andI-I != [0,0]in general.
 |  | LIB "rootisolation.lib";
ring R = 0,(x,y),dp;
ideal I = 2x2-xy+2y2-2,2x2-3xy+3y2-2;
box B = list(bounds(7/8, 9/8), bounds(-1/10, 1/20));
ivmat J = evalJacobianAtBox (I, B); J;
==> [69/20, 23/5],[-61/40, -27/40]
==> [67/20, 24/5],[-159/40, -93/40]
==> 
list result = ivmatGaussian(J);
ivmat Jinv = result[2];
Jinv;
==> [1060/4273, 620/187],[-1220/561, -180/4273]
==> [2680/12819, 1280/187],[-920/187, -3680/12819]
==> 
Jinv * J;
==> [-7657037/799051, 12073835/799051],[-3961862/799051, 6773506/799051]
==> [-18293234/799051, 73172936/2397153],[-7807572/799051, 15513600/799051]
==> 
 | 
 
 |