lab6ex5

clear

A=load('A.txt');
b=load('b.txt');
n=size(A,1);
x=zeros(n,1);

if size(A)==size(A')
    for i=1:n
        for j=(i+1):n
            if A(i,j)~=0
                disp('stop')
                return
            end
        end
    end

    if det(A)==0
        disp('hey')
        return
    end

    if n~=size(b,1)
        return
    end

    x(1)=b(1)/A(1,1);

    for i=2:n
        sum=0;
        x(i)=(1/A(i,i));
        for j=1:(i-1)
            sum=sum+A(i,j)*x(j);
        end
        x(i)=x(i)*(b(i)-sum);
    end
end