The tag has no wiki summary.

learn more… | top users | synonyms

44
votes
3answers
1k views

When I can assume that all decimal digits returned by Mathematica are provably correct?

How to Control the Precision and Accuracy of Numerical Results Arbitrary-Precision Numbers Mathematica works with exact numbers and with two different types of approximate numbers: ...
12
votes
2answers
295 views

Why do NumberForm and Round apparently use different tie-breaking methods?

When rounding numbers (for example, rounding a real number to the nearest integer), the "round to nearest" rule is usually used. For example, 1.4 is rounded down to 1 and 1.6 is rounded up to 2. ...
8
votes
2answers
180 views

Error computing sum of sum of digits

I've defined a function that computes the sum of the base-b digits of n: DigitSum[n_, b_] := Total[IntegerDigits[n, b]] Then I defined a function that computes ...
7
votes
1answer
269 views

Confused by (apparent) inconsistent precision

$$ e^{\pi \sqrt{163}} \approx 262537412640768743.99999999999925 $$ E^(Pi Sqrt[163.0]) N[E^(Pi Sqrt[163.0]), 35] NumberForm[E^(Pi Sqrt[163.]), 35] returns ...
6
votes
3answers
221 views
6
votes
3answers
1k views

Strategies to solve an oscillatory integrand only known numerically

