%% GAUSS elimination 

A = input("Enter matrix") ;
B = input("Enter source vector");
N = length(B);
X = zeros(N,1);
Aug = [A B];

for j = 1:N-1
    for i = j+1:N
        m = Aug(i,j)/Aug(j,j);
        Aug(i,:) = Aug(i,:) - m*Aug(j,:);
    end
end

Aug 

X(N) = Aug(N, N+1)/Aug(N,N) ;

for k = N-1:-1:1
    X(k) = (Aug(k , N+1) - Aug(k, k+1:N)*X(k+1:N))/ Aug(k,k);
end 

X


%% GAUSS SEIDEL 

A = input("Enter matrix") ;
B = input("Enter source vector");
N = length(B);
X = zeros(N,1);
n = 10 ;
Y = zeros(N,1);
P = input('initial giess vector ');
e = 0.0001;


for j = 1:n
    for i = 1:N
        X(i) = (B(i)/A(i,i)) - (A(i, [1:i-1, i+1:N]) * P([1:i-1, i+1:N])) / A(i,i);
        P(i) = X(i);
    end 
    if abs(Y-X)