about summary refs log tree commit diff
path: root/usth/MATH2.4/labwork/3/subst.m
blob: 5c6e783fa7ebf5f4305cb860b370bffca3cef7cd (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
function x = subst (aug)
  [m n] = size (aug);
  x = zeros (m, 1);
  if (istril (aug))
    x(m) = aug(1, n) / aug(1, m);
    for k = 2 : m
      x(k) = (aug(k, n) - aug(k, k+1:m)*x(1:k-1)) / aug(k, k);
    end
  elseif (istriu (aug))
    x(m) = aug(m, n) / aug(m, m);
    for k = m - 1 : -1 : 1
      x(k) = (aug(k, n) - aug(k, k+1:m)*x(k+1:m)) / aug(k, k);
    end
  else
    error ('aug must be a triangular matrix');
  end
end