GaussMatlab

%% gauss con funcion
function y=FunGauss(A,B)
A=[2 -3 -1;3 2 -5;2 4 -1]
B=[3;-9;-5]
[n n]=size(A);
A=[A';B']';
x=zeros(n,1);
for p=1:n
    for k=[1:p-1,p+1:n];
        m=-A(k,p)/A(p,p);
        A(k,:)=A(k,:)+m*A(p,:)
    end
end
x=A(:,n+1)./diag(A)

% Grafica
[x,y]=meshgrid(-5:05);
z=2*x -3*y -1
surf(x,y,z);
xlabel('X');
ylabel('Y');
zlabel('Z= C');
colorbar
hold on