From: Mohammad on
Here is a simple program that does Gaussian elimination:

for j=2:r,
for i=j:r,
mat(i,:) = mat(i,:) - mat(j-1,:)*mat(i,j-1)/mat(j-1,j-1);
end
end

I works, but not to the extant I want it. It considers values:
0.0204327942426589
-0.000401333398864821
0.00660003568422868
0.00856120376245412

to be equal to zero, such that when I do substitution, my answers are a little of when I use a few variables (like ten). But when I substitute for 230 variables, I get answers way off. Is there a way to make this program much more accurate?

Thanks!

Cordially, Mohammad