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)

0
votes
0answers
94 views

Complex Convolution

I am attempting to integrate a convolution variable using the following code. However, the program is taking too long to complete the integration. Does anybody have any coding tips that may make ...
0
votes
1answer
86 views

Defining a function that differentiates another function with rule replacement

This question is related to the question I posted here. In that question, I thought I have simplified my original problem into an equivalent concise version, but I found that it is not. So I decided ...
-1
votes
1answer
130 views

Finding the maximum of a gradient vector

Im trying to find the maximum of my gradient vector G[x,y], I've tried several options including FindMaximumValue, FindMaximum etc. but i couldn't find it. The full function is shown below, any help ...
1
vote
1answer
163 views

Plotting a multivariable function [closed]

m[x_, y_] := 0.9*Exp[-((x - 1)^2 + (y - 1)^2)] + 0.5 Exp[-(3^2 ((x - 2.5)^2 + (y - 1.5)^2))] I'm trying to create a 3d plot for the function above. I tried ...
0
votes
1answer
80 views

Forcing an integral to be solved in separate terms

Using Mathematica I want to solve an integral over a function that contains both determined and undetermined parts which looks like this: ...
2
votes
2answers
235 views

Finding Expectation of function of a Log-normal distribution

Say $Y=g(X)$ and $p_X = \frac{e^{-\frac{(\mu -\log (x))^2}{2 \sigma ^2}}}{\sqrt{2 \pi } x \sigma }$ is Log-normal density function: [Wiki] Find E[Y]? Since $E[Y] = \int_0^\infty y f_Y \ dy = ...
3
votes
1answer
225 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$ ...
1
vote
1answer
115 views

Generate conditions seems to not work [closed]

I am trying to compute the following integral Integrate[E^(I*k*Omega*t), {t,0,T}, GenerateConditions->True] for which Mathematica returns ...
0
votes
1answer
97 views

Integrate and NIntegrate yield different results for double integral

Evaluating a double integral with bivariate normal distribution yileds widely different results depending on the method used. I define a bivariate normal distribution with {10, 3} and {8, 1.5} as ...
1
vote
3answers
172 views

Answer when integrating by partial fractions

Evaluating $$\int \frac{x^2-x+5}{(x-2) (x-1) (x+3)} \, dx$$ Version 9.0.1 gives the following answer $$\frac{1}{4} (-5) \log (1-x)+\frac{7}{5} \log (2-x)+\frac{17}{20} \log (x+3)$$ It seems to ...
4
votes
2answers
234 views

Function given exact arguments returns hugely different value than it returns when given equivalent inexact arguments

