g=[5 -6 4 3; 3 -3 2 2; 4 -5 2 1]; wektor=g; n=3; for i=1:(n-1) for j=(i+1):n mn=g(j,i)/g(i,i); for k=i:(n+1) g(j,k)=g(j,k)-(g(i,k)*mn); end end end a x(1,n)=g(n,n+1)/g(n,n); for i=(n-1):-1:1 suma=0; for j=(i+1):n suma=suma+(g(i,j)*x(1,j)); end x(i)=(1/g(i,i))*(g(i,n+1)-suma); end x %%Wektor rozwiązań wektor(:,1:3/wektor(:,4)) %% Macierz po eleminacji