gaus

function[W] = met_gaus(X,b)
[n,~]=size(X);
W=zeros(n,1);
for i=1:n-1
    m=X(i+1:n,i)/X(i,i);
    X(i+1:n,:)=X(i+1:n,:)-m*(X(i,:));
    b(i+1:n,:)=b(i+1:n,:)-m*(b(i,:));
end
     W(n,:)=b(n,:)/X(n,n);
   for i = n-1:-1:1;
       W(i,:)=(b(i,:)-X(i,i+1:n)*W(i+1:n,:))/X(i,i);
       end
end