Metoda Choleskiego Banachiewicza

function x=ChB(A,b) %funkcja pobiera macierz A i wektor wyrazów wolnych b, a zwraca wektor rozwiązań x n=size(A); %pobranie rozmiaru macierzy

3for	k= : n s=);
¥	for p-

end

FJfor k=^ :n %pętla po kolumnach macierzy L i wierszach macierzy L*

3=s+L(k,p)^A2;

L(k,k)=sqrt(A(k,k)-s); %obliczanie k-tegc wyrazu głównego z użyciem wyżej obliczonej sumy (wg. wzoru z wykładu) for i=k+l :n

s=0; Iwyzerowanie wartości zmiennej s

for p= :k-- %pętla sumująca iloczyny wyrazów macierzy L s=s+L(i,p)*L(k,p);

end

L (i, k) = (A (i, k)-3)/L (k, k) ; %cbliczenie wyrazu (nie z przekątnej głównej) z użyciem policzonej sumy (wg. wzoru z wykładu)

end

end

y=ForwardS(L,b); %wywołanie funkcji podstawiania w przód w celu obliczenia równania L*Y=3 => Y (bc L jest dolnotrójkątna) x=3ackS(L',y); %wywołanie funkcji podstawiania wstecz w celu obliczenia równania U*X=Y => X (bo U jest górnotrójkątna) end