trojGor

function [X] = troj_gor(X,b)
[n,~]=size(X);

for i=1:n-1
    m=X(i+1:n,i)/X(i,i);
    X(i+1:n,:)=X(i+1:n,:)-m*(X(i,:));
    b(i+1:n,:)=b(i+1:n,:)-m*(b(i,:));
end
for i=1:n
X(i,:)=X(i,:)/X(i,i);
b(i,1)=b(i,1)/X(i,i);
end
for i = n-1:-1:1
       X(i,:)=X(i,:)-X(n,:)*X(i,n);
       b(i,1)=b(i,1)-b(n,1)*X(i,n);

end
for i
end