Questions on the use of numerical functions NIntegrate and NDSolve.

learn more… | top users | synonyms

5
votes
1answer
106 views

Finding minimum fly-by radius between Mars and spacecraft from interpolating function

I've written an interplanetary trajectory solver/plotter that plots the path taken by a spacecraft on an Earth-Mars mission, but have run into a little trouble when the spacecraft actually reaches ...
1
vote
1answer
60 views

Invariant error plot for an arbitray system of ODE

i try to use InvariantErrorPlot command for a generic system of ODE but i find only example with pre-building equations from Mathematica Packages. This is one of ...
0
votes
0answers
156 views

Mapping the imaginary axis of the complex (s) plane

In this problem is very difficult to analyze the roots s. Firstly I thought that the problem is in precision, but I am not sure. The criterion which I want to obtain is that one of these roots has a ...
0
votes
0answers
60 views

trying to speed up sequential use of NDSolve for large system of ODEs

I don't have a basic example just yet but I'll describe the problem for now. I have a large system (>100) of coupled first order ODEs. Where the output from one solution becomes the ICs for the next ...
3
votes
1answer
185 views

Convergence in NIntegrate vs Integrate

I am faced with this situation that for a certain integration, $\int _0 ^\infty \frac { \tanh (\pi \sqrt{x} )} {\sqrt{x+10} } dx$ - the command Integrate returns ...
0
votes
1answer
68 views

Strange result using NIntegrate with Max and PDF

This is really trivial and I appologise in advance if I have made a stupid error somewhere. I have been chasing down some strange results when using NIntegrate with ...
1
vote
1answer
655 views

Euler's method for system of differential equation

I need to program Euler's method to solve a system of two diffferential equations of first order. Fist, I have programmed the Euler's method for just one differential equation: ...
3
votes
1answer
138 views

How to do this complex integration on the real line?

$m, r$ are parameters in the following integral: Integrate[z Exp[I z r]/Sqrt[z^2 + m^2], {z, -∞, ∞}] How to do this integration directly? The result should be ...
0
votes
1answer
97 views

NDSolve with arrays and Tables of Equations and WhenEvent

I've been toying around with NDSolve for a while, and read through the website. By doing so I discovered that I could use it for vectors and arrays with much pleasure. So I wanted to write a simple ...
0
votes
1answer
108 views

Multiple NIntegral with singularity

I need to compute the following integral FF[0]. When evaluating FF[0] Mathematica simply takes too much time. The problem may ...
0
votes
1answer
164 views

Initial and boundary value errors

NDSolve::ivone: Boundary values may only be specified for one independent variable. Initial values may only be specified at one value of the other independent variable. >> I keep getting this, ...
-6
votes
1answer
138 views

Poincaré section for 3BP

I try to plot the Poincaré section for the 3BP i use the following code that works very well: ...
0
votes
2answers
165 views

How to numerically integrate this integral?

I want to integrate the function funcin (spherical coordinates). ...
1
vote
1answer
127 views

Integrate with FunctionInterpolation differs from NIntegrate

In working with interpolating functions and integration, I noticed a lot of differences between NIntegrate and Integrate in the ...
1
vote
2answers
212 views

intersection between function and plane

I try to solve a ODE system of equations and i want to plot the intersection of the 3D solutions with a costant plane. To do this i use the follow code: ...
3
votes
0answers
104 views

Non-linear integral equation

I'm trying to solve with Mathematica an integral equation. I found this excellent answer (How to solve a non-linear integral equation?) solving with a collocation method a problem which can be ...
11
votes
3answers
884 views

Plot Matlab icon