I have an integrand that looks like this: the details of computation are complicated but I only know the integrand numerically (I use NDSolve to solve second ...
6
votes
2answers
203 views

A problem about fixed point iteration theory

Description Recently, I have been learning a couse called "Numerical Analysis". The fixed point iteration theory was introducted to solve the ...
6
votes
2answers
146 views

Large Number times Tiny Number

A (maybe noob) question: Let a = 1234567891234567889998.5 b = 1234567891234567889999.5 Mathematica (v8) yields for ...
6
votes
2answers
270 views

Cannot Get Numerical Results to Match

I try this numerical summation (in two parts) ...
6
votes
0answers
59 views

Is there a way to get BSplineFunction beyond MachinePrecision?

I have been using BSplineFunction (Mathematica 8) to generate a smooth representation of some data, which I might then do some further processing on. As I ...
5
votes
2answers
1k views

How to convert a string of number into real numbers without precision loss

I searched and did not find similar answers which can resolve the 'no precision loss' issue. A sample string is: ...
5
votes
1answer
465 views

ReadList can't read from file containing high-precision Reals

I want to read a file which contains seven columns of real numbers. Mathematica seems to have a problem with the high precision I need. My code: ...
5
votes
1answer
121 views

Numerical error in Mathieu functions

Consider the MathieuCharacteristicA function, which is a piecewise function according to the documentation. The discontinuity happens at integer number. ...
5
votes
1answer
55 views

Mysterious behavior of Precision for complex arrays

Mysterious behavior of Precision: {{1.0+I*0.0},{0.0+I*0.0}} // SetPrecision[#,30]& // Precision // Print; ...
4
votes
6answers
1k views

Why is this Mandelbrot set's implementation infeasible: takes a massive amount of time to do?

The Mandelbrot set is defined by complex numbers such as $z=z^2+c$ where $z_0=0$ for the initial point and $c\in\mathbb C$. The numbers grow very fast in the iteration. ...
4
votes
2answers
184 views

Why does N[Re@f] give complex result?

Consider this code: ...
4
votes
3answers
184 views

Increasing the precision of a calculation

Clearly I have misunderstood how to do this: I seem unable to understand how to get more digits out of a calculation and I can't see why - the following piece of code illustrates the issue. ...
4
votes
1answer
96 views

Issues with $MachineEpsilon

I'm attempting to add $MachineEpsilon to numbers I am pulling from the domains of 3 interpolating functions. I pulled 3 numbers, and for one of these numbers adding ...
4
votes
1answer
158 views

Padding plot ticks with zeros on the right

I would like to know how can I make a real-valued plot tick pad with a zero to the right of the decimal point on integer values. This is what I have: ...
4
votes
2answers
56 views

numerical precision in output file

I want to write the numerical results of a computation to a file which is to be read in Fortran. The data looks like: ...
4
votes
1answer
194 views

How to find derivative of a numerical solution, where precision is ambiguous?

I am trying to take the derivative of a numerical solution. I am concerned that the way I'm doing this may be problematic due to numerical error; I think there must be a better way but I'm not very ...
4
votes
1answer
269 views

ParallelTable and Precision

I'm using ParallelTable[] to calculate a function over a range of my parameters , ($\omega,\ell$). This seems to be working well (in terms of speed increase) except ...
4
votes
1answer
721 views

How to set the working precision globally? $MinPrecision does not work

I want to increase the precision globally to 50 at least, but $MinPrecision does not work. ...
4
votes
3answers
151 views

Importing High Precision (22 digits) data from a file

I want to import list from a file and I am using ReadList but my data is in format like this 8.2457409790900004e+08 but it i being imported like $8.24574*10^8$. Is ...
3
votes
4answers
169 views

AccountingForm w/o significant figures

Basically I want to build a chart like something I would do in Excel. It should look like this: Here is the formula I'm using and what I've done so far ...
3
votes
1answer
41 views

Is it possible to change the MinPrecision in command “Compile”?

I'm using "Compile" to do some numerical calculation and I need high Precision. I was wondering if it is possible to change the Precision in command "Compile" to some $MinPrecision?
3
votes
1answer
214 views

Numerical Error with Large Matrices

I am writing a Finite Element Analysis program in Mathematica. The code involves handling a large matrix with large entries. I get an error when I try to use Mathematica's "LinearSolve" to solve a ...
3
votes
2answers
104 views

Increase precision for a entire notebook

I am quite a newbie and it is taking me so much time to understand what I need to do to achieve something really simple and I still haven't found it. I run a dynamic model of a metabolic pathway in ...
3
votes
1answer
38 views

Deferent results depending on the input form of the matrix [closed]

Why is it so, N[Eigenvectors[{{1., 0.5, 4.}, {2., 1., 2.}, {0.25, 0.5, 1.}}]] returns: ...
3
votes
0answers
90 views

Precision problem while evaluating RegionMember or Element in V10

I'll try to give you a background on where the problem occured. The problem itself is defined near the end of the post. I'm running a sort of a simulation of a particle in a box. I finally got my ...
3
votes
0answers
90 views

Tools for Stability/Automatic Error Analysis in Mathematica

I have a longer analytic expression in several variables containing special functions and others. Does Mathematica bring tools - or are there any packages - to examine the stability when I evaluate ...
2
votes
4answers
132 views

Precision problem with numerical solution of a differential equation

I want solve $$ 2\sqrt{|\gamma|}x = \int_{1}^{t} dy \sqrt{\frac{1+2|\gamma_2|y}{y^2(1-y)}} $$ where $0<t\leq $1. I'm using a for cycle to evaluate t, calculate the integral and the assign the $x$. ...
2
votes
2answers
132 views

NIntegrate - WorkingPrecision

I am trying to NIntegrate this expression. It is very close to zero, and so gives all possible types of errors. More precisely, ...
2
votes
1answer
108 views

Why the EuclideanDistance of these 2 exact points generates the warning N::meprec?

The two points are: pts = {{1/4 (-1 - Sqrt[5]), Sqrt[5/8 - Sqrt[5]/8]}, {1/4 (-1 - Sqrt[5]), -Sqrt[5/8 - Sqrt[5]/8]}}; EuclideanDistance @@ pts ...
2
votes
2answers
834 views

Problem with setting working precision in NIntegrate

I want to obtain a good numerical approximation (up to 10 decimal place would be ok for me) to an integral: $$ \int^{\infty}_{0} f(r)r^2dr $$ I am using the function $f(r)$, which is related to the ...
2
votes
1answer
51 views

Working precision and pure functions in FindMinimum

I have this function: fun = (# - 5.0)^2.0 & Then, I want to find the minimum: ...
2
votes
1answer
57 views

Range and decimal digits [duplicate]

When I set Plot[Sin[x], {x, -2, -1}, Frame -> True, Axes -> False, FrameTicks -> {{Automatic, Automatic}, {Range[-2, -1, 0.2], Automatic}}] I ...
2
votes
0answers
119 views

how to calculate with just double floating point precision [duplicate]

I want to calculate with simple double floating point precision. ...
1
vote
2answers
80 views

SetPrecision — different values subtracted from unity gives same value

How are my results possible? CODE: SetPrecision[FibersFailedNew, 40] 3.612171938234261799295174910010284760488*10^-16 ...
1
vote
1answer
106 views
1
vote
1answer
147 views

ListPlot with high precision

I have a list of numbers, e.g. {{1.,0``27.698970004336022},{2.,0``27.698970004336015}} but usually much longer, suitable for ...
1
vote
2answers
67 views

ListPlot3D with small (but acceptable) precision

Consider the following minimal example: ...
1
vote
1answer
111 views

NDSolve fails for certain choices of parameters and solve range

I'm trying to solve a pair of coupled ODEs with NDSolve. I know roughly what the solution should look like (both should give periodic functions, pi/2 out of phase, the amplitude of which damp towards ...
1
vote
1answer
172 views

Numerical Instability?

Current Status Please skip ahead to Update 2: the key questions are now: 1/ Why does Mma generally fail to find the minimum for a well behaved function when the function does not have infinite ...
1
vote
1answer
181 views

Output precision

I've solved some equations using FindRoot and then computed some values. Now when I print the output, I only get a certain precision {{0.01, 496.983, 61.3147, ...
1
vote
1answer
53 views

Increase Precision in Numerical Integration

I need to generate a table of Chebyshev expansion coefficients of trigonometric functions (in this case Cos[2 Pi t] to very high accuracy. Code is: ...
1
vote
1answer
80 views

How to suppress a warning message

If I change $MinPrecision before running ListPlot[{{1, 2}}], there will be a warning in Mathematica 10. But there is no warning ...
1
vote
1answer
64 views

Inaccuracy of Piecewise function? [duplicate]

I have the following issue and can't find a solution so far. ...
1
vote
1answer
56 views

ReadList error with negative very small numbers

When I try to import data with ReadList from a file containing very small numbers (~10^-318), positive data is handled well, ...