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)

5
votes
2answers
389 views

How can I calculate the limit without using the L'Hopital's rule

I need to prove this limit without using the L'Hopital's rule: $$\lim_{x\to 0} \frac{(1+a\,x)^{1/4} - (1+b\,x)^{1/4}}{x} = \frac{a-b}{4}$$ How can I do it in Mathematica?
5
votes
2answers
347 views
+500

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} = ...
5
votes
3answers
377 views

Bug in complicated Limit

The code Limit[Log[2 - Sin[x]*Cos[x]], x -> Infinity] outputs Interval[{0, Log[3]}] in Mathematica 10.0.2.0 . It should ...
5
votes
2answers
1k views

How to find the non-differentiable point(s) of a given continuous function?

For example, the non-differentiable point of the function $f(x)=|x|$ is at $x=0$. How to find the non-differentiable points of a continuous function that is defined numerically?
5
votes
2answers
618 views

How to receive the result of the integral $ \int\frac{1}{x}dx$ is $\ln|x|$?

When I input \[Integral]1/x \[DifferentialD]x in Mathematica, I get Log[x]. How can I make the result of the integral ...
5
votes
2answers
186 views

Unexpected Integration Constant [closed]

Why does this integration Integrate[15-30x+6x^2-1/(x+5),x] return 700 as integration constant? If we make a slight modification ...
5
votes
2answers
2k views

How do I evaluate a double or triple integral over a region?

Say I need to evaluate the integral $\iiint_W f(x,y,z) dx dy dz$ and $W$ is a region given to me like $W = \{ (x,y,z) : 1 \leq x^2 + y^2 \leq 4, 1 \leq z \leq 5\}$. I don't how to do this with a ...
5
votes
3answers
347 views

ComplexInfinity/Indeterminate error when evaluating derivative

Motivation: I want to find the coefficients for the polynomial that is obtained when one adds together the first $n$ natural numbers to the power of $a$; that is, when you consider $1^{a} + 2^{a} + ...
5
votes
1answer
415 views

Finding ranges of a parameter for which a function is always positive

I have a complex function of a single variable expressed in analytical form, which also depends on a parameter. I would like to have Mathematica show me for which values of the parameter the real part ...
5
votes
1answer
309 views

How to solve this probability symbolically or numerically?

I am trying to calculate the following probability $$\mathbb{P} \big(\sum_{i=1}^{m} (A_i + S_i) \le L < \sum_{i=1}^{m+1} (A_i + S_i) \big)$$ where, $$A_i \sim \exp(\lambda), \quad S_i \sim ...
5
votes
1answer
211 views

How to find this limit correctly?

How to find the limit Limit[n*Sin[2*Pi*Exp[1]*n!], n -> Infinity] ? Mathematica 10 outputs ...
5
votes
2answers
3k views

How to find (numerical) value of a derivative at point?

I have the following function: f[0, 0] = 0 f[x_, y_] := Exp[-(x^2 + y^2)^(-1)] How do I find its partial derivatives at any given point, including $(0,0)$? This ...
5
votes
1answer
452 views

Asymptotic expansion

I wanted to expand a function of $x$ about $x=\infty$ and see its coefficients as a function of the parameters $m,n,q,y$ and I wrote this - but it didn't work! It gave me back a very complicated ...
5
votes
1answer
1k views

Calculating volume between two surfaces of revolution

Once again, I'm very new to Mathematica and winging it as I go. I probably don't have the math education at this point to even justify buying it in the first place, but that's off topic. I wanted to ...
5
votes
1answer
200 views

Simplification of integrals depending on a parameter [duplicate]

Assuming[Element[n, Integers], Integrate[Sin[x]*Sin[n*x],{x,0,Pi}]] returns 0, which is obviously wrong for n=1. ...
5
votes
1answer
234 views

How can I Integrate large trigonometric expressions more efficiently?

I have a requirement to integrate some large symbolic expressions containing products of trigonometric functions. A single integration is averaging 7 minutes, and I have 2628 expressions to integrate, ...
5
votes
1answer
106 views