I started to explore this on a whim and hasn't succeeded yet… Some introduction for the icon is found here but I can't understand it very well. (I admit that, though playing with ...
0
votes
0answers
74 views

Speeding up Integrate and the numerical solution

I want to compute the drag of a plate Below is the code I am using to integrate the function ...
0
votes
1answer
176 views

Infinite expression 1/0.^(1/4) encountered in NDSolve

I am trying to solve very simple system of coupled differential equations using the NDSolve function, but I am getting ...
4
votes
2answers
252 views

How to use NIntegrate when there are symbolic constant coefficients

I would like to numerically integrate an equation such as the one below in which there are symbolic constant coefficients. I used a very simple code but it doesn't work in general, that tried to deal ...
0
votes
0answers
42 views

MathLink wrapper for NIntegrate

I wonder if I can benefit from the Mathematica numerical quadrature routines in a way which integrates completely into an existing C-program. To be more specific I have C-program which does numerical ...
0
votes
3answers
323 views

Solving an Integral Numerically

I have been trying to solve the integral equation below, but cant seem to find a way out of this. Can someone please help me out with suggestion? $f(t)=\int_0^{\infty}\frac{K_1a(t)}{a(t)+K_2}\,dt$ ...
0
votes
0answers
83 views

Can this function be numerically integrated

I need to integrate the function, $f(a,b)=-g(a,b)ln(g(a,b)) = -1/ Pi(( (E^(-((a^2 - E^(-t [Gamma]) [Alpha]^2)/( 1 + 2 [Sigma]^2)))) [Pi] ([Sigma]^2) )/((1 + 2 [Sigma]^2) (E^(2 ...
-1
votes
1answer
192 views

Multiple NIntegrate (again)

Basically I have the same question as here: Multiple NIntegrate but since I don't have enough "reputation" I cannot comment there. I want to solve the following multiple Integral numerically at given ...
0
votes
1answer
226 views

Using NDSolve for three variables as function of time

I want to be able to plot z, and r. Should I be using DSolve, NDSolve or something else? I've tried to clear it up a bit Thanks, Will ...
3
votes
1answer
207 views

Integrate equations of motion when force depends on position

I'm fairly new to Mathematica. I'm trying to test my C++ implementation of a fourth order Runge Kutta method for Newton's equations of motion. I want to test my integrator when the applied force is ...
1
vote
1answer
182 views

Boundary condition is not specified on a single edge of the boundary of the domain

I am attempting to solve a PDE as prescribed in a paper I am reading; however, Mathematica doesn't seem to be happy with the way I am posing the problem. Im my case, r(t,100) is an approximation of ...
0
votes
2answers
211 views

Double integral with variable in the limits

I'm trying to solve this integral: $$ ...
2
votes
1answer
218 views

Numerical Integration with Variable Parameters

So I want to numerically compute the integral of a long complicated expression over a specified domain (in this case an ellipse). I know how to use a Boole function to sample within the ellipse, but I ...
4
votes
1answer
684 views

Is it possible to calculate a Lebesgue integral in Mathematica?

As the title says, I wonder if it is possible to calculate a Lebesgue integral in Mathematica, especially when the domain of integration is $\mathbb{R}^N$, or in other words multivatiate Lebesgue ...
0
votes
0answers
179 views

Using Mathematica to Solve for Electric Field from finite cylinder

For a cylinder that is of length L and radius s and possesses a charge distribution of sigma across only the cylinder part and not the endcaps, I am trying to calculate the Electric-field. The ...
1
vote
2answers
156 views

How to make a single-variable out of a NDSolve solution?

So I've got this line that contains the solution of a partial, non-extraordinary differential equation (because Mathematica doesn't handle extraordinary partial differential equations): ...
2
votes
0answers
127 views

Global vs local adaptive integration: lattice propagator

I was doing the first exercise in the paper Lattice QCD for Novices. This is the expected result: With the default "GlobalAdaptive" method for NIntegrate it threw errors saying that the error had ...
1
vote
0answers
128 views

Choosing PrecissionGoal and AccuracyGoal for decomposed integral with NIntegrate

Assume I have a complicated integral $I$ which can feature all kinds of difficulties like infinite interval, rapidly oscillatory integrand, integrable singularities and numerically almost singular ...
1
vote
1answer
224 views

Error Function Integral (Erf)

Any idea how to solve analytically this integral Integrate[(a Erf[a Sqrt[b/(a^2 + b)] c])/(a^2 + b)^(3/2), a] I tried substitution u=a^2 + b, but it didn't work. ...
1
vote
3answers
630 views

Solve for variable in integral limit when using NIntegrate

I am trying to solve for a variable which appears in the limits of an integral, rather like in the following; ...
0
votes
0answers
94 views

Convergence problems with numerical integration of specific functions

I want to numerically integrate two functions that have several poles inside the integration region. These are the functions: ...
7
votes
1answer
260 views

Is there a way to see the result of NIntegrate's symbolic preprocessing?

NIntegrate can do a number of different types of symbolic preprocessing on the integrand before starting the numerical calculations, including changes of variables. ...
8
votes
2answers
280 views

Why do I get a different value when I change the order of integration?

I think the following two-dimensional integrals should be equal, since they both integrate the function over the half plane defined by $t>\tau$. $$\int_{-\infty}^\infty \mathrm{d}t ...
4
votes
5answers
414 views

NIntegrate::slwcon Problem

I have a problem with numerical integration of this function. Integral value is zero, but NIntegrate[] needs a lot of time to calculate this. Is there any way to ...
4
votes
1answer
224 views

How to choose MaxStepFraction for optimal speed of NDSolve

I'm trying to use NDSolve to solve a 1D Schrodinger's equation, and it seems that MaxStepFraction has huge effect on the ...
3
votes
1answer
300 views

A Bessel & Struve functions related integral

I try to numerically compute this integral and I don't figure out why on earth Mathematica is not able to do it. Is my input correct? Does it possibly have a closed form? ...
0
votes
1answer
274 views

How to use NIntegrate in a function using parameters from a list

I would like to use NIntegrate in a function with some parameters from a list. I simplified my problem for this forum. The list of parameters is as follows: ...
0
votes
0answers
84 views

ParametricNDSolve: Parameter inside UnitStep leads to longer integration times

I have found one of the most annoying issues in ParametricNDSolve. When a parameters occurs inside UnitStep (and Piecewise and If) the integration takes much longer time. I illustrate this by an easy ...
3
votes
1answer
314 views

Question on accuracy and precision of NIntegrate

As a Mathematica newbie, I was testing the accuracy/precision of NIntegrate (9.0.1.0 on Mac) and have obtained a very peculiar result. ...
0
votes
0answers
64 views

PDE syntax question

I am trying to solve a PDE of a reaction diffusion equation, and while I'm getting a decent looking solution, I am also getting inconsistent BC / initial conditions warnings and wanted to check if my ...
0
votes
0answers
126 views

Differential equation with random variable

How can I derive analytically or compute numerically the solution to following differential equation $$ dy/dt = y\cdot X\cdot (y\cdot X - g(y,X))\cdot X $$ where X is a random variable (e.g. from a ...
0
votes
2answers
163 views

NDSolve boundary query / Extracting values from solution

I have a function $O_{2}$ with boundary conditions $O_{2}(r_{o}) = p_{o} $ and $O_{2}(g) = O_{2}'(g) = 0$. I plot it using the solver code below; ...
2
votes
0answers
128 views

Slow integration of interpolating function

I have a large file with irregularly spaced data in the format: ...
5
votes
2answers
201 views

Triple fractional part-related integral

The evaluation with Maple suggests the triple integral is around $1$, but Mathematica tells it's $0.0958758$. When using the code ...