% Script zgausstest.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 = zgauss(A,b)
A*x - b

% This example has a zero pivot

A = [ [1,-1,2,-1]; [2,-2,3,-3]; [1,1,1,0]; [1,-1,4,3] ];
b = [-8; -20; -2; 4];

x = ngauss(A,b)
x = zgauss(A,b)
A*x - b