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
1answer
1k views

How to expand a function into a power series with negative powers?

Is there any way to expand this expression a+b(1-Exp[-T/(b c)]/(z-Exp[-T/ (b c)]) (where a, ...
6
votes
2answers
2k 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?
6
votes
2answers
345 views

Force integration to be linear over sum?

Is there an obvious way to force Mathematica to separately integrate the terms in a sum? According to the docs, When part of a sum cannot be integrated explicitly, the whole sum will stay ...
6
votes
1answer
395 views

What is the mathematical meaning behind D[f]?

y = a + b x; I can understand this output of the ordinary differentiation of y w.r.t. x ...
6
votes
1answer
310 views

How do I compute the entropy of the beta distribution?

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

Curvature Application

I found these amazing animations at https://courses.engr.illinois.edu/tam212/avt.xhtml#avt which I think will be wonderfully useful when teaching multivariable calculus next semester. I've started to ...
6
votes
1answer
355 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 ...
6
votes
2answers
281 views

An integral with a fractional part in 3 dimensions

The evaluation with Maple suggests the triple integral is around $1$, but Mathematica tells it's $0.0958758$. When using the code ...
6
votes
3answers
196 views

Speed up derivative evaluation

I'm trying to calculate the normal vectors and the tangent vectors at the discrete points of a suface. For example: ...
6
votes
1answer
407 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. ...
6
votes
4answers
1k views

Can't compute definite integral

Consider a scalar field (in polar co-ordinates), $f(r) = l-r$. Now I want to evaluate the field integral over a circular region of radius $b$, centered at a distance of $x$ from the origin. According ...
6
votes
1answer
91 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 ...
6
votes
1answer
181 views

Defining a function implicitly and calculating its derivatives

I'm interested in using Mathematica's symbolic manipulation to obtain, for a particular function $f$, derivatives of arbitrary order evaluated at zero. Normally I'd use the ...
6
votes
1answer
280 views

Find exponential generating function from the first few terms

The function FindGeneratingFunction computes the ordinary generating function of a sequence, given a sufficient number of initial terms. I have two questions in ...
6
votes
2answers
189 views

Series resulting in “No more memory available.”

Is there any way to force Mathematica to come up with the closed form to Sum[1/(i^18 (i^2 + j^2)), {i, 1, Infinity}, {j, 1, Infinity}] or ...
6
votes
3answers
295 views

Using Grad, Div, etc in polar coordinates

I'm trying to use the newish functions Grad and Div to work with a PDE with circular or spherical symmetry and getting what looks to me like a wrong result. However, I'm not sure of this--perhaps I ...
6
votes
2answers
310 views

Is there a built-in function which detects singularities in a function?

Given a function f[x] and a region M in the complex x-plane, how can I find singularities of f in this region, i.e., issue a ...
6
votes
1answer
195 views

Traces of the level surface $z=4x^2+y^2$

I came up with this method to plot the traces of the surface $z=4x^2+y^2$, in this case for $z=1$, 2, 3, and 4. ...
6
votes
1answer
200 views

Working with a system of differential equations that cannot be solved explicitly

I have to work a lot with three functions $\;o_1(t), o_2(t), o_3(t)\;$ that are solutions to the certain system of differential equations: ...
6
votes
1answer
95 views

Is it possible to find a function from first few terms in the expansion

Is it possible to find a function if first few terms of the expansion is known. For example if I have this series $f(x)=\frac{k^3 x^2}{6}-\frac{k^5 x^4}{120}+\frac{k^7 x^6}{5040}-\frac{k^9 ...
6
votes
1answer
1k views

Convolving/integrating problems

Let's say I want to convolve two functions (f and g), a gaussian with a breit-wigner: ...
6
votes
1answer
242 views

Solving $\frac{dx}{dt} = A \frac{ (1-x)}{(t-t^2)} - \frac{(B*x -C*x^2)}{(t-t^2 )*(t-x)}$

I would like to solve the following equation: $\frac{dx}{dt} = A \frac{ (1-x)}{(t-t^2)} - \frac{(B*x -C*x^2)}{(t-t^2 )*(t-x)}$ I have already posted it here, but still it doesn't work for me (I ...
6
votes
1answer
240 views

Incorrect evaluation of integral involving a DiracDelta, whose argument has infinitely many zeros

Let's say $X>0$ is a random variable with probability density $p_X(x)={\rm e}^{-x}$. Define the random variable $Y=\sin(X)$. From the transformation theorem for probabilities we know that its ...
6
votes
1answer
131 views

Is it possible to circumvent a bug inside SeriesCoefficient?

As far as I can tell, there seems to be a bug in SeriesCoefficient: ...
6
votes
1answer
111 views

Examining the function $f(x,y)=xy(x^2-y^2)/(x^2+y^2)$

Consider the function $$f(x,y)=\begin{cases} xy\dfrac{x^2-y^2}{x^2+y^2},&(x,y)\ne(0,0)\\ 0,&(x,y)=(0,0) \end{cases}$$ I can show that $f_x(0,0)=f_y(0,0)=0$. ...
6
votes
1answer
258 views

Integral of x^p

Can anyone explain why Mathematica does not return a conditional expression that handles the case of p=-1 for Integrate[x^p,x]? ...
6
votes
4answers
381 views

Conditional Expectation — How can Mathematica find a more general closed form?

Mathematica can't find a solution to this Expectation. ...
6
votes
1answer
221 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,$$ ...
6
votes
1answer
164 views

Symbolic integration of elliptic functions

Is there a clever way to integrate products of elliptic functions like $\wp(z;g_2,g_3)$ or $\zeta(z;g_2,g_3)$ in Mathematica? Mathematica seems to be able to integrate functions like ...
6
votes
1answer
116 views

Unexpected Behavior from Integrate with PrincipalValue -> True

Bug introduced in 9.0 or earlier and persisting through 10.3 or later The following code, taken from my answer to 97024, ...
6
votes
1answer
164 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
155 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
2answers
486 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
164 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
323 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
980 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
2answers
567 views

Order of a pole

Is there a simple way to determine the order of the poles of a rational function? I have a difficult function where I need Mathematica to find the poles. It would be interesting to also know what the ...
5
votes
3answers
274 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
500 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
4answers
2k views

Numerical differentiation methods

Is it possible to write code in Mathematica that implements various differentiation methods (like forward, central, extrapolated, etc.)?
5
votes
3answers
507 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: ...
5
votes
5answers
119 views

How to find the derivative value at $(\pi,0)$ for this implcit function $n$ times?

I am trying to take the implicit derivative at $\sin(x+y)+\sin(x)=y$ and substitute $x=\pi$ and $y=0$ at least 6-7 times since I need to find the Taylor series for this function. Since I barely ...
5
votes
3answers
109 views

Differentiate a numerically defined function

My function is f[a_, b_] := NIntegrate[Sqrt[(Cos[t] - a)^2 + b^2], {t, 0, Pi}] I want to calculate g[1,1] where ...
5
votes
3answers
236 views

Limit of Nested Radical $\sqrt{1+\sqrt[2!]{2^2+\sqrt[3!]{3^3+…}}}$

I want to find limit of infinite nested radical $\quad \quad \sqrt{1+\sqrt[2!]{2^2+\sqrt[3!]{3^3+...}}}$ but I don't know how to define this expression in Mathematica. How can I define it and ...
5
votes
3answers
406 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
4answers
207 views

Limit of partial sums involving inverse squares

Consider the finite sum rs[x_, n_] := x/n Sum[n^2/(i + (n - i) x)^2, {i, 1, n}] Is there a way to bring Mathematica to calculate the limit for ...
5
votes
2answers
350 views

Replace rule with function? Derivatives don't evaluate

Say I have an expression (call it expr) involving a function, f[x]. I'd like to be able to evaluate that for a particular choice of ...
5
votes
2answers
246 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 ...