Linsys esimerkki, tapaus (b)

> with(linalg):
> A:=matrix(3,3,[1,1,-1,3,1,1,1,-1,4]);
                                  [1     1    -1]
                                  [             ]
                             A := [3     1     1]
                                  [             ]
                                  [1    -1     4]

> b:=vector([1,9,8]);
                                b := [1, 9, 8]

> Ab:=augment(A,b);
                                [1     1    -1    1]
                                [                  ]
                          Ab := [3     1     1    9]
                                [                  ]
                                [1    -1     4    8]

> gausselim(Ab);
                             [1     1    -1    1]
                             [                  ]
                             [0    -2     4    6]
                             [                  ]
                             [0     0     1    1]