Questions related to the calculus and analysis branches of Mathematica, including, but not limited to, limits, derivatives, integrals, series, and residues.

learn more… | top users | synonyms (3)

1
vote
1answer
37 views

NIntegrate::inumri: error

I am getting the following error: NIntegrate::inumri: The integrand c (1-0.6 (1-(1+Times[<<2>>])^0.5)-0.4 (1- (1+Times[<<2>>])^0.5)) (1-0.6 (1-(1+Times[<<2>>])^0.5)-0.4 (1- ...
2
votes
1answer
217 views

Existence of an analytical form for integral

Can the following integration be performed with an analytical output? ...
2
votes
2answers
132 views

Symbolic vs numerical integration

This is related to a question that I asked before but didn't resolve quite well. I want to obtain the solution of the following expression analytically: ...
0
votes
0answers
51 views

Product of two Meijer's Function

I want to evaluate an integral $I_1$ defined in $Eq.(1)$ as \begin{align} I_1=\int_{0}^{\infty}\frac{x\exp(-\beta x)K_1(\alpha x)}{1+x}dx\tag{1} \end{align} Where $\alpha\geq0$, $\beta\geq0$, and ...
1
vote
0answers
29 views

Hidden assumption in Mathematica integration [duplicate]

If I do Integrate[Exp[I*omega*t],t], it gives me $\frac{1}{\omega} e^{i \omega t}$, as expected in most of the cases, except when $\omega$ is 0. Mathematica seems ...
0
votes
1answer
36 views

Plotting a function based on complicated integral

I have this function : f[Lambda_] := K Integrate[x Exp[- x^2-Lambda x] HypergeometricU[-Lambda,1/2,(x+ Lambda/2+2)^2],{x,0,Infinity}]; where ...
3
votes
1answer
65 views

Euler-Maclaurin summation

I want to compute asymptotic approximations to partial sum of harmonic series in Mathematica, using Euler-Maclaurin summation formula. ...
0
votes
0answers
62 views

How to evaluate this integral with a bessel function

I have to evaluate this integral, Sin[θ] 2 π BesselJ[0, k Sin[θ] (μ[j] - μ[i])] Exp[2 k^2 σ^2 Cos[θ]^2], wrt $\theta$ between the interval [0, $\pi$]. The ...
0
votes
1answer
55 views

How to calculate the derivative at a point of a non-integrable multivariable function?

I want to integrate this function with respect to θ1 and then take the derivative of that function with respect to n. Since the ...
2
votes
1answer
41 views

Why do Wolfram Alpha and Mathematica give different answers? [closed]

I've just got my hands on Wolfram Mathematica and tried to see if I can figure out how to use it by myself. I've plugged in an integral $\int \frac{dx}{(1-x^2)^{2/3}}$, knowing that the result is ...
2
votes
1answer
105 views

Problem with Volume [closed]

A bug probably. This should give the volume of the unit sphere: vol = Volume[{1. r, th, ph}, {r, 0, 1}, {th, 0, Pi}, {ph, -Pi, Pi}, "Spherical"] but it gives: ...
18
votes
3answers
2k views

How to deal with complicated Gaussian integrals in Mathematica?

As we know, for most Gaussian integrals, we can get the analytical result. Now I have many Gaussian integrals to treat, which have the following general form, ...
16
votes
3answers
360 views

How to efficiently find moments of a multinormal distribution?

Update: Starting from V10.0 the build-in Moment is fast enough for practical use. I have a multinormal distribution with covariance matrix $\sigma$ and zero ...
2
votes
2answers
131 views

Confirmation of a nice closed form

Can Mathematica 10 confirm a nice closed form for ...
1
vote
2answers
110 views

How to get the $0$ result for this limit?

Consider the following series: $$\sigma (\text{x$\_$})\text{:=}\sqrt{\sum _{k=0}^{x-1} \frac{(x! x!) (-1)^{x-k}}{(x-k) (k! (2 x-k)!)}+\frac{\pi ^2}{12}}$$ From the following codes, it seems that the ...
4
votes
6answers
870 views

Randomly Multiply OR Divide by a Random Number?

I want that Mathematica decides randomly if a given variable, gets multiplied or divided by a random number in a given range. The "naive" idea to just use ...
4
votes
1answer
86 views

Finding Euler spiral convergence point for powers other than 2

I'm playing around with Euler spirals, using a variable power to see how the plot changes, $ \left\{\begin{matrix} x = \int_0^s \cos \left ( s^n \right ) ds \\ y = \int_0^s \sin \left ( s^n \right ) ...
1
vote
1answer
59 views

Wrong Limit[…] for inexact expressions

I noticed that Limit can return nonsense when using inexact parameters. In the following code, a is my (exact/inexact) ...
0
votes
1answer
72 views

Integrating high oscillatory function - triple integral - error message

I am working in a triple integral, but the second integral (qQ) and the first too ( aA ) is a oscillatory funcion. The third integral ( os11 ) is like a gaussian, and is easy to integrate. i am ...
9
votes
1answer
129 views

Non-rectilinear integration of InterpolatingFunction

Bug introduced in 7 or earlier and persisting through 10.3.1 [CASE:3487737] I have a 2-d array that I would like to resample into a different coordinate system and integrate along one of the new ...
2
votes
2answers
121 views

A simple question: Input problem

I have a simple problem. First I define the value of variables, then I calculate the expression as the following: Code 1: ...
0
votes
1answer
324 views

Integration of linked PDFs over probability simplex

I'm trying to integrate the following expression: f[a1,p1+p2+p3+p4]*f[a2,p2+p3+p4]*f[a3,p3+p4]*f[a4,p4] Where ...
14
votes
3answers
1k views

How to evaluate the $0/0$ type limit in Mathematica?

When I use Limit to evaluate the limit $$\begin{align}\lim_{k \to 0} \frac{ (k+2)(\alpha^2 - \sqrt{\alpha^4 + k}) + k}{\alpha^2 - \sqrt{\alpha^4 + k} + 2 ...
0
votes
0answers
28 views

How do I numerically integrate the result of NDSolve if it depends on several variables?

I am trying to numerically integrate solution of NDSolve. I am trying to integrate my result of NDsolve in line 97 of my notebook, but am getting the error ...
4
votes
1answer
86 views

Pertubation series of given equations

I would like to define a function that does the following. If I have an equation, for example: $$f(x)+[1-f^2(x)]+f^{\prime\prime}(x)+f(x)\ f^\prime(x)=0$$ and I am given a value for $f$ to expand ...
1
vote
1answer
51 views

The assumption in integration does not work?

I am new to mathematica. Here is the question. I try to do the following integration with the assumption that k1<0, k2<0 and k1,k2 are real number. However, the result includes \sqrt{k1}? ...
4
votes
2answers
325 views
0
votes
0answers
15 views

Integration with IntervalMemberQ

I define an interval and its associated index function. ...
1
vote
2answers
137 views

Is there a way to do a cyclic convolution analytically in Mathematica?

Let's say I want to convolve a Gaussian function Pulse[t_] = Exp[-(t - Mid)^2/(2 sig^2)]; with an exponential broadening function ...
0
votes
1answer
54 views

Integrate over FindRoot solutions

I have a function of bivariate normal PDF and its marginals defined as ...
0
votes
2answers
70 views

Expression which should be real-valued appears to be complex-valued [closed]

I want to plot the answer of the following integral: Λ = 10; Ω = 10; t = 1; Integrate[ω E^(-ω/Λ) Sin[(ω + Ω)t/2]^2/(ω + Ω)^2, {ω , 0, ∞}] The answer of ...
6
votes
2answers
281 views

An integral with a fractional part in 3 dimensions

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

Gradient of interpolated 3D data in mathematica 9

tt = Flatten[Table[{x, y, z, btot[x, y, z]}, {x, -1, 1, 0.1}, {y, -1, 1,0.1}, {z, -1, 1, 0.1}], 2]; ff = Interpolation[tt] Till here it is working fine as it is ...
4
votes
1answer
101 views

Products of Differential Operators

I have a differential equation defined as the product of operators which I want to expand out into a polynomial in powers of $z\frac{d}{dz}$ $\qquad \prod_{n=1}^p(z\frac{d}{dz}+a_n)$ However when I ...
0
votes
2answers
126 views

How can I solve this integral analytically?

I want to solve the following integral: Integrate[( E^(-ω/Λ) Sin[((ω + Ω) t)/2] Sin[(ω t)/ 2])/((ω + Ω)/4) , {ω, 0, ∞}] But unfortunately Mathematica cannot ...
2
votes
0answers
86 views

Why Can't Mathematica Integrate this? [closed]

I am trying to integrate Absolute-value functions from Wikipedia: ...
5
votes
2answers
188 views

Results in terms of a function

Let's assume I have a function $f(x)=\exp(10x^2)$, i.e. f[x_]:= Exp[10*x^2]. When I differentiate $f$ by evaluating the expression ...
1
vote
1answer
46 views

Differential Equations and Unit Notation

What do I get the Syntax::sntxi: Incomplete expression; more input is needed . error when I try to use the math palette for derivatives? What does the dot in 10. Farad Micro Volt mean? ...
1
vote
1answer
92 views

Use Meijer-G function to represent elementary functions

I want to represent these elementary functions: $x^{2}\sqrt{x}$, $\sin{4x}$, and $x\ln{x}$ as cases of MeijerG. What arguments should I give to ...
3
votes
3answers
118 views

Limit n->Infinity of recursive sequence

I have defined a recursive sequence a[0] := 1 a[n_] := Sqrt[3] + 1/2 a[n - 1] because I want to calculate the Limit for this ...
5
votes
5answers
119 views

How to find the derivative value at $(\pi,0)$ for this implcit function $n$ times?

I am trying to take the implicit derivative at $\sin(x+y)+\sin(x)=y$ and substitute $x=\pi$ and $y=0$ at least 6-7 times since I need to find the Taylor series for this function. Since I barely ...
-1
votes
1answer
54 views

Taylor and Maclaurin series [closed]

I need to calculate limit of sin(x)*cot(tanx) at 0 using Maclaurin series I took the usual approach but did not get anything of help.
0
votes
0answers
62 views

How to evaluate this integral

I am trying to evaluate the following integral in Mathematica: ...
0
votes
1answer
43 views

Same integral gives different results with assumption

b = 1 k = 12560 r = 500 res = Integrate[Exp[I*k*x*Sin[o]]/(r + x*Sin[o]), {x, -b/2, b/2}] resab = Abs[res]^2 Plot[resab, {o, -Pi/20, Pi/20}, PlotRange -> Full] ...
0
votes
0answers
15 views

How to solve this system of Partial Differential Equations [5] [migrated]

$${\partial z \over \partial x}f(x,y,z)+{\partial z \over \partial y}g(x,y,z)=h(x,y,z) $$ $${\partial z \over \partial x}F(x,y,z)+{\partial z \over \partial y}G(x,y,z)=H(x,y,z) $$ ...
1
vote
1answer
113 views

DensityPlot showing artifacts [closed]

This is my first post on Mathematica.SE, and I'm still a bit new to Mathematica. I believe I'm probably doing something silly in my code, but I have not been able to find where. I'm trying to plot a ...
0
votes
0answers
33 views

Integration of a Gaussian function in kernel density estimation

I came recently across an odd result of mathematica, which so far I cannot explain. In a standard kernel density estimation (example http://www.ssc.wisc.edu/~bhansen/718/NonParametrics1.pdf) the ...
5
votes
2answers
147 views

Symbolic integration of SphericalBesselJ

Backslide introduced in v10 and persisting through v10.3.1. Consider the following integral ...
3
votes
0answers
80 views

Inconsistent results for symbolic trigonometric integral

I am trying to evaluate (on Mathematica v.9) the integral \begin{equation} \int_0^x [\sin (x) \sin (2 y)-\sin (2 x) \sin (y)]^t \,\mathrm{d}y, \end{equation} where $t$ is an even, positive integer. ...
5
votes
3answers
507 views

Plotting the direction field of a differential equation

I have to sketch the direction field for the following differential equation: $$\frac{dy}{dx}=\frac{-0.02 y +0.00002 xy}{0.08 x-0.001xy}$$ This is the code I used, which gives an incorrect plot: ...