Gaussian elimination: Difference between revisions
Content deleted Content added
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
<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 (in this case, 5% of 1,000 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. */
|