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)

10
votes
2answers
149 views
10
votes
1answer
533 views

Why does Mathematica return a Fourier transform for a function for which it is not defined?

The following function $$g(x) = (1 + x^{1/a} )^a $$ should NOT have a Fourier transform, as far as I am aware, for any real values of $a$ since $g(x)$ is not nice in the sense of decays quickly ...
10
votes
1answer
510 views

Is it possible to find a limit of a sequence given by its recurrence relation?

I need to calculate a limit of a sequence given by its recurrence relation. I tried the following: ...
10
votes
2answers
288 views

Asymptotics of $\frac{\sum _{i=0}^{\lfloor n/2 \rfloor} {2(n-2i) \choose n-2i} {n \choose 2i} {4i \choose 2i}}{2^{3n - 1}}$

I am fairly sure that asymptotically $$\frac{\sum _{i=0}^{\lfloor n/2 \rfloor} {2(n-2i) \choose n-2i} {n \choose 2i} {4i \choose 2i}}{2^{3n - 1}} \sim \frac{2}{\pi n}.$$ I tried ...
10
votes
1answer
460 views

Symbolic Integration along contour: branch cut problem?

Context Following this question on path integrals in the complex plane, having defined again a numerical and symbolic integrator along a path as ...
10
votes
1answer
258 views

Contour Integration along a contour containing two branch points

I need to compute following contour integrations: $$f(u)=\oint_\alpha dz \sqrt{z^3+z+u} \qquad ; \qquad g(u)=\oint_\beta dz \sqrt{z^3+z+u}$$ In which $\alpha$ and $\beta$ are two contours in ...
10
votes
2answers
258 views

Mathematica integration failure - new or old?

Bug introduced in 7.0.1 or earlier, persists through 10.1 Consider the following integral: Integrate[Log[a Cos[x]^2 + b Sin[x]^2], {x, 0, 2Pi}] This takes a ...
10
votes
2answers
485 views

How to take derivative of parameterized coordinate?

Suppose I have a vector in $\mathbb{R}^n$ but $n$ is not known in advance. I want to be able to write functions which operate on the components of that vector, and then I'd like to be able to take ...
10
votes
1answer
274 views

Dirichlet coefficients as limits: wrong

Perhaps I should have included the word "bug" in my question. Here we go There is a bug in this Limit Limit[3^s (-1 - 2^-s + Zeta[s]), s -> ∞] (* 0 *) which ...
10
votes
1answer
144 views

How to predict the degree of the first series coefficient?

Given an expression f that is a function of x and a number x0, what is the least integer ...
10
votes
1answer
319 views

Limit of sequence of functions behaving strange

I'm trying to determine the limit of the sequence of functions $$f_n(x)=\left(\frac{1}{\pi}\arctan(n x) + 1/2\right)^n. $$ I define ...
9
votes
4answers
3k views

How to plot and find the volume of a solid?

How to plot and find the volume of the solid enclosed between the paraboloid z=5(x^2+y^2) and z=6-7x^2-y^2 And the answer of ...
9
votes
5answers
5k views

Notation of partial derivative

I want to write partial derivatives of functions with many arguments. Why is it that when I type f[x,y] ctrl+6 (0,1) it turns out to be bad syntax? The output of ...
9
votes
6answers
796 views
9
votes
3answers
7k views

How to find the nth derivative?

This question is not the same as my last one. How do you find the $n$-th derivative where $n$ is a variable? For example, you can find the nth derivative for a specific $n = 3$ ...
9
votes
2answers
233 views

Volumes of Revolution with detailed diagram

All, In calculus I have to do images such as the following in helping explain technique to students. This one is by rotating $y=\sqrt x$ about the x-axis, an image copied from Stewart's Calculus ...
9
votes
3answers
865 views

Inverting a function in a certain region

InverseFunction works well for globally invertible functions, like f = 2*# + 2 &; InverseFunction[f] ...
9
votes
2answers
841 views

Find closed form expression for series expansion coefficients [duplicate]

Is there a built-in function that will find a general expression for the coefficient of the series expansion of a function? Series will only give the explicit ...
9
votes
2answers
725 views

