Tagged Questions

The tag has no usage guidance.

69 views

Error In polynomial Root Finding

I have a polynomial in y like so: 2.00855*10^20 + 6.89796*10^20 x y + (5.17347*10^20 + 5.92241*10^20 x^2) y^2 - 1.4806*10^21 x y^3 + 7.77316*10^20 y^4 == 0 I ...
176 views

Is manual adjustment of AccuracyGoal and PrecisionGoal useless?

This is a problem confusing me for years. AccuracyGoal and PrecisionGoal are two options that I never truly understand and, to ...
29 views

Use of AccuracyGoal & MachinePrecision in NIntegrate

I just learned about the Inverse Symbolic Calculator and it seems like a very useful tool. For example, to find an analytic solution to  \int_0^{\infty } \frac{1-e^{-x} (x+1)}{\left(e^x-1\right) \...
156 views

Observe that ...
744 views

Why is Mathematica destroying this graph?

Here I have a picture of a function I graphed: reg[x_,y_]:=(x^2+y^2)Cos[4ArcTan[y/x]]; Plot3D[reg[x,y],{x,-2,2},{y,-2,2},AxesLabel->Automatic] And here is ...
26 views

Implementing more accurate boundary conditions in NDSolve

I want to solve a second order differential equation numerically. The boundary conditions are needed to be imposed at z=0 and at ...
193 views

What's the proper way to 'cut out' an irrelevant part of an integral, so I don't run into problems with precision?

I'm doing a numerical integration that involves an integrand with a few exponentials whose values over the integration region range enormously. I'm really only interested in the part of the range that ...
50 views

My NIntegrate expression returns a wildly inaccurate value

I am trying to integrate a function using NIntegrate: ...
98 views

Funny behavior when computing dot product of coefficients with high-order polynomials

I have a similar problem to Funny behaviour when plotting a polynomial of high degree and large coefficients. However, the thing being evaluated is not just a polynomial but a dot product of some ...
77 views

WorkingPrecision and AccuracyGoal for using with NIntegrate

I'm modeling a quantum tunneling barrier. Because I want to express my quantities in mks units, that means that some of the variables I'm using are pretty small (10^-9). I use NIntegrate for part of ...
131 views

Why has Version 10.3 precision reduced?

In version 7.0.1.0 and versions 10.0 and 10.1 the following is produced: ...
79 views

Improve Accuracy of FindRoot

This is a follow up question to an earlier questions I asked. The link is: Product of $N\ 2 \times 2$ matrices and subsequently solving an equation dependent on the product As I am new here, I am not ...
64 views

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

Why is it so, N[Eigenvectors[{{1., 0.5, 4.}, {2., 1., 2.}, {0.25, 0.5, 1.}}]] returns: ...
41 views

Maintaining working precision in my program

I am badly stuck with this program I am writing. I want output of program below to be precise up to 25 digits after decimal point. I am supplying all inputs with precision 25. But after each iteration ...
105 views

Unexpected behavior from Accuracy

This is my code Table[With[{x = 10^n + 1/17}, N[x, {Infinity, 5}]], {n, 0, 5}] // Column Or like this ...
81 views

Exact calculating the multiply of two real [closed]

As the machine calculate,the exact calculating with real number have some problem to puzzle me long time.We know the 0.142*0.36523=0.05186266 is exactly equalized ...
581 views

more numerically accurate inverse matrix

I encountered the following matrix mat = {{2, 2.161209223472559 + 1.682941969615793 I}, {2.161209223472559 - 1.682941969615793 I, 2}} and ...
73 views

Handling Accuracy and Simplifying

I have the following problem. I have a system of non-linear equations that I log-linearize around a certain point, let's call it point A, using a function that I ...
81 views

$\tt PrecisionGoal$ and $\tt AccuracyGoal$ making the solution less accurate?

I am having this differential equation: ...
66 views

Accurate Integration of 2D interpolating function with weighting

I want to integrate a 2D interpolating function received from NDSolve[]. To illustrate my problem I will give a simple example: ...
270 views

NIntegrate fails to converge under almost any PrecisionGoal, MinRecursion etc. How can I trust the result?

I have been getting some ideas by reading other related questions in the forum, but the integral I have to do is not converging in many cases. The integrand is of the form: ...
60 views

Problem with precision of fraction numbers

I have tried to take a series of harmonic numbers using Mathematica and its precision but there has been an issue. So far when I computed the sums at a whole numbers using a precision of 100 digits I ...
71 views

Problem with calculating Harmonic Numbers

I have tried to take a series of Harmonic Numbers using mathematica but there have been issues in calculations. So far when I computed the sums at a small value range such as ...
76 views

In NMinimize, how can I efficiently handle NIntegrate errors for non-integrable functions?

I am using NMinimize to find parameters that minimise the integral of a function in the least squares sense: ...
617 views

Setting the Accuracy of calculations

I need to optimize an expression that involves a number of trigonometric functions and Exp[]. How do I make sure that all my calculations have an accuracy of 120-...
80 views

Is there a way for accuracy backtick notation to accept variable inputs?

Backtick notation uses 1.020 to give 1 with 20 places of accuracy after the decimal point. But what I want is ...
333 views

3D Gamma function

I try to plot the following function: Plot3D[Gamma[1+0.5*(n+m)]/Sqrt[Gamma[1+n]*Gamma[1+m]],{n,0,1000},{m,0,1000}] I expect that for m=n and near to it the value ...
95 views

NDSolve returns asymmetric solution to a symmetric equation

I am trying to find the potential of a conducting cylindrical electrode by solving the Laplace's equation. Both, the boundary conditions and the equation are symmetric w.r.t. the change $r\to-r$. ...
163 views

Accuracy of Grad-Shafranov PDE's NDSolveValue over implicit region

Context I need to solve for the toroidal flux of the magnetic field above an accretion disc. For this purpose, I define the region over which the flux is non zero, ...
269 views

How to add AccuracyGoal and PrecisionGoal options in a numerical function?

There are many built-in functions that contain AccuracyGoal, PrecisionGoal and ...
171 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 ...
433 views

Precision and accuracy in NDSolve and NMinimize

I use NDSolve quite a lot and have noticed that setting values for AccuracyGoal and ...
54 views

If statement for evaluating small values

I have a simple calculation which I am trying to use if condition to evaluate summation of some values, here's a code; ...
7k views

Understanding differences between Maple and Mathematica in examples picked by Maplesoft

I am reading the document How Maple Compares to Mathematica. On page 15 there is an example where Mathematica produces wrong results. Does anybody know why? MAPLE: MATHEMATICA: Also on page 17 ...
220 views

Precision while calculating Fourier transform

I'm trying to understand how precision works in Mathematica. Particularly I'm calculating discrete Fourier transform using the Fourier function and calculating it "...
227 views

Precision of NIntegrate

At the moment I am considering a "difficult", highly-oscillatory integral in Mathematica. It calculates the integral without any complaints. However, I am also trying out a numerical method with which ...
77 views

Negative accuracy numerics ( 0-128 notation )

What does it mean to have a negative accuracy number? I understand 0128 to mean "the number zero to 128 decimal points". Mathematica corroborates this: ...
439 views

How trustworthy is NMaximize?

Suppose I solve a constrained optimisation problem using NMaximize. How confident can I be of the accuracy of the result? For concreteness, suppose that F,G are (...
284 views

How to get the digits after the decimal point

I have N[Sqrt[2], 8] which outputs 0.4142136 How do I get the 4142136 part out of it? I ...
221 views

Numerical Integration different in Mathematica version 9 and 10 with same options

I have noted that the same function with the same settings gives me different results in Mathematica version 9 and 10. This involves integrating numerically interpolating functions and so on. Here a ...
924 views

Floating point operations — division vs multiplication; how do they affect accuracy?

Given two functions: f[x_] := x ((x + 1)^(1/2) - x^(1/2)) g[x_] := x/((x + 1)^(1/2) + x^(1/2)) Which one is more accurate? Side note: If you could explain why, ...
108 views

Option: number of digits

The command NestList[Cos, 1.5, 10] produces the list: ...
704 views

How to improve the solution of a double pendulum system

First, the code: ...
88 views

User defined functions and decimal places displayed

I have a couple of complicated user-defined functions, one which calls the other, and I wish to control the number of decimal places displayed when I call the "mother" function. I have tried the ...
154 views

the default option of NDSolve

How should I find out what the default value of Accuracy and Precision is in the NDSolve? I have a pde that when I asked it to solve for me the equation with a special Accuracy and precision it gave ...
92 views

How to solve the warning problem and obtain real roots without imaginary part?

I am trying to solve a equation with Newton's method via FindRoot, and the codes are: Define the functions: ...
59 views

How to check the decimal accuracy of result? [duplicate]

I am trying to check if the result is enough accurate to use. I need an accuracy of 10^-6 from some numerical method. Something like root finding algorithm. But how do I check that in mathematica? I ...