function X $=$ myGaussJordan$($A$,$ B$)$

n $=$ size$($A$,1)$;

X $=[$A$,$ B$]$;

for i $=1$:n

$[$~$,$ max$\_$row$]=$ max$($abs$($X$($i:n$,$ i$)))$;

max$\_$row $=$ max$\_$row $+$ i $-1$;

X$([$i$,$ max$\_$row$],$ :) $=$ X$([$max$\_$row, i$],$ :$)$;

pivot $=$ X$($i$,$ i$)$;

if abs$($pivot$)<1$e$-10$

error$(\text{'}$Matrix is singular or poorly conditioned.$\text{'})$;

end

X$($i$,$ :) $=$ X$($i$,$ :) $/$ pivot;

for j $=1$:n

if i ~$=$ j

X$($j$,$ :) $=$ X$($j$,$ :) $-$ X$($j$,$ i$)*$ X$($i$,$ :$)$;

end

end

end

X $=$ X$($:$,$ n$+1$:end$)$;

end

A $=[2,1,-1$;

$1,3,2$;

$1,0,3]$;

B $=[4$;

$3$;

$5]$;

X $=$ myGaussJordan$($A$,$ B$)$;

disp$(\text{'}$Solution:$\text{'})$;

disp$($X$)$;