Questions on the numerical functions of Mathematica, implementing numerical methods and numerical computing with Mathematica.
-1
votes
1answer
44 views
How to extract data of large function
Here I have one problem with command for I guess, but I am not sure. If I have function f[x] but not enough ram memory to save it in closed form, how to make code to compute it for chosen values. It ...
1
vote
3answers
84 views
Making a calculation with high precision
I would like to make the following calculation:
1/Sqrt[1 - (150^2 10^(-4))/(9 10^16.)] - 1
Mathematica 8 returns 0. The result is obviously not 0, but my ...
-2
votes
0answers
52 views
2
votes
0answers
55 views
Speeding up a numerical constrained quadratic optimization
I'm trying to solve a quadratic optimization problem in 35 variables, $\vec{α} = \left< α_1, \ldots, α_{35}\right>$:
$$
\begin{aligned}
&\operatorname*{maximize}_\vec{α}&&1.0\cdot ...
10
votes
2answers
350 views
Is it possible to use the LevenbergMarquardt algorithm for fitting a black-box residual function?
I have a black-box multiargument multiparametric function of the type SRD[dataPoint_List,params_List] which accepts experimental data along with the parameters of ...
2
votes
2answers
87 views
28
votes
1answer
487 views
0
votes
1answer
134 views
DAE - varying initial conditions
I want to solve a DAE-system and I want to vary more than one initial conditions and to manipulate them. I looked here:
Putting NDSolve into ParametricPlot
But it does not work:
...
2
votes
0answers
31 views
NIntegrate/NSum with parameters [duplicate]
I'm trying to calculate a continuous integral within a discrete integral.
Something similar to this (yet more complex):
...
4
votes
2answers
105 views
10
votes
1answer
219 views
Why can't I change the value of MaxRecursion in NIntegrate when integrating BesselJ?
I am trying to evaluate this integral numerically
$$
\int_0^{\infty } J_0(q R) \tanh(q) \, \mathrm{d}q
$$
for large values of $R$. This makes the integrand oscillate more quickly and Mathematica ...
0
votes
1answer
68 views
0
votes
0answers
62 views
7
votes
0answers
140 views
Numerically solve 2nd order differential equation with singularity
Consider a second order differential equation with a potential that diverges at some generic value in the variable. For example:
$$-y^{\prime\prime}(s)+\frac1{\mathrm{cn}{(s\mid k^2)}}y(s)=0$$
where ...
1
vote
1answer
94 views
why there is a small imaginary part [closed]
I encountered a problem. I have a eigenvector eigvsI[1]
...
1
vote
1answer
167 views
Why is arithmetic faster for inexact arithmetic?
I have been trying to compute eigenvalues of a rather sizable matrix A, about $500 \times 500$ (but sparse). I asked Mathematica to compute ...
4
votes
1answer
166 views
NDSolve for a large system of simple ODEs
I am solving a system of many (more than 100) ODEs.
It is the kind of standard rate equation encountered in semiconductor physics.
Here is the system:
...
4
votes
2answers
293 views
Any ideas on how GeneralMiniMaxApproximation is implemented?
GeneralMiniMaxApproximation is used to construct minimax approximations of parametrically defined functions. I am curious about how ...
4
votes
1answer
73 views
Minimize failing on a polynomial
Calling:
Minimize[{-0.4877 - 0.1190 r^2 - 0.1885 r^4 + 2.9703 z - 0.5531 z^2,
0 <= z <= 3.5 ∧ 0 <= r <= 1.75}, {r, z}]
returns ...
1
vote
0answers
114 views
Adapting NDSolve to circumvent NDSolve::bdord: error for 1-D Euler Equations
I attempted to use NDSolve for the 1-D isentropic unsteady flow equations with low subsonic inflow velocity and prescribed inflow total enthalpy; along with a ...
6
votes
2answers
213 views
Numerical partial derivative
For a one-variable numerical function it's simple to calculate the derivative at a point with Derivative as @Szabolcs has pointed out before:
...
3
votes
0answers
119 views
FindRoot gives a wrong solution which obviously should not be there
I got stuck on FindRoot and I didn't see any similar problem posted, so let me explain what I am trying to do and what problem I meet here.
I try to find roots of a particular function, which in the ...
5
votes
1answer
122 views
FindMaxValue specifics
I'm using FindMaxValue to study the distribution of maxima of Abs[RiemannSiegelZ[t]] between consecutive values of ...
0
votes
0answers
41 views
How to force evaluation/numerical result of a function? [closed]
I defined a function m[x] using
...
1
vote
0answers
73 views
Why is FindRoot initial value far from the specified one?
I am trying to numerically find the root of a function that looks a bit like: 1/x - (SchurDecomposition[A[x]][[2]])[[1]], where ...
7
votes
2answers
382 views
How do you force a decimal output? [duplicate]
I have some very small values such as 2.601519253*10^-8. I'd like to output these values to CSV for another program to work with. I've tried N[value, 50], but Mathematica still insists on producing ...
0
votes
0answers
153 views
Rounding values [duplicate]
In Mathematica 8 when I enter 1 - 0.99 - 0.01 I get 8.67362*10^-18 instead of zero. How do I fix this problem?
I am getting ...
13
votes
1answer
271 views
Optimizing a Numerical Laplace Equation Solver
Laplace's Equation is an equation on a scalar in which, given the value of the scalar on the boundaries (the boundary conditions), one can determine the value of the scalar at any point in the region ...
10
votes
6answers
2k views
About multi-root search in Mathematica for transcendental equations
I have some questions for multiroot search for transcendental equations. Is there any clever solution to find all the roots for a transcendental equation in a specific range?
Perhaps ...
-1
votes
0answers
56 views
Why is NonlinearModelFit calling the model function with symbolic arguments?
I have a model to fit which is made of numerical functions and so it can give a value for any numerical value of the arguments, but creates an error when called with symbols. ...
6
votes
1answer
461 views
Handling failed FindRoot calls
I want to handle FindRoot calls which did not converge (e.g "thrown" error message FindRoot::cvmit)
...
1
vote
1answer
66 views
Plot FindRoot for non-trivial function
I would like to plot the results of FindRoot over certain range of inputs. I tried to do this with the code:
...
3
votes
1answer
75 views
Find point at which equation stops having roots (if it exists)
I am interested in the roots of this function:
f[M_, b_] := 1 - (2 M Gamma[2, 0, (1/M + b M)/Sqrt[b]])/(1/M + b M)
for fixed values of b. In particular I want ...
1
vote
2answers
127 views
How can I use FindRoot on an expression from NDSolve?
I have a second order ODE that I can only solve numerically using NDSolve, but I then need to use the solution in FindRoot and am running into errors. A simplified but analogous problem is the ...
34
votes
10answers
976 views
Can Mathematica propose an exact value based on an approximate one?
Sometimes, I use Mathematica to do some hypothesis on homeworks to make the question easier. For instance, when I have to compute big sums when $n\to\infty$ and Mathematica can't give the exact ...
19
votes
2answers
1k views
Why round to even integers?
According to the Mathematica help:
Round rounds numbers of the form x.5 toward the nearest even integer.
For example:
Round[{0.5, 1.5, 2.5, 3.5, 4.5}]
...
0
votes
1answer
80 views
How to guess initial complex value for FindRoot
I have to solve a transcendental equation for a parameter, say $\beta$. Now, the $\beta$ has a range from $ik$ to $k$ where $i$ is the usual imaginary root $\sqrt{-1}$ and $k$ is a real number. ...
0
votes
0answers
62 views
FindMaximum inconsistency
The code below seems to work for n<11. But for n=11, and above, it outputs newa then just outputs "beep" sound.
WhyTheBeep says "The kernel Local has quit ...
3
votes
1answer
179 views
Tutorial for basic numerical methods for PDEs
I'm afraid this is probably not going to be a "good" question, but I'd like to use Mathematica to learn about basic numerical schemes for solving pdes. For example, I'd like to compute the solution of ...
14
votes
5answers
768 views
The difference between 0. and 0
I have a function for which 0 is a special case:
f[A___, 0, B___] := 0
But since I am doing numerics, sometimes in the course ...
1
vote
0answers
62 views
Using Root for numerical evaluation [closed]
In my expressions I have a lot of roots to calculate. For example:
Root[#1^3 + #1^2 + #1 &, 1]
Since most of them cannot be solved analytically, I would like ...
2
votes
2answers
139 views
Position function not always retuning an answer even with no apparent problems
I'm having some problems with Position.
Sometimes it will give an empty list instead of the actual position of the element I am looking for when that element is ...
3
votes
1answer
99 views
Strange behaviour of PolyLog Function
I discovered some strange behaviour of the PolyLog[] Function in Mathematica which seems to me like a bug in the function implementation.
It looks like ...
0
votes
0answers
304 views
Jacobi eigenvalue and eigenvectors algorithm
Suppose we have a symmetric matrix with dimensions n x n. I need to find the eigenvalues and eigenvectors of this matrix using the Jacobi method. I wrote this code:
...
0
votes
0answers
84 views
Speeding up numerical computations
So, after seeing this and this post, I was thinking whether there were any general rules for speeding up numerical computations (not only NIntegrate or ...
2
votes
1answer
71 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 ...
-1
votes
1answer
130 views
Minimization problem using FindMinimum. Functional value “is not a real number”? [closed]
This question is a follow-up of another one I asked a few days ago. I followed the instructions given in the answer provided by the user that responded. I modified that answer to solve another ...
3
votes
2answers
135 views
7
votes
0answers
281 views
Dual complex integral over implicit path using contour plot
Context
I am interested in doing double contour integral over paths which are defined implicitely.
For the sake of debugging, let's assume its
$$\oint_{\cal C}\oint_{\cal C} \frac{1}{u\, x} d u d x$$ ...
1
vote
1answer
98 views
What do these error messages mean when numerically solving differential equations?
I am using Mathematica to numercially solve the following equations:
...