Why aren't these additions of integrals and summations equal?

I have the following code: Simplify[Integrate[f[x] + g[x], x] == Integrate[f[x], x] + Integrate[g[x], x]] To test: $$\int{\left(f(x) + ...
9
votes
2answers
502 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 ...
9
votes
2answers
719 views

Problem with numerical evaluation of analytically solved integral, solution way off

The following command in Version 9.0.1: N[Integrate[x^50*Sin[x], {x, 0, 1}]] gives $1.4615\times 10^{48}$ which is way off from the correct solution which is ...
9
votes
1answer
9k views

Quick Hessian matrix and gradient calculation?

I am absolutely new to Mathematica and I actually want to try implementing a little optimization method . Long story short assuming I have a predefined two-variable function f(x,y) I want to ...
9
votes
2answers
178 views

Integrating a BesselJ integrand to obtain the same result as Maple 16

I would like to check the following integration: Integrate[y*Integrate[1/x*BesselJ[1,x*Exp[I*π/4]]*BesselJ[1,x*Exp[-I*π/4]],{x,0,y}],{y,0,r}] Mathematica 9.0 is ...
9
votes
1answer
343 views

Laplacian and DiracDelta

It is known that: $\nabla^2 \dfrac{1}{|r-r'|} = - 4 \pi \delta^3(\vec{r}-\vec{r}')$. If I do that with Mathematica, I find: ...
9
votes
1answer
164 views

When to use GenerateConditions -> True

Many functions, usuallly those involving integration, take a GenerateConditions option which often defaults to False, or at ...
9
votes
1answer
297 views

Mathematica: computing a difficult integral

I am trying to compute the following integral: Integrate[Exp[Sum[-((cw λ - b[i])^2/(2 σ^2)), {i, 1, n}]], {cw, 0, 1}] And currently Mathematica outputs ...
9
votes
1answer
381 views

How to calculate this integral? Integrate[BesselJ[0, x - BesselJZero[0, 1]]/x, {x, -Infinity, Infinity}]

I tried to calculate the following integral, but it returned unevaluated. ...
9
votes
1answer
190 views

Wrong Limit with LaguerreL

Bug introduced in 7.0 and fixed in 10.2.0 Limit[(n! LaguerreL[n,-1]/(n^(n+1/4)/Sqrt[2]/E^(n-2 Sqrt[n]+1/2))-1) Sqrt[n], n->∞] Mathematica (wrong) output ...
9
votes
1answer
255 views

Why don't products of Dirac deltas integrate correctly?

Bug introduced in 10.0.0 and fixed in 10.0.1 The integral $\int \int \ \delta(x) \delta(y) \ dx dy=1$ evaluates to 0 in Mathematica ...
9
votes
1answer
294 views

Expanding derivatives of hypergeometric functions

Sometimes Mathematica expresses results of integration or summation in terms of symbolic derivatives of Hypergeometric2F1 function, and cannot further simplify ...
9
votes
1answer
576 views

Symbolic scalar-by-matrix derivative

Let's say I want to calculate the following scalar-by-matrix derivative $$\frac{\partial}{\partial A} \text{tr} \left[(\vec X^T A)^T (\vec X^T A)\right],$$ with $\vec X$ and $A$ being a $n \times 1$ ...
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 ...
9
votes
1answer
231 views

Suspected bug in Integrate

Bug introduced in 9.0.1 or earlier and fixed in 10.1.0 In version 10.0: ...
9
votes
1answer
568 views

Integral with HeavisideTheta takes too long to evaluate

I tried to compute $$\int_{-1}^1 d x_1 \int_{-1}^1 d x_2 \int_{-1}^1 d y_1 \int_{-1}^1 d y_2 \theta(x_1 x_2 + y_1 y_2)\,$$ where $\theta$ is Heaviside's step function, by using ...
9
votes
0answers
238 views

Is there a way to teach integrate new solutions?

I have an integral which I can solve, but integrate cannot: ...
8
votes
3answers
2k views

Using D to find a symbolic derivative

I need to do the following: Define a function Take the derivative of this function and have a look at the symoblic representation Substitute in some values With the bonus that I want to use the ...
8
votes
7answers
1k views

Trying to prove that $x\sin(\frac{\pi}{x})\ge\pi \cos(\frac{\pi}{x})$ for $x\ge 1$

Consider the function f[x_] := x Sin[Pi/x] I want to prove that this function is increasing for $x\ge 1$. This can be done with the first derivative. We have to ...
8
votes
5answers
579 views

simplifying $\frac{\log x^a}{a} = \log x$

If one makes the assumptions $x>0,a>0$, then $\frac{1}{a}\log x^a = \log x$. Thus, in Mathematica, Simplify[1/a*Log[x^a], {a > 0, x > 0}] returns ...
8
votes
2answers
3k views

How do I find line integrals?

For example, how can I calculate $$\int_{\left | z \right |=1}\frac{dz}{z}$$ I know that the answer is $2\pi i$ but how do I do it using Mathematica?
8
votes
3answers
668 views

How can I calculate the perimeter of an equation-defined curve?

This is how the curve looks like: ...
8
votes
2answers
618 views

Weird plot with SphericalPlot3D

Taking the equation $x^2-y^2-z^2=1$ and using ContourPlot3D: ContourPlot3D[ x^2 - y^2 - z^2 == 1, {x, -3, 3}, {y, -3, 3}, {z, -3, 3}] Yields the proper image. ...
8
votes
2answers
661 views

Probability: Calculating a multiple integral

Find the value of $P[\Pi_{i=1}^{10}X_i > C]$ for $C=2,5$, where $X_{10\times 1}$ is a random vector with $10$ dimensional Cauchy Distribution having location parameter $\mu_{10\times 1} = ...
8
votes
2answers
321 views

About an infinite product

I have a curiosity as regards the infinite product below. I wonder why Mathematica v.8.0. says the limit is $1$. This is not true. ...
8
votes
3answers
562 views

Biot-Savart Law/Magnetostatics Solution

I am working on solving the Biot-Savart Law equation for the magnetic field around a charged ring of uniform current density. The expression that Mathematica gives is rather nasty, as expected. Notice ...
8
votes
2answers
272 views

Why does Mathematica say $\int_0^1\int_0^1\int_0^1\frac{1.0}{xyz}\,dz\,dy\,dx=0$?

Mathematica 9 says that $\int_0^1\int_0^1\int_0^1\frac{1.0}{xyz}\,dz\,dy\,dx=0$ and $\int_0^1\int_0^1\int_0^1\frac{1}{xyz}\,dz\,dy\,dx=0$. ...
8
votes
4answers
395 views

Definite integral incorrectly giving a nonreal value

fixed in 10.1 (windows and Mac OS X) In Mathematica 10.0, when I enter N[Integrate[Sqrt[1+x^3],{x,-1,3}]] I get a nonreal value (i.e., the imaginary part is ...
8
votes
2answers
716 views

Double integral over a parallelogram

Problem: Evaluate the following double integral $$\int \!\!\! \int_D (14 x^2+61 xy+42 y^2)^{3} \, dxdy$$ where $D$ is a parallelogram between these four lines: $2x+7y=6$, $2x+7 y=-6$, $7x+6 y=6$ and ...
8
votes
4answers
2k views

How do I obtain the correct double limit?

The command Limit[(Sin[x^2] + Sin[y^2])/ (x - y) /. x -> 0, y -> 0] (* 0 *) I think that Mathematica finds the iterated limit instead of the double ...
8
votes
1answer
834 views

Hankel Transform integrals won't work in Mathematica

I'm trying to do this integral, which is shown on the Wikipedia page on the Hankel transformation: $$\int_0^{2\pi}\mathrm d\varphi\;e^{\mathrm im\varphi}e^{\mathrm ikr\cos(\varphi)}$$ The answer is ...
8
votes
1answer
767 views

Integral of the Sinc product

Let as consider the following integral $$ B_n = \int_0^\infty \prod_{k=1,3,5,\dots}^n\frac{\sin (x/k)}{x/k}dx $$ By definition, Sinc[x]=Sin[x]/x therefore ...