Controlling the Precision of Results
In doing numerical operations like NDSolve
by default uses machine numbers. But by setting the option WorkingPrecision->n
you can tell it to use arbitrary-precision numbers with n
This does a machine-precision computation of a numerical integral.
This does the computation with 30-digit arbitrary-precision numbers.
When you give a setting for WorkingPrecision
, this typically defines an upper limit on the precision of the results from a computation. But within this constraint you can tell Mathematica
how much precision and accuracy you want it to try to get. You should realize that for many kinds of numerical operations, increasing precision and accuracy goals by only a few digits can greatly increase the computation time required. Nevertheless, there are many cases where it is important to ensure that high precision and accuracy are obtained.
Options for controlling precision and accuracy.
This gives a result to 25-digit precision.
Given a particular setting for WorkingPrecision
, each of the functions for numerical operations in Mathematica
uses certain default settings for PrecisionGoal
. Typical is the case of NDSolve
, in which these default settings are equal to half the settings given for WorkingPrecision
The precision and accuracy goals normally apply both to the final results returned, and to various norms or error estimates for them. Functions for numerical operations in Mathematica
typically try to refine their results until either the specified precision goal or accuracy goal is reached. If the setting for either of these goals is Infinity
, then only the other goal is considered.
In doing ordinary numerical evaluation with N[expr, n]
automatically adjusts its internal computations to achieve n
-digit precision in the result. But in doing numerical operations on functions, it is in practice usually necessary to specify WorkingPrecision