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