Proszę o sprawdzenie, ewentualnie wytłumaczenie błędów. I mam do tej metody napisać wykres zbieżności, tylko nie wiem jak, podejrzewam, że będzie to coś w stylu
Epsilon<=norm(xn-xp) ??
Moja "twórczość":)
Kod: Zaznacz cały
A1=[8, 1, 2, 3; 1, 6, 1, -1 ;2 , 1, 18, 2; 3, -1 ,2 ,40]
B=[28, 12, 66, 167] A=[8, 1, 2, 3, 28;1, 6, 1, -1 , 12;2, 1, 18, 2, 66;3, -1, 2, 40, 167]
x=[1;1;1;1]
xp=x
Vp=xp/norm(xp)
for k=1:30
x(1)=(A(1,5)-A(1,4)*x(4)-A(1,3)*x(3)-A(1,2)*x(2))/A(1,1)
x(2)=(A(2,5)-A(2,4)*x(4)-A(2,3)*x(3)-A(2,1)*x(1))/A(2,2)
x(3)=(A(3,5)-A(3,4)*x(4)-A(3,2)*x(2)-A(3,1)*x(1))/A(3,3)
x(4)=(A(4,5)-A(4,1)*x(1)-A(4,2)*x(2)-A(4,3)*x(3))/A(4,4)
xn=x
xp=xn
Vp=xp/norm(xp)
end