Metoda Jacobiego

(function x=jacob(A,b,xO,e) %funkcja pobiera macierz A, wektor b, początkowy wektor przybliżeń i warunek zbieżności e n=size(A); %pobranie rozmiaru macierzy

ek=l; %przypisanie początkowych wartości dla warunku zbieżności ek, wektora x i licznika iteracji x=x0 ; it=0;

while ek>e %pętla sprawdzająca warunek zbieżności

xold=x; %przypisanie wartości wyliczonej w poprzedniej iteracji pod zmienną xold

Tfor i=l:n %pętla po kolumnach

s=0; %początkowa wartość sumy

for j=- :n %pętla sumująca iloczyny wyrazów macierzy A i wektora x z pominięciem aktualnie liczonego wyrażenia if j!=i %sprawdzenie pcv;yższegc warunku s=s+A(i,j)*xold(j); %suma

end

end

x(i)=(b(i)-s)/A(i,i); %obliczanie i-tego wyrazu v;ektora x it=it+l; %aktualizacja licznika iteracji

end

ek=3qrt(dot((x-xold),(x-xold))); %obliczenie liczby warunkującej zbieżność

if it>500 %sprav;dzenie czy nie zestala przekroczona liczba dozv7olonych iteracji disp('proces rozbieżny’); break; %przerv;anie działania pętli

end

end