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

5 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..
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.