gauss1

% http:/Vwww.algorytm.org/index.php?option=com_contentóitask=viewd(id=82<&rtemid=28 function X=gauss(A,b) % A - macierz rownan dla X-niewiadome, równanie: Ax=b

if(size(A,l) = =Size(A,2) & Size(A_/l) = = Size(b,l)) %spr. poprawnego rozmiaru macierzy A=[A b]; % łączenie macierzy

for(k=l:size(A,l)) % ile wierszy, tyle zerowań, dopóki to nie osatni wyraz w wierszu współczynników x

1=1;

for(j=(k+l):size(A,l)) % indeks kolumny for(i=l :size(A,2)) % szukanie elementu różnego od zera if(A(k,i)~ = 0 & A(k,k)) % jeśli różny od zera - wykonuj ciąg dalszy A(j,:)=A(j_/:)-A(k,:)*A(j_/k)/A(k,k); % odejmnij przemnożonożony współczynnik %A

break

elseif(i==size(A,2)) disp('Nieskończenie wiele rozwiązań.'); k=0; X=k; return end end end

if( i<Size(A,2) & k==Size(A,l)) % etap drugi - "postępowanie odwrotne"

% dispCpart 2');

X=A(: ,Size(A,2)); % wyrazy wolne ¹ b¹ X(k)=X(k)/A(k,k);

1=0;

whilef k=k-l)_