Systems of linear equations Our Matlab function for naive Gaussian elimination looks like this: function x = naiv_gauss(A,b) n = length(b) x = zeros(n,1) for k=1:n-1 forward elimination for i=k+1:n xmult = A(i,k)/A(k,k) for j=k+1:n