errors

errors, a Fortran90 code which demonstrates how reasonable computations can produce numerical nonsense. This illustrates that the programmer must not assume that a numerical algorithm that seems reasonable will always produce correct and reliable results.

The computations include polynomial evaluation and root finding, linear system solution, minimization, and Taylor series approximation.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

errors is available in a Fortran90 version and a MATLAB version and an Octave version.

Related Programs and Data:

errors_test

Reference:

1. U Kulisch, C Ullrich, Editors,
   Wissenschaftliches Rechnen und Programmiersprachen,
   (Scientific Computing and Programming Languages),
   Berichte des German Chapter of the ACM,
   (Reports of the German Chapter of the ACM),
   Volume 10, Teubner Verlag, 1982.
2. Cleve Moler, Charles Van Loan,
   19 Dubious Ways to Compute the Exponential of a Matrix, 25 Years Later, SIAM Review,
   Volume 45, Number 1, pages 3-49, March 2003.
3. Yves Nievergelt,
   Numerical Linear Algebra on the HP-28, or How to Lie with Supercalculators,
   The American Mathematical Monthly,
   Volume 98, Number 6, June-July 1991, pages 539-544.
4. Siegfried Rump,
   Wie Zuverlaessig Sind die Ergebnisse Unserer Rechenanlagen?
   (How Reliable are the Results of our Computations?)
   Jahrbuch Ueberblicke Mathematik 1983, pages 163-168.

Source Code:

• errors.f90, the source code;
• errors.sh, compiles the source code;