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)

8
votes
2answers
327 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
625 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
895 views

Simplifying the derivative of $|x|$

Context In[855]:= D[Abs[x], x] /. x -> 1 Out[855]= Derivative[1][Abs][1] In[856]:= D[x, x] /. x -> 1 Out[856]= 1 Question Why is ...
8
votes
2answers
278 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
2answers
801 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
882 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
834 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 ...
8
votes
1answer
391 views

ReplaceAll[] and Limit[] don't give correct results for this expression under extreme variables [duplicate]

Possible Duplicate: Funny behaviour when plotting a polynomial of high degree and large coefficients 1/x^2 + (3 + x)/(6 (1 - Exp[x] + x)) ——This is a ...
8
votes
2answers
283 views

Make Sinc'[0] return 0 instead of Indeterminate

I want Sinc'[0] to return 0, but instead it returns Indeterminate. I've tried Unprotect[Sinc] Unprotect[Derivative] Derivative[1][Sinc][0] ^= 0 But it doesn't ...
8
votes
1answer
205 views

Strange result of parameter-dependent definite integral

Inspired by this question on Math.SE, I wanted to see what Mathematica could do with the integral $$ \int_0^{\pi}e^{k\cos t}\cos(k\sin t)\,dt. $$ Thus, I entered into Mathematica (tried on both 10.3 ...
8
votes
1answer
269 views

Strange bug: Product[1/(1-1/2^(n-1)),{n,2,Infinity}] = 0?

This infinite product is correct: Product[1/(1 - 1/2^n), {n, 1, Infinity}] (* 1/QPochhammer[1/2, 1/2] *) Numerical value of this result is not equal to zero ...
8
votes
2answers
243 views

Can your Mathematica do the following integral?

Backslide introduced in v10 and persisting through v10.3.1. I was trying to do the following integral in Mathematica: ...
8
votes
1answer
224 views

Curious solutions of x' = Sqrt(x), x(0)=4

Consider the solution of $x'=\sqrt{x}$, $x(0)=4$ using DSolve. ...
8
votes
2answers
866 views

Symbolic integration in real domain only ( assumptions and ComplexExpand don't work)

Integrate[m^2/((x - m^2)^2 + y^2), m] mathematica gives me a complex-valued reuslt, but maple 17 gives me what I want. I tried using assumptions, but it doesn't ...
8
votes
1answer
658 views

Limit not evaluating

Why is the input Limit[IntegerPart[Sin[x]/x], x -> 0] not being evaluated?
8
votes
1answer
689 views

Is there any automatic differentiation package?

I'm wondering if an automatic differentiation package exists for Mathematica. This is what I mean by automatic differentiation.
8
votes
1answer
146 views

LogisticSigmoid residue bug?

Bug introduced in 10.0 and fixed in 10.2.0 On version 10.0 for arm: Residue[LogisticSigmoid[z],{z,I Pi}] gives 0 not 1. ...
8
votes
2answers
919 views

Bug in Integrate for Mathematica

Consider the following: ...
8
votes
2answers
808 views

Numerical contour integrations in the complex plane - contour deformation gives different answer for analytic kernel

I am trying to do a contour integration in Mathematica numerically. In particular, I'm checking the identity: $$ H_m^{(1)}(z) =\frac{i^{-m}}{\pi}\int_{-\pi/2 + i \infty}^{\pi/2 - i \infty} \exp[i m ...
8
votes
0answers
84 views

Slot number cannot be filled

I got a strange message with the derivative of a two-argument function defined by a definite integral. For example if "f" is defined as: ...
8
votes
0answers
100 views

Strange kernel segmentation fault in Integrate

Bug introduced in 8.0 or earlier and persisting through 10.4 Define tointegrate = D[V[y1[t], y2[t]], y2[t]] D[V[y1[t], y2[t]], {y1[t], 2}] and try to find its ...
8
votes
0answers
73 views

Finding simplifying substitutions for an integral involving limits and integrand

[The following is based on a William Lowell Putnam Mathematical Competition problem.] Consider the definite integral: $I = \int\limits_2^4 \frac{\sqrt{\log (9-x)}}{\sqrt{\log (9-x)}+\sqrt{\log ...
7
votes
3answers
1k views

Region bounded by x^2+y^2=1, y=z, x=0, z=0, in first octant

I need to draw (pencil and paper) the region bounded by $x^2+y^2=1$, $y=z$, $x=0$, and $z=0$ in the first octant. So the first assistance I asked of Mathematica is: ...
7
votes
5answers
2k views

Find function inverse

I'm trying to find the inverse of a function: (30*x^2 (1 - x)^2) (* where 0<x<1 *) I tried all the following options: 1. ...
7
votes
3answers
399 views

Generating a polynomial that's accurate to within an error of no more than 1/10^5

I'm currently stuck on a question for class that asks... "Find a polynomial p[x] that you can use to calculate 6 ArcTan[x] to ...
7
votes
3answers
302 views

How to calculate this sum?

I want to find the sum $$S=f\left(\dfrac{1}{2012} \right) +f\left(\dfrac{2}{2012} \right) +\cdots + f\left(\dfrac{2011}{2012} \right), $$ where $$f(x) = \dfrac{4^x}{4^x + 2}.$$ I tried ...
7
votes
4answers
470 views

How to find a tangent line with 2 points of tangency for a curve?

Say I have a function like this: f[x_] := 4 x^4 - 9 x^3 - x^2 + 10; Plot[f[x], {x, -1, 2}] It's obvious that there's a tangent line with 2 points of tangency ...
7
votes
3answers
364 views

Where did I go wrong with my implementation of the trapezoidal rule?

One method for doing quadrature, called the trapezoidal rule, improves accuracy by connecting the points on the curve corresponding to the points of subdivision with line segments, forming trapezoidal ...
7
votes
2answers
377 views

How much time should one give Mathematica for an integral evaluation?

Sometimes when I do integrals in Mathematica (M), it keeps thinking and thinking and I have no idea what is going on inside M. For how long should one wait or how does one know whether M has not got ...
7
votes
2answers
416 views

Need help with the Lagrange multiplier method

Well I am new to Mathematica and got really stuck solving this problem. In class, I was ask to use the Lagrange multiplier method to find the maximum and minimum value of $f(x,y) = x^2+y^2$ which ...
7
votes
3answers
401 views

Negative integral of a positive function

fixed in 10.1 (windows) For a parameter $t\in (0,1)$ $Assumptions = t ∈ Reals && t > 0 && t < 1 I define an obviously positive function ...
7
votes
4answers
1k views

Mathematica gives wrong answer for an integral

When I execute the following, Integrate[ Exp[-w^2 + I w^3], {w, -∞, ∞}] I get (2 E^(2/27) BesselK[1/3, 2/27])/(3 Sqrt[3]) ...
7
votes
2answers
413 views

Check homework integration in Mathematica

Given in my homework I have to compute (by hand) $$\iint\limits_{x^2+y^2\leq 1}(x^2+y^2)\,\mathrm dx\,\,\mathrm dy.$$ My solution so far: Let $f(x,y)=x^2+y^2$ and $K=\{(x,y)\in\mathbb{R}^2:x^2+y^2\leq ...
7
votes
2answers
8k views

Simple ways to evaluate a derivative at a point?

The contrast in behavior between, say, f[x_] = Sin[x^2]; f'[2] vs. u[x_, y_] = Cos[x + y^2]; has always bothered my ...
7
votes
3answers
357 views

Limit of a real function

Evaluating Limit[2^x Log[x], x -> 0, Direction -> -1] (right limit) gives correctly the answer -∞ but the ...
7
votes
2answers
2k views

How to change coordinates of integration?

Is there some built-in routine, some easier method to change variables of integration, without/before solving the integral? Say I have ...
7
votes
2answers
283 views

How to simplify the integration result of the following program?

In the following program, the fec[x] is a function depending on $x$ while q3[t] and q5[t] ...
7
votes
1answer
1k views

How to find the smallest root

I have a continuous, differentiable, monotonic, bounded function called F[t]. If t -> Infinity then ...
7
votes
1answer
137 views

Why is RegionMeasure so slow when calculating intersection area of a 2D and a 3D object?

If you think this description is too long, you can read the problem directly I know normally when one wants to calculate an region, this guide is useful. However, when it comes to calculating an ...
7
votes
1answer
255 views

How to compute the residue of $e^{z-\frac{1}{z}}$ at z=0?

I've tried the following but it didn't work: Residue[Exp[z - 1/z], {z, 0}] not even this: Residue[Exp[1/z], {z, 0}] ...
7
votes
1answer
157 views

Old fashioned region method?

Just stumbled across this idea in http://users.rowan.edu/~hassen/Mathematica/Volume%20III/Chapter%2015.pdf. ...
7
votes
1answer
94 views

Integrate yielding a ConditionalExpression but I don't think the condition is necessary

Suppose I take the PDF of the LogNormal distribution with parameters m and s evaluated at x. I obviously get an expression involving m. I now want to integrate that expression not with respect to x ...
7
votes
2answers
331 views

Integrate fails for generic parameters but works for any specific values (v. 10.0.x bug)

Bug introduced in 9.0 and fixed in 10.1 I was trying to get the following integral ...
7
votes
3answers
335 views
7
votes
1answer
521 views

Mathematica 9 can't integrate this function but earlier versions could

Integrate[ ArcTan[x]/(1 + x) Log[(1 + x^2)/2], {x, -1, 1}] I used Mathematica 9.0.1 on Windows7 32bit, Mathematica 9 cannot compute this, but Mathematica 8 gives ...
7
votes
1answer
329 views

Series expansion for small ratio of variables

I have a messy expression of variables that I would like to simplify under the assumption that certain ratios of the variables are small. For example, consider $\sqrt{x+y}\;$ expanded for small ...
7
votes
2answers
939 views

Integration with vector coefficients

I asked this same question in Mathematics, and it was suggested I might try here. I'm more comfortable with Maple, but if I can get Mathematica to do what I'm after, so much the better. Basically ...
7
votes
1answer
119 views

How to calculate residues using different branches of the logarithm

My goal is to be able to compute the residues of something like $z^{\alpha} R(z)$ where $0 < \alpha < 1$, $R(z)$ is rational, and the exponential can be chosen with any branch. I found a way to ...
7
votes
2answers
525 views

Integrating a function over a surface integral

From a first principles bandstructure calculation I get an energy scalar field in three dimensions $E(x,y,z)$. It's now easy to plot a constant energy (contour)-surface for dedicated values ...