diff options
author | Nguyễn Gia Phong <vn.mcsinyx@gmail.com> | 2020-02-16 14:26:55 +0700 |
---|---|---|
committer | Nguyễn Gia Phong <vn.mcsinyx@gmail.com> | 2020-02-16 14:26:55 +0700 |
commit | 82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7 (patch) | |
tree | f7b7ae0bce69070c47a1b31a85bd2bc69dfecf09 /usth/MATH2.2/labwork/3/subst.m | |
parent | 29d1001e2e21eff289bff23412e284c8b2e44595 (diff) | |
download | cp-82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7.tar.gz |
[usth] Numerical Method is MATH2.4
Diffstat (limited to 'usth/MATH2.2/labwork/3/subst.m')
-rw-r--r-- | usth/MATH2.2/labwork/3/subst.m | 17 |
1 files changed, 0 insertions, 17 deletions
diff --git a/usth/MATH2.2/labwork/3/subst.m b/usth/MATH2.2/labwork/3/subst.m deleted file mode 100644 index 5c6e783..0000000 --- a/usth/MATH2.2/labwork/3/subst.m +++ /dev/null @@ -1,17 +0,0 @@ -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 |