function [x] = Gs_100()

num=100;


A=zeros(num,num);
b=zeros(num,1);


b(1)=2;
b(num)=2;

for i=2:num-1
    b(i)=1;
end

for i=1:num
    A(i,i)=3;

end

for i=2:num
    A(i-1,i)=-1;
end


for i=1:num-1
    A(i+1,i)=-1;
end


A

b

d=zeros(num,1);
ep=zeros(num,1);


d(1)=A(1,1);
ep(1)=b(1)/d(1);


for i=2:num-1
    
    d(i)=A(i,i)-A(i-1,i)*ep(i-1);
    ep(i)=b(i)/d(i);
    
end

d(num)=d(num)-A(num-1,num)*ep(num-1);


L=zeros(num,num);
U=zeros(num,num);


for i=1:num
 L(i,i)=d(i);
end

for i=1:num-1
    L(i+1,i)=A(i+1,i);
end

for i=1:num
    U(i,i)=1;
end

for i=1:num-1
    L(i+1,i)=A(i+1,i);
end


for i=2:num
    U(i-1,i)=A(i-1,i);
   
end

fprintf('Aνάλυση LU του πίνακα Α.O U θα είναι άνω τριγωνικός με μονάδες στην κύρια διαγώνιο και ο L κάτω τριγωνικός\n')

U
L

y=GaussSeidelFunction(L,b);

x=GaussSeidelFunction(U,y);

end