Eliminacja Gaussa

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

DHD

for i= :n %pętla po kolumnach

for j=i+ :n %pętla po wierszach (poniżej aktualnego elementu przekątnej głównej)

m(j,i)=A(j,i)/A(i_fi); %oliczenie mnożnika dla eliminacji w kroku odpowiadającym wyrazowi A(i,i) for k=i+l:n %pętla wykonująca odejmowanie wierszy A(j,k)=A(j,k)-m(j,i)*A(i,k);

end

b(j)=b(j)-m(j,i)*b(i); %obliczanie nowego wektora b(wektor wyrazów wolnych)

end

end

x=3ackS(A,b) %wywoianie funkcji podstawiania v;stecz end