Gaussian elimination: Difference between revisions

m
→‎version 3: added whitespace to the REXX section header wording.
(Added Wren)
m (→‎version 3: added whitespace to the REXX section header wording.)
Line 5,102:
This is the same as version 2, but in addition, it also shows the residuals.
 
Code was added to this program version to keep a copy of the original   '''A.i.k'''   and   '''B.#'''   arrays   (for calculating the residuals).
<br>residuals).
<br>Also added was rounding the residual numbers to zero if the number of significant decimal digits was
 
less or equal to 5% of the number of significant fractional decimal digits &nbsp; (in this case, 5% of 1,000 digits for the decimal fraction).
<br>Also added was the rounding the residual numbers to zero &nbsp;if&nbsp; the number of significant decimal digits was &nbsp; <big> &le; </big> &nbsp;5%&nbsp; of
less or equal to 5% of <br>the number of significant fractional decimal digits &nbsp; (in this case, &nbsp;5%&nbsp; of &nbsp;1,000&nbsp; digits for the decimal fraction).
<lang rexx>/*REXX program solves Ax=b with Gaussian elimination and backwards substitution. */
numeric digits 1000 /*heavy─duty decimal digits precision. */