ls
#gauss.m#	backsolv.m	backsolv_c.m	cmup		index.html
Oct01		backsolv0.m	cmdown		gauss.m		xplSys.m

xplSys


A

A =

     1    -3    -4    -3
    -1     6     4     1
    -2     6    10     7
    -3     3    14    15

b

b =

    -5
    16
     9
    -9

x = gauss(A,b)

A =

     1    -3    -4    -3    -5
    -1     6     4     1    16
    -2     6    10     7     9
    -3     3    14    15    -9


A =

     1    -3    -4    -3    -5
     0     3     0    -2    11
     0     0     2     1    -1
     0    -6     2     6   -24


A =

     1    -3    -4    -3    -5
     0     3     0    -2    11
     0     0     2     1    -1
     0     0     2     2    -2


A =

     1    -3    -4    -3    -5
     0     3     0    -2    11
     0     0     2     1    -1
     0     0     0     1    -1


x =

     1
     3
     0
    -1

x = gauss(A,b)

A =

     1    -3    -4    -3    -5
    -1     6     4     1    16
    -2     6    10     7     9
    -3     3    14    15    -9


A =

     1    -3    -4    -3    -5
    -1     3     0    -2    11
    -2     0     2     1    -1
    -3    -6     2     6   -24


A =

     1    -3    -4    -3    -5
    -1     3     0    -2    11
    -2     0     2     1    -1
    -3    -6     2     2    -2


A =

     1    -3    -4    -3    -5
    -1     3     0    -2    11
    -2     0     2     1    -1
    -3    -6     2     1    -1


x =

     1
     3
     0
    -1



L = [1 0 0 ; 2 1 0; -1 1 1]

L =

     1     0     0
     2     1     0
    -1     1     1

U = [2 2 -1; 0 -1 1; 0 0 3]

U =

     2     2    -1
     0    -1     1
     0     0     3

L*U

ans =

     2     2    -1
     4     3    -1
    -2    -3     5

exit