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)

3
votes
1answer
281 views

Sum and NSum gives different solutions

I'm working on a Mathematica lab for Calc. 2, and I ran into a problem last night. I was trying to calculate the midpoint approximation of the definite integral of ...
3
votes
2answers
261 views
3
votes
2answers
1k views

How to substitute numeric values in a symbolic Jacobian matrix?

I have a multi-variate function from $\mathbb{R}^n\to\mathbb{R}^n$. Choosing any desired initial vector, we can produce the corresponding function value, which is a vector as follows. The main problem ...
3
votes
1answer
120 views

Gateaux (directional) derivatives and higher order differentials of a functional

I would like to calculate the Gateaux derivative of a functional (i.e. a function depending on functions). A simple example for the Dirichlet functional: $L(u(x))=\int_0^1 \frac{1}{2} (u'(x))^2 dx$ ...
3
votes
2answers
106 views

About the wrong evaluation of an integral

Here is an integral I've been studying in my research and I've just realized that Mathematica $8.0$ is unable to correctly compute it. I have 2 simple questions to ask: $1)$. Is my code below ...
3
votes
1answer
75 views

Problem with limit that requires L'Hôpital's rule to compute

Consider the following limit. Limit[(a - Sqrt[a^2 + x])/(a^2 - a*Sqrt[a^2 - x]), x -> 0, Assumptions -> {a > 0}] Mathematica 9.0.1.0 gives ...
3
votes
1answer
261 views

n-fold symbolic integral in Mathematica

I am trying to compute symbolically a n-fold integral (n is a parameter of a function) over, say, the cube [0,a]^n. My code looks like this ...
3
votes
1answer
83 views

How to calculate the residue of $1/f(z)$ at a numerical approximation to a root of $f(z)$?

The input Residue[1/DirichletL[19,10,s],{s,s0}] gives 0 even when s0 is a root. For ...
3
votes
1answer
112 views

Making time differentials look like the textbook [duplicate]

I need to have time differentials to look like the 'textbook'. My code is Dt[x y^2] /. {Dt[x] -> dx/dt, Dt[y] -> dy/dt} which gives the output ...
3
votes
1answer
281 views

Turn list of edges into a polygon function

I have a list of coordinates that define the edges of a polygon and I would like to get a function defining the area Inside out if it (The polygon is convex and the points are in order) So that for ...
3
votes
1answer
291 views

Integral of Lorentzian yields different results depending on when parameter assigments are made

I'm evaluating the integral of a Lorentzian, which I know equals one. First I define the function and evaluate the integral in two slightly different ways. Surprisingly, I do not get the right answer ...
3
votes
1answer
99 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 ...
3
votes
1answer
500 views

Normal and Tangent of Acceleration in 3D

I need to figure out how to find the Normal and Tangent of acceleration. I know the formula for the tangent of acceleration is $((Acceleration . Velocity)/(Velocity.Velocity))*Velocity$ and the normal ...
3
votes
1answer
317 views

Multi-dimensional integral in the complex plane with poles and essential singularity

I've passed the last week searching a way to numerically integrate this multi-dimensional integral in the complex plane at the poles and avoiding the singularity at z=0: $$ \oint_{C}\oint_{C\ auound\ ...
3
votes
1answer
339 views

Integrating over Bessel Function erroreous? (Hankel Transform)

The Hankel Transform is given by Integrate[f[x] x BesselJ[0,x t],{x,0,Infinity}] It is self-inverse, so ...
3
votes
0answers
60 views

Using Mathematica to help obtain correct analytic formula for logarithm integration

