Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement