% Script sgausstest.m

% This example has no zero pivots

A = [ [6,-2,2,4]; [12,-8,6,10]; [3,-13,9,3]; [-6,4,1,-18] ];
b = [ 16; 26; -19; -34];

x = sgauss(A,b) % doesn't change A, b
A*x - b

% This example has a zero pivot

A = [ [1,1,3,1]; [2,2,-1,2]; [2,3,42,0]; [2,3,4,4] ];
b = [-1; 5; 3; 3];

x = sgauss(A,b)
A*x - b