// Asier Aranzabal
// OCW-EHU
// http://ocw.ehu.es/irakaskuntza-teknikoak/zenbakizko-kalkulua-ingeniari-kimikoentzat/Course_listing
// Zenbakizko kalulua

function x=gauss(A,b)
//
// Ax=b ekuazio sistema lineala ebazteko Gauss-en ezabapen metodoa erabiltzen du.
//
// Sarrerako baliak:
//   A= Aldagaien koefizienteak dituen matrizea.
//   b= Emaitzak biltzen dituen bektorea.
//
// Irteerako balioak
//   A= Goiko matrize triangeluarra
//   b= emaitza bektorea
//   x= ebatzitako bektorea.
//
//
//
//ezabapena
l=size(A);n=l(1);
for k=1:n-1;
   for i=k+1:n;
      z=A(i,k)/A(k,k);
      for j=k:n;
         A(i,j)=A(i,j)-z*A(k,j);
      end
      b(i)=b(i)-z*b(k);
   end
end
x(n)=b(n)/A(n,n);
//
// atzerako ordezkapena
for i=n-1:-1:1;
   sum=0;
   for j=n:-1:i+1;
      sum=sum+A(i,j)*x(j);
   end
   x(i)=(b(i)-sum)/A(i,i);
end
x=x'; // ilara bateko bektorea zutabe bihurtzen du

endfunction