Let consider the family of linear systems defined by the matrix equality

where is a set of unknowns, a square matrix whose elements are functions of and a dimensional vector whose elements are also functions of . The

The enclosure can be computed using the `LinearBound` procedure
whose syntax is

LinearBound(A,B,Derivative,Vars,Init)where:

`A`is a square array whose elements are function of the unknowns in`Vars``B`is an array whose elements are function of the unknowns in`Vars``Derivative`is an integer. If set to 0 the procedure uses the classical interval Gaussian elimination scheme to calculate the enclosure. If set to 1 it will use an`ALIAS`specific version of the Gaussian elimination scheme that uses the derivatives of the elements of`A, B`to improve the enclosure calculation`Vars`: a list of unknowns names`Init`: a list of ranges for the unknowns

with(ALIAS): with(linalg): A:=array([[x,y],[x,x]]): B:=array([x,y]): VAR:=[x,y]: LinearBound(A,B,0,[x,y],[[3,4],[1,2]]); LinearBound(A,B,1,[x,y],[[3,4],[1,2]]);which returns the enclosure [[.76923076923077, 10], [-13, -.076923076923077]] if