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
138 views

Why is the output of my limit expression an interval?

Why is the output of the limit below an interval? It should be precisely $1$. ...
3
votes
2answers
269 views

integral, explicit form

I want to estimate the ratio of integrals: $$ \frac{\int \frac{4 a T^3}{\frac{4 a T^4}{3}+\Lambda_1 \left(\frac{4}{T^3}+1\right)}}{\int \frac{4 a T^3}{\frac{4 a T^4}{3}+\Lambda_0 ...
3
votes
1answer
195 views

One Integral — Two Results

I compared Mathematica's default integral evaluation to what you get by computing the indefinite integral and then plugging in the limits. (I did this because the latter approach was much faster for ...
3
votes
1answer
82 views

Why does Mathematica list the variables in a double integral backwards? [closed]

I am an advanced novice Mathematica user and have done a fair amount of single-variable calculus things with it. This semester I am using Mathematica for multi-variate calculus and when I tried to ...
3
votes
1answer
271 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? ...
3
votes
2answers
255 views

Integrating special functions

I would like to integrate the following function with Legendre polynomial and Gamma function. I am open to suggestions. ...
3
votes
1answer
996 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 ...
3
votes
1answer
474 views

Evaluating double Integral

While trying to evaluate the integral $\int_{y=0}^{x_2}\int_{x=0}^{min(x_1,y)} n (n - 1) (1 - y)^{(n - 2)}dxdy $ , Mathematica does not seem to yield any results. ...
3
votes
2answers
301 views

How do I represent a system dynamics feedback loop?

In System Dynamics, if I want represent the relationship between Speed and Distances, I create a Flow (Speed) and a Stock (Distance) as you can see in this Insight Maker Sample . Heres an Image of how ...
3
votes
1answer
292 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
114 views

Fredholm integral equation of the second kind with kernel containing Bessel and Struve functions

I need to solve this Fredholm integral equation of the second kind: f[s]+integrate[f[t] K[s,t],{t,0,1}]=s where ...
3
votes
1answer
133 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
107 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
79 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
264 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
84 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
114 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
288 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
2answers
311 views

Algorithm for parts integration

Sorry if this is a duplicate, I've searched how to do this to no avail. What I'd like to do is a function that integrates by parts $n$ times, i.e $$ \int u(x) v(x) dx = u ...
3
votes
1answer
310 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
104 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
517 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
336 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
354 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
55 views

Symbolic integration of elliptic functions

Is there some clever way to integrate products of elliptic functions $\wp(z;g_2,g_3)$ or $\zeta(z;g_2,g_3)$ in Mathematica? Mathematica seems to be able to integrate functions like ...
3
votes
0answers
69 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
0answers
80 views

Derivative of generating function (Example from documentation)

In the documentation for GeneratingFunction, the following example is given under Examples -> Properties & Relations -> Derivative: ...
3
votes
0answers
128 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
133 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
191 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
532 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
317 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
1answer
138 views

Integration of a rational function

I am trying to solve the following rational integral (I got assistance here): $$\int_{-\infty}^{\infty}{\frac{a+x}{b^2 + (a+x)^2}\frac{1}{1+c(a-x)^2}}dx$$ where $\{a, b, c\}\in \mathbb{R}$. I would ...
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
2answers
93 views

Integrate returns unexpected result

Consider the following function $$g(x,y):= \frac{1}{( (1+y)^2+x^2 )( 1+ax^2y^2 )^2}$$, where I assume that $y\geq 0$ and $a\in (0,1]$ is a parameter. When I try to evaluate the integral $\int ...
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
397 views

Quickly differentiate and evaluate a function of several variables

How can I differentiate a function with respect to several variables and evaluate it at the same time ? I want to specify also the variable index that I want to differentiate and the number of times I ...
2
votes
1answer
170 views

Mathematica policy for correctness of results [closed]

Does Mathematica provide any kind of warranty that their calculations are correct? Say I'm running a billion dollar company, and I relied on Mathematica to do calculations, like aerodynamics or car ...
2
votes
2answers
261 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
481 views

Speeding up a slow indefinite integral

I have found an alternative way to express an old solution: ...
2
votes
1answer
943 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
589 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
202 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
199 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
86 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
100 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 ...