I was trying to compute the probability that a coin is from a particular underlying distribution given that a particular set of tosses was observed. (I know this can be done in a different way, but I ...
1
vote
1answer
101 views

Using Mathematica to find Derivatives

I am given a system first order differential equations: $x'=y$ and $y'=6x^2-a/2$, where $a$ is a constant and $'$ denotes $t$-derivatives. I then make the substitution $(x,y)=(x_1y_1,y_1)$. This ...
3
votes
5answers
121 views

Using Mathematica to find a Derivative

I am given a system first order differential equations: $x'=y$ and $y'=6x^2-a/2$, where $a$ is a constant and $'$ denotes $t$-derivatives. I then make the substitution $(x,y)=(x_1y_1,y_1)$. This ...
8
votes
8answers
1k views

Finding maximum or minimum of implicit functions

is there any built in function that can be used to find maximum or minimum of implicit functions? For example, if we have the equation $$x^2 + y^2 = (2 x^2 + 2 y^2 - x)^2,$$ then we can visualize the ...
0
votes
0answers
30 views

How to calculate this derivative? [duplicate]

Can Mathematica do the following simple calculation? $$F(x_0, x_1, x_2, ...) \equiv \sum_n x_n^2$$ I want to do the symbolic calculation: $$ \frac{\partial F}{\partial x_m} = \sum_n 2 x_n ...
1
vote
1answer
76 views

Integrating a function forcing parameters to be real and positive

I'd like to integrate the function $f(x)=x^{a-1}e^{-b x}$ over the interval $[0,\infty)$. f[x_] := x^(a - 1)*Exp[-b*x] Then: ...
0
votes
0answers
60 views

Computing integrals with differential operator

I've to compute this expression $$ \hat{H} = \frac{1}{4}g_2\int d^3R\int d^3r\ \bar{\Psi}(\vec{R}+\frac{\vec{r}}{2})\bar{\Psi}(\vec{R}-\frac{\vec{r}}{2})\left[ \delta(\vec{r})\nabla_{\vec{r}}^2 ...
-3
votes
1answer
69 views

Finding integral bounds

i have this integral shown below is equal to 1, and i need to find "a" on mathematica, but I'm not sure how. integral (from 0 to a) sqrt((-50.8938 sin(8.4823 t))^2+(4-11.3097 sin(11.3097 t))^2) dt = ...
0
votes
1answer
68 views

Integrating function can't be epressed? [closed]

Hi I have a function in a form of $$\frac{a}{x}exp(\frac{b}{\sqrt{3}-x})$$ The code is Integrate[a/x*E^(b/(3^0.5 - x)), {x, 1, 3^0.5}] But it can't work ...
1
vote
2answers
137 views

A regularized hypergeometric function related question

I'm interested in finding a way (if possible) of expressing this specific value of the regularized hypergeometric function in terms of known constants. How might I use Mathematica to check this ...
2
votes
3answers
307 views

Marking specific points on a 3D curve

Here is the curve below. (written as a position vector) r = {Cos[7 Pi*t], Cos[6.2*Pi*t], 2.5*t} ParametricPlot3D[r, {t, 0, 2}] I want to know how I can mark ...
0
votes
1answer
218 views

Is this function convex?

How can I determine convexity of the function f = Log[ x, 1 + (x^a - 1) (x^b - 1)/(x - 1)] with the parameters $a,\,b$ belonging to the interval $(0,1)$ in ...
1
vote
2answers
342 views

How to make Mathematica use the chain rule?

Lets say I have the following PDE: $$x^2 u_{xx} - u_{yy} + u_y = 0$$ And I have the following change of variables: $$ s(x,y) = x e^y \, \, \, , \, \, t(x,y) = x e^{-y}$$ How can I use Mathematica ...
0
votes
2answers
60 views

Problem with using the chain rule [closed]

I have the functions w[x_] := x*Cosh[x] and x[a_] := n/k*a. I take the derivative ...
0
votes
0answers
67 views

Evaluating an indefinite integral at the boundaries

I have the following expression i need to integrate over $\alpha$ from 0 to 1 and $x$ from $x1$ to $\infty$ ...
0
votes
1answer
97 views

Help with Computation-heavy triple integration

I'm Trying to use Mathematica to symbolically integrate a function with respect to x, y, and z. For some reason, this computation either takes about, 60s and returns the input, or takes many hours ...
1
vote
0answers
87 views

Help wanted in evaluating an integral

Trying to evaluate an integral. To be specific, I'm trying to find the expectation value of $x^2$ for the ground state of hydrogen. Here is the code I used: ...
0
votes
1answer
52 views

Derivative of the composition of two functions

Background, this seems to work (This is manually transcribed from the program, beware of errors): J1[f_, i_] := DifferenceDelta[f[i], i]^2 D[J1[f + t v, i], t] ...
5
votes
2answers
229 views

Working with abstract vectors

I often need to compute derivatives or integrals involving N-dimensional vectors (where the dimension could be equal to 2 or 3 but is not particularly relevant for the sake of the derivation). The ...
0
votes
2answers
253 views

How to calculate an improper integral

How can I find the range of the parameter a such that the improper integral Integrate[(x^(a - 1))/(1 + x), {x, 0, ∞}] ...
8
votes
1answer
237 views

Why is this infinite series wrongly computed by Mathematica?

Could you let me know if Mathematica (newer versions) is able to correctly compute this one? Sum[(-1)^(n + 1) Cos[3^n x]^3/3^n, {n, 1, Infinity}]
0
votes
1answer
139 views

Integration of linked PDFs over probability simplex

I'm trying to integrate the following expression: f[a1,p1+p2+p3+p4]*f[a2,p2+p3+p4]*f[a3,p3+p4]*f[a4,p4] Where ...
0
votes
1answer
91 views

Numerical vs Symbolic Integration: Loss of precision

I am trying to integrate the following expression over the time interval $0\leq t \leq \text{period}$. ...
2
votes
4answers
397 views

How to calculate this integral?

I am trying to Integrate the following Integral : $\int_1^{\infty } \dfrac{\left(x^2-1\right)^{13/2} e^{-ax} }{x^{10}} \, dx \,\, \,\,\,\,\,\,\,\,\,\,\,\,(a=\textrm{real>0})$ Mathematica didn't ...
2
votes
2answers
334 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!!!
20
votes
3answers
362 views

Negative probability?

Bug introduced in 9.0.1 and fixed in 10.0.2 I am trying to get the sum of the squares of seven random variables, all uniformly distributed. This is what I tried. ...
2
votes
0answers
88 views

Integration of UnitStep or HeavisideTheta

I am using Mathematica 9. If I evaluate Integrate[Piecewise[{{1 - r, r > 0 && r < R}}, 0], r] then I get ...
0
votes
2answers
177 views

Numeric mixed derivatives

In short I need NumericCalculus`ND extended to mixed partial derivatives. This can be done by nesting them, but care has to be taken to evaluate underling ND only when the parameters become numeric. ...
10
votes
3answers
324 views

Correct way to integrate a certain function

I want to integrate $$ \int_{0}^{2\pi} dt_1 \frac{a + b*\cos(t_1 - t_2)}{c + d*\cos(t_1 - t_2)} $$ where $a, b, c, d, t_2$ are real numbers and $c + d > 0$ & $0 \leq t_2 \leq 2\pi$. I used ...
2
votes
2answers
115 views

Graphical solution of a a two-variable system containing an equation and an inequality

I'm currently studying atmospheric-gravity waves and their dispersion relations. I want to have look at values of $(x,y)$ that solve the following system : \begin{equation}\label{}(1) \left\{ ...
0
votes
0answers
89 views

Complex integrals and residues

I have a question. I try to integrate this function: Integrate[Exp[I k z]*Exp[-Sqrt[k^2 - A^2]x]*Exp[-k^2*L^2/4]/Sqrt[k^2 - A^2], {k, -Infinity, -A}] But ...
0
votes
1answer
68 views

Defining arbitrary order derivatives function [duplicate]

I have a function $\phi(t)$ and I want to define a function $\phi D(t,r):=\phi^{(r)}(t)$, in other words evaluate the $r$'th derivative of $\phi$ at $t$. The naive code that I have is: ...
25
votes
2answers
1k views

Bug in mathematica analytic integration?

I found Mathematica provides me a wrong answer for a relatively simple analytically solvable integral: ...
11
votes
3answers
515 views

Determine frequency of oscillations

I am wondering how I could determine the frequency of oscillations of a differential model equation? How could I find the frequency from this example given in Mathematica Documentation: ...
1
vote
1answer
100 views

Two graphs in the same plot [closed]

I have the following plots for two graphs: ...
2
votes
3answers
116 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 ...
8
votes
2answers
715 views

Is Mathematica really getting this limit wrong?

I'm trying to calculate: $$ \lim_{\beta \to \infty} \tanh \left( \beta A \right) = \mathrm{sgn} \left( A \right) $$ for $A \in \mathbb{R}$, $\mathrm{sgn} \left( \# \right)$ being the sign function. ...
0
votes
0answers
123 views

A power series expansion

Consider the function, $f(z) = z tanh(\pi z) log (z^2 + a^2)$ for some $a>0$. Now I am considering 3 different situations, $z = i(n+0.5) - i\epsilon + \delta - it$ for $n \in \mathbb{Z}$ and ...
1
vote
1answer
100 views

About doing an sum

I want to explain this summation to Mathematica, For $a >0$ define $n_u, n_d \in \mathbb{Z}$ such that if $\{ a \} \leq 0.5$ then $n_u = [a]-1$ and $n_d = - [a]$ and if $\{ a \} > 0.5$ then ...
6
votes
1answer
116 views

Is it possible to circumvent a bug inside SeriesCoefficient?

As far as I can tell, there seems to be a bug in SeriesCoefficient: ...