When to use GenerateConditions -> True

Many functions, usuallly those involving integration, take a GenerateConditions option which often defaults to False, or at ...
5
votes
2answers
220 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 ...
5
votes
3answers
327 views

Integration gives ComplexInfinity when it shouldn't

A simple example of my problem is that Assuming[k ∈ Reals, Integrate[Exp[l]/((1 + Exp[l])*Exp[k*l]), l]] /. {k -> 2} gives ...
5
votes
1answer
227 views

Definite Integral over a path

What does Integrate[f[z], {z, a, b, c, d}] exactly calculate? Is it $$\int_a^b f(z)\, \mathrm{d}z +\int_b^c f(z)\, \mathrm{d}z +\int_{c}^d f(z)\, \mathrm{d}z ...
5
votes
1answer
631 views

Integration over region given by inequality [duplicate]

Sometimes we need to integrate over some bounded region given by inequality/inequalities. Consider the following simplest example of area of an ellipse (LaTeX code): $\int_S dx dy$, where $S = \{ ...
5
votes
2answers
351 views

How can I get the inverse of my function?

The original function is $$f=\frac{1}{e^{\frac{\phi }{k t}}-1}-\frac{m+1}{e^{\frac{(m+1) \phi }{k t}}-1}$$ I want to express $e^{\frac{\phi }{k t}}$ about $m$ and $f$, so I tried: ...
5
votes
1answer
244 views

Prove an identity in quantum harmonic oscillator

Problem: In the context of quantum harmonic oscillator the eigenfunctions are given by: $ u_n(x) = (N_n/\sqrt{b}) H_n(x/b) \exp\left[-x^2/(2b^2)\right] $, where $N_n$ is the normalization factor: $ ...
5
votes
1answer
2k views

finding poles of a function

Is there a command to find the poles of a function $f=f(z)$? example: let $$f(z) = \frac{1}{z^2-1}$$ then we know that the poles are at $z=\pm 1$ but is there a special command in mathematica to do ...
5
votes
2answers
334 views

Mathematica cannot calculate a limit

When I evaluate Limit[E^(-n)*Sum[n^k/(k!),{k,0,n}], n -> ∞] Mathematica gives me the result ...
5
votes
1answer
1k views

Symbolic integration in the complex plane

Context While answering this question, I defined (symbolic and numerical) path integrations as follows ...
5
votes
2answers
662 views

Mathematica 10 cannot solve definite integral

Bug introduced in 10.0.0, solved in 10.0.2 Mathematica 10 fails to solve the following integral, saying that it does not converge. ...
5
votes
2answers
181 views

Why does Integrate return a solution that is not defined at a particular point when it actually is well defined at that point?

I am trying to compute Integrate[Sqrt[x^4 + (y - y^2)^2], {x, 0, y}] Mathematica 8 gives ...
5
votes
3answers
236 views

Calculate integral for arbitrary parameter n in infinite square well problem

I'm continuing[1,2] the study of an infinite square well in the context of quantum mechanics. Ultimate goal is to calculate the product $\Delta x\Delta k$, for various eigenstates, that is for ...
5
votes
1answer
145 views

Integration over a convex combination of a region: $\int_{\Omega} (w_1 z_1 + w_2 z_2)^{1-\sigma} d (z_1, z_2)$ where $\Omega = \{ z_1 + z_2 = 1\}$

Take $w,z\in R^{n}$. I am interested in integrating (as generically if possible) $$\int_{\Omega}(w \cdot z)^{1-\sigma} d z$$ Where the domain of $\Omega$ is $1$ dimension, and includes the convex ...
5
votes
2answers
251 views

Working with abstract vectors

I often need to compute derivatives or integrals involving N-dimensional vectors (where the dimension could be equal to 2 or 3 but is not particularly relevant for the sake of the derivation). The ...
5
votes
0answers
88 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 ...
5
votes
0answers
191 views

Integrate wrong for absolute value of trig function

I was trying to get $\int_0^1 \lvert \cos(2 \pi k x) \rvert \,\mathrm{d}x$ for $k \in \mathbb{Z}$, and was surprised by the result (using Mathematica 10.0.1.0): ...
5
votes
0answers
194 views

Strange behaviour of MMA in derivatives of some standard functions

There are some peculiar things to be discovered in derivatives of some standard functions in MMA: Strange behaviour Example 1: Abs We have ...
5
votes
0answers
101 views

Convoluting inverse square root with Gaussian

I would like to convolute the inverse square root on the interval [0,inf] with a Gaussian function, like so: ...
5
votes
0answers
178 views

Strange Integrate behavior (a bug!)

The following two calculations should give the same result. After all, integration is a linear operation. I have pasted the code below in case you want to play with it. ...
5
votes
0answers
126 views

Calculating a limit with a result that is discontinuous in the parameters

The following limit is left unevaluated (Edit: added the assumption that $\epsilon$ is real thanks to the comment below): ...
4
votes
4answers
618 views

Visualizing Line Integrals

I have a plane curve $C$ described by parametric equations $x(t)$ and $y(t)$ and a function $f: \mathbb{R}^2 \rightarrow \mathbb{R}$. The line integral of $f$ along $C$ is the area of the "fence" ...
4
votes
2answers
1k views

Chain rule while differentiating

I am trying to find the derivative of a function defined in polar coordinates with respect to $x$ and $y$. My function is defined as follows: $ v_x(r, \theta ) = v_r \cos (\theta ) - v_{\theta }\sin ...
4
votes
2answers
378 views

Define product derivative

How do I define the $n$th product derivative of a function in Mathematica? The first product derivative $f^\ast$ of a function $f$ is $$ f^\ast(x)=\exp\left(\frac{f^\prime(x)}{f(x)}\right) $$ The ...
4
votes
3answers
380 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: ...
4
votes
2answers
120 views

Integration of a Rational Function returns RootSum[]

I am trying to solve an integral given below: Integrate[r^2/(-α r^3 + r^2 - 2 m r + Q^2), r] but since coefficients of this cubic polynomial are as parameters, ...
4
votes
2answers
224 views

Limit not giving expected result

I am taking the limit Limit[Sin[π Sqrt[4 n^2 + n]], n -> ∞] the returned answer is Interval[{-1, 1}] I think the right ...
4
votes
1answer
853 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 ...
4
votes
3answers
1k views

Definite integral takes a very long time

I'm trying to integrate the following function: $$ r_e=\int_0^z\frac{dz}{\sqrt{\Omega_m(1+z)^3+\Omega_\Lambda}} $$ where $\Omega_\Lambda=1-\Omega_M$ and both $\Omega_\Lambda$ and $\Omega_M$ are real ...
4
votes
2answers
118 views

Well-defined symbolic integral leading to ConditionalExpression

I would like to determine a closed-form expression for the following symbolic integral $$ \int_{-1/2}^{1/2} \!\!\!\! \mathrm{d} x \int_{-1/2}^{1/2} \!\!\!\! \mathrm{d} y \, \frac{1 + b x + c y}{1 + e ...
4
votes
2answers
240 views

Using Mathematica to confirm Bernoulli's inequality

I have several challenges that I want to confirm is true. I have chosen this one because it is rather simple (proof by induction). There are times when I do not want to spend ages trying find proofs. ...
4
votes
2answers
212 views

Can Mathematica (or its extensions) do integration following Risch algorithm?

I wonder whether there are option for indefinite integration in Mathematica that alow to choose the algorithm? Is there an option to use this algorithm in Mathematica?
4
votes
1answer
235 views

How do I compute the entropy of the beta distribution?

I tried Expectation[-q*Log[q], q \[Distributed] BetaDistribution[a, b]] and got ...
4
votes
2answers
870 views

Sum of series with log in each term [closed]

I was solving recurrence relation of Introduction to Algorithms by CLRS, 3rd. edition. Problem 4-3 (i) $$ T(n) = T(n-2) + \frac{1}{lg \; n} $$ I tried few ways, like expending with iteration method. ...