I need to program into my Mathematica code the analytic form for the result of the integral: $$I(a,b)=\int ^1 _0 dx \frac{\ln(x-a)}{x-b}$$ that is valid for all complex $a$ and $b$ (but $\text{Im } ...
3
votes
3answers
310 views

Simplifying expression involving total derivatives [closed]

I'm trying to simplify expressions involving partial derivatives of an abstract functions. To obtain the human readable result, I would like to gather terms of total derivatives and hold on ...
3
votes
0answers
127 views

Is this result wrong because of calculation time? (and more questions about Assumptions)

I am confused with the Integrate Given by Mathematica. First let's see a one-dimensional case: ...
3
votes
0answers
132 views

A triple sum related question

I'm trying to compute the triple sum Sum[1/(i! j! k! ), {i, 1, Infinity}, {j, i + 1, Infinity}, {k, j + 1, Infinity}] but Mathematica doesn't return any value. ...
2
votes
4answers
189 views

How to evaluate this indefinite integrate $\csc(4x)\sin(x)$

I tried to integrate the following integral using Integrate[Sin[x]Csc[4x],x] and I am getting a strange result. $$\frac{1}{8 \sqrt{2}}\left(-2 i ...
2
votes
2answers
512 views

3N-dimensional integral

I'm trying to compute a multidimensional integral with a variable number of dimensions. The integral is as follows: $$ \int d^{3N}\!p~e^{-\frac{\beta}{2m}\vec p^2}. $$ I have tried this ...
2
votes
2answers
309 views

Write a function that returns the logarithmic derivative

How can we write a function that if we input an expression f, it returns the log derivative $\frac{1}{f} \frac{df}{dx}$. We have to use a conditional or pattern test so that the function accepts any ...
2
votes
1answer
160 views

How can I get the limit of a certain infinite product?

A first look at it suggests the fact the limit is precisely $1$. To check that I tried Mathematica, but no output so far. Most probably it converges very slowly. Is there any way to get the limit? ...
2
votes
2answers
190 views

Triple series - evaluation delayed

Trying to figure out if the infinite triple series has a nice closed form. It seems Mathematica is unable to help us here. Numerically, things remain the same, no response. Could you help? ...
2
votes
3answers
166 views

Manipulate fails with f[x]=x^k/(1+x^k)

I have this code: Clear[f] f = Function[x, x^k/(1 + x^k)]; Manipulate[Plot[f[x], {x, 0, 5}], {k, 1, 10}] But nothing draws as I move the slider. Am I missing ...
2
votes
2answers
253 views

How to compute this triple integral? [closed]

In Mathematica how can I compute this integral:$$ \iiint_{D}\sqrt{(1-9z^2)(1-4y^2-9z^2)}\,dx\,dy\,dz$$ where D is the domain: $$D: x^2 +4y^2+9z^2\le1$$ Please I need help!!!
2
votes
1answer
468 views

Speeding up a slow indefinite integral

I have found an alternative way to express an old solution: ...
2
votes
3answers
568 views

How to integrate a function over a 3D planar polygon?

I am trying to integrate a function over a planar polygon in 3D. In 2D, this is fairly straightforward, using either answer from this question (I use the second answer). If we use an equilateral ...
2
votes
1answer
899 views

Laurent series expansion

Can someone share how to find a Laurent series expansion of $$f(z)=\frac{1}{(z^2-1)(z^2-4)}$$ centered at zero on the annular disk $1<|z|<2$?
2
votes
1answer
565 views

How to find all roots and maxima and minima and graph them

I need to find all the roots of this function. I know it is simple to find by hand, however, I wish to learn how to do it so I can apply it later. The problem is that it only gives me one root. ...
2
votes
1answer
196 views

Assumptions with D

I want to give assumptions for the D function. Say I want to calculate $\frac{d|x|}{dx}$ for $x>0$, which is 1. I write ...
2
votes
1answer
198 views

Checking if first derivative is zero at a given point

I have a function V[a_, h_, tau_] := (Sqrt[3]/2) h (Sqrt[3] a + h/3 tau)^2 now I want to check, whether its first derivative is zero in a given point. Say for ...
2
votes
3answers
84 views

Converting HeavisideTheta[]s and Sign[]s functions to a single Piecewise[]

I have the product of some Heaviside functions or sign functions, a simple example could be: HeavisideTheta[q^2 - 4]*Sign[q - 3] Another slightly more ...
2
votes
1answer
85 views

Prove (or check) an expression is positive given constraints on variables?

I'm trying to figure out whether an expression is always positive given positive parameters. When the expression is complicated, I can't do this by eye. Is there any way to make Mathematica prove ...
2
votes
2answers
168 views

Another fractional part related question

Having known the fact that the evaluation of the integrals containing fractional parts often contains errors, I wonder if the integral below possibly evaluates to $0.7$ or this is just a coincidence ...
2
votes
2answers
235 views

integrating slightly complicated expression which includes Abs

I have been waiting for Mathematica to give me something for the following integral, errors welcome, but it has been "running" for almost 30 minutes now. ...
2
votes
1answer
812 views

Lie-Bracket of two vector fields

I'm new to Mathematica (installed couple of hours ago) and I need to compute a few Lie brackets between two vector fields $f$ and $g$. $$ f\left(\mathbf{x}\right) = \left( \begin{array}{c} ...
2
votes
1answer
586 views

Finding residues of multi-dimensional complex functions

Say I have a function $f$ of $n$ complex variables, $\{ z_i \}_{i=1}^{i=Nc}$. And then I want to contour integrate the expression such that for each $z_i$ its an integration on an unit circle about ...
2
votes
2answers
188 views

Calculate information entropy integral in infinite square well problem

In the context of information theory, entropy is a measure of uncertainty of a random variable. In quantum mechanics, the uncertainty principle states that $\Delta x\Delta k \ge 1/2$. The same can be ...
2
votes
2answers
122 views

Limit of integral gives incorrect output

I'm trying to evaluate $\lim_{e\to 0} \, \frac{i}{e}\int_{\pi }^0 \frac{1-\exp (i e \exp (i \theta ))}{\exp (i \theta )} \, d\theta$ with Mathematica 9.0.1.0 on OS X. However, I get "Undefined" for ...
2
votes
1answer
94 views

Use Mathematica to determine the falling law

We have a one-variable equation $\rho(R)$ where ρ = (14656.4+277.526*R^2)/(45.9225+R^2)^{5/2} + 0.370036/(R*(0.25+R)^3) This equations describes the evolution of ...
2
votes
3answers
907 views

Multiple assumptions in integral

I want to add two assumptions, so I can get this probability density function to equal 1, though I can't get a solution. ...
2
votes
1answer
205 views

Integrate[(1 + x/n)^n*Exp[-x], {x, 0, Infinity}]

I evaluated Integrate[(1 + x/n)^n*Exp[-x], {x, 0, Infinity}] . thinking the answer should be approximately $\sqrt{\pi n/2}$. Mathematica gave me ...
2
votes
1answer
128 views

How can I calculate a limit with a free variable?

For example, when I evaluate Limit[(1 + x^n + (x^2/2)^n)^(1/n), n -> Infinity] Mathematica does not output any result. When I evaluate ...
2
votes
1answer
311 views

Time derivative of unit vector in spherical coordinates

Is it possible to take a time derivative of a vector given in some curvelinear coordinate system (i.e. spherical)? Mathematica would need to take into account the time dependence of the basis vectors. ...
2
votes
2answers
76 views

ConditionalExpression Limiting Domain

I came across a problem when evaluating the following integral: $\int_0^t \sqrt{9 x^4+1} \, dx$ Now, when I evaluate that, I get ...
2
votes
1answer
159 views

A question about derivatives (cubics)

Given an equation: f(z) = a (x + I y)^3 + b (x + I y)^2 + c (x + i y) + d Which can be re-written as: ...
2
votes
1answer
449 views

How to take the derivative w.r.t. an arbitrary function?

EDIT: The answer posted by @Jens works for built-in functions only. However I'm starting to wonder if this is related to my Mathematica install somehow. I have version 7. Here is an example in the ...
2
votes
1answer
272 views

How can I handle curve singularity in this NIntegrate integration?

Yesterday I asked a question about the non converging integral. Woods told me that it is due to the function which has a singularity along a line which passes through the integration region. (Why ...