# hello i need base code (triu) for gauss elimination

9 views (last 30 days)
masoud mansouri on 28 Nov 2015
Answered: Geoff Hayes on 29 Nov 2015
I write body code and i have problem with code..i need good code for triu..i havent use triu code..

Geoff Hayes on 28 Nov 2015
masoud - if you have written the code and it isn't working as expected, then please attach it to your question so that we can comment on what may need to be corrected.
masoud mansouri on 28 Nov 2015
well..i wrote this but i don't like use code triu to solve matrix ...i solve gausse elimination only with loop.(my teacher said don't use triu(x) you can write only for e.g A(a11.a12)..it's my project so thanks for saw my problem ...if you skill Matlab please send me email(<mailto:masoudmansouri8@gmail.com masoudmansouri8@gmail.com>) thanks
masoud mansouri on 28 Nov 2015
clear all
clc
c = input ('please enter matrix :'); %matris elim zarib
disp(c);
b=input('please enter matrix :'); %ansewr matrix
disp(b);
b = b'
A = [c b]
...????? if you run this code after last sentence you can write triu but i don't like used them

Geoff Hayes on 29 Nov 2015
masoud - triu returns the upper triangular part of a matrix so it is should be straightforward on how to implement this function. If X is your input matrix, then the output matrix U has to be the same size as X, so
U = zeros(size(X));
Now if you are only concerned with those elements above and including the main diagonal, then you need to copy these elements from X into U. But which ones are they? Look at the diagram drawn in the description for triu. If we start with the first row, then we copy all column elements i.e. those elements from columns 1,2,3,...,*m* where m is the number of columns. When we move to the second row, we copy all those element from the second column onwards. So what is the pattern?
for r=1:size(X,1)
for c=???:size(X,2)
% do the copy from X to U
end
end
Since this is homework, I will leave it to you to determine the pattern and what c should begin iterating over.