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)

6
votes
2answers
154 views

Different results of the integrate $\int \frac{\sin (\pi s)}{(s-0.5)^2+90} \,ds$?

Following results are given of big different values. I think NIntegrate's one is right but others are not explained. ...
6
votes
1answer
94 views

Why should the spatial derivative order of the ODE *not* exceed two?

Following this question I came across this strange behaviour. Let me define a 1 D interval implicitely ...
6
votes
1answer
151 views

How to do automatic differentiation?

Does Mathematica have any AD functionality or does it only support symbolic differentiation? If not, are there any packages or other third party implementations available? Edit (J. M.) Here is an ...
6
votes
2answers
299 views

Integral does not converge (when it should) [duplicate]

Bug introduced in 10.0 and fixed in 10.0.2 I was looking at integrals like: Integrate[HermiteH[50, x]*Exp[-x^2], {x, 0, Infinity}] which gave me a "does not ...
6
votes
1answer
534 views

Integrating over Bessel Function erroreous? (Hankel Transform)

Bug introduced in 8.0.4 or earlier and persists through 10.0.2. The Hankel Transform is given by Integrate[f[x] x BesselJ[0,x t],{x,0,Infinity}] It is ...
6
votes
1answer
158 views

Problem overloading D[] using UpValues

First I define an entry of UpValues overloading D[] for expressions with head ftest: ...
5
votes
4answers
2k views

How to find the sum all even numbers of this sequence?

I have a sequence $(u_{n})$ $$u_1= 1, \quad u_2 = 2, \quad u_3 = 3, \quad u_{n}= -u_{n-3} + 3u _{n-2} +2 u_{n-1}, \quad \forall n \geqslant 4.$$ I want to list the first $20$ terms of this sequence ...
5
votes
5answers
3k views

Using implicit differentiation to find a line that is tangent to a curve at a point

A few days ago I asked about using differentiation to find a line that is tangent to a curve at a given point. J.M. provided a very elegant way to solve these kinds of problems in Mathematica. Now ...
5
votes
3answers
314 views

How to force Dt[ ] to recognize dependencies

I want to differentiate a function f[] for which I don't have a specific expression. f[] depends on ...
5
votes
3answers
897 views

recursive integration

I am trying to do multiple integrations recursively. For instance, I would like to do the following equation for arbitrary integer $n$: $\displaystyle R_n(t) = \int_0^t \mathrm dt' R_0(t-t') ...
5
votes
3answers
225 views

Derivative of contour in ContourPlot

Sorry to bring this question up again, since there are many similar questions on the site. I use to think this is a easy job to do, because for the worst case I can follow the answers on this site. ...
5
votes
2answers
418 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
3answers
394 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
1answer
262 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 ...
5
votes
2answers
237 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 ...
5
votes
1answer
1k 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 ...
5
votes
2answers
668 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
193 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
1answer
267 views

How do I compute the entropy of the beta distribution?

I tried Expectation[-q*Log[q], q \[Distributed] BetaDistribution[a, b]] and got ...
5
votes
3answers
357 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
475 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
2answers
182 views

Double improper integral

I am trying to compute the following double integral: Integrate[(r BesselJ[0, r])/( 2 (r^2 + z^2)^(1/2)), {z, -Infinity, Infinity}, {r, 0, Infinity}] Which ...
5
votes
1answer
305 views

Lagrange Multiplier

We are asked to maximize and minimize $f(x,y)=4xy$, given the constraint $4x^2+y^2=8$, using the Lagrange Multiplier method. First, I enter the functions f and g. ...
5
votes
1answer
220 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
4k 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
147 views

Visualize transformation used in multiple integration

In class, we will sketch the region $D$ in $R^2$ bounded by the curves $xy=1$, $xy=4$, $y=x$, and $y=x+2$. Then we will use the change of variables given by $u=xy$ and $v=y-x$. We need to show that ...
5
votes
3answers
197 views

Extrema of a function of three variables

To identify and classify the critical points of the function $f(x,y,z)=x^3+xz^2-3x^2+y^2+2z^2$, I used the Hessian matrix method. ...
5
votes
1answer
590 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
154 views

Fitting a Parametric Curve

I apologize if this question is obvious or simple. I'm in high school and have started playing around with Mathematica only recently. For my Calculus Final, I am supposed to use rotation about the ...
5
votes
1answer
271 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
204 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
77 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 ...
5
votes
2answers
238 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
342 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
248 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
701 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
375 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
263 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
2answers
398 views

How can I numerically solve for fractional functions and fractional derivatives?

I would like to plot fractional functions. Say, $f(x)=\sin^{(1/2)}(x)$. By that, I mean that $f(f(x)) = \sin(x)$. Similarly, I can define a half-derivative to be an operator $H$ such that ...
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
378 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
192 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
1answer
114 views

ParametricPlot3D with tube and arrow

I have this: ...
5
votes
1answer
66 views

Integrate gives Simplify::infd “simplified to Indeterminate” message

Using version 10.1 on windows 7 Integrate gives this Simplify::infd Kernel message. What causes it, is this a bug? ...
5
votes
1answer
81 views

Second derivative of SmoothKernelDistribution PDF?

Does anyone know why I can't seem to get a second derivative of this interpolating function f? In this example, h[x] is zero ...
5
votes
3answers
252 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
166 views

Confirming conservation laws for complex valued functions

Consider the nonlinear Schrödinger equation (I would like to do this for a more complicated set of equations, but to gain understanding I'll consider this simplified case) $$A_t+iA_{xx}+i|A|^2A =0,$$ ...
5
votes
1answer
207 views

How to plot the result of this singular integral?

Please I open a new post here after this one : http://mathematica.stackexchange.com/a/59203/10158 Now I want to plot the function $f(a,b)$ as a function of $b$ for different values of $a$ : $a=0.5$ , ...