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)

16
votes
5answers
421 views

Mismatch between numerical and analytic evaluation of an integral

I evaluated the following integral $$\int_0^1 \sqrt{r} \left | \cos \left(\left(k+\frac{1}{2}\right) \pi r\right)\right | dr$$ ...
2
votes
1answer
330 views

Solve ODE $d^2u/dx^2 + u/A = 0$

How can I solve following ODE with Mathematical: $$d^2u/dx^2 + u/A = 0 \quad (\text{or } \ C),$$ with the conditions: $$ \left.\frac{d^2u}{dx^2}\right|_{x=0} = 0, $$ $$u(x=0) = B$$ and ...
4
votes
2answers
272 views

Integrating a periodic function

I have a periodic function ff: ff := Function[x, Piecewise[{{ff[x - 1], x >= 1}, {2 x, 0 <= x < 1}, {ff[x + 1], x < 0}}]] Plotting it works fine: ...
1
vote
2answers
152 views

Integral with unreliable result

I want to calculate $\int_R^1 \sqrt{r} |\cos((k+\frac{1}{2})\pi r)|dr $ and I get a result from Mathematica. Then I try to check the result putting the value of $k$ and $R$, (k=1 and R=0.5) in the ...
3
votes
1answer
119 views

a mysterious error in multiple symbolical integration

I am trying to do the following integration: ...
2
votes
0answers
346 views

Symbolic matrix calculus: What's new in Version 9

I have seen some of the related posts, which ask about doing matrix calculation on tensors unknown dimensions. One of the posts mentioned that version 9 has some capability, but it was concerned about ...
2
votes
1answer
186 views

Check for holomorphy of a function

Given a (rather complicated) function H(z), what is the best approach to check symbolically whether it is holomorphic? What I tried is checking explicitly the ...
8
votes
1answer
240 views

Expanding derivatives of hypergeometric functions

Sometimes Mathematica expresses results of integration or summation in terms of symbolic derivatives of Hypergeometric2F1 function, and cannot further simplify ...
6
votes
4answers
338 views

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

Mathematica can't find a solution to this Expectation. ...
2
votes
1answer
106 views

Finding a root of a parameterized integral

I have a function given as a parameterized definite integral: f[a_] := Integrate[BesselJ[0, x - a] BesselJ[0, x + a], {x, -∞, ∞}] I suspect it has a root near ...
0
votes
1answer
87 views

A problem with “RiccatiSolve::ncsol: Cannot compute a stabilizing solution to the matrix equation ”

I have a SSM ( State Space Model ) and I want to calculate its optimum gains unisng LQOutputRegulatorGains[] command. When I run the code, it solves the equations ...
0
votes
1answer
232 views

Integral over geometric region [duplicate]

I'd like to calculate this integral $$ \int_E y\ dydz $$ where $E = \{ (x,y,z) \in R^3 : z^2+6 < y^2 < 5z \}$ By hand i've got $\frac{1}{12}$ but i'm not sure, and i'd like to verify this ...
2
votes
3answers
794 views

Implicitly differentiate an equation, then solve the resulting equation

Suppose I have an extremely tedious equation to differentiate and want mathematica to help do the differentiation and solve. Consider a less tedious equation: $$y (x,z) = sin \left(\frac{1}{x} ...
0
votes
0answers
108 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
90 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
147 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 ...
0
votes
1answer
113 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
347 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
355 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$ ...
2
votes
1answer
132 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 ...
1
vote
1answer
129 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
199 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
251 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 ...
3
votes
6answers
1k views

Compute Triple Integral on spherical coords

I need to compute: $\int \int \int z dxdydz$ over the domain: $\{x^2+y^2+z^2\leqslant 16,z\geqslant 0\}$ Im trying to use spherical coords as: $$\int_{0}^{2\pi} \int_{0}^{\frac{\pi}{2}} ...
0
votes
1answer
896 views

Second derivative implicit differentiation using Wolfram Alpha input?

How would you perform second derivative implicit differentiation using Wolfram Alpha input? The reason that I'm using WA input is that it gives you step-by-step solutions and I'm a first year ...
1
vote
1answer
110 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
138 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 ...
10
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 ...
1
vote
1answer
106 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
74 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
89 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 = ...
2
votes
3answers
411 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
354 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
605 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
1answer
107 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
89 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
57 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] ...
1
vote
1answer
697 views

Solving a complex-valued differential equation with NDSolve

I am trying to solve $dx/dt=\sqrt{1+(ix)^{1.8}}$ for initial condition $x[0] =-0.9877 + i 0.1563$, where $x$ is a complex variable. I would like to plot the imaginary part of the solution versus the ...
5
votes
2answers
272 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
672 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, ∞}] ...
12
votes
1answer
348 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
247 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
99 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
464 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
430 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
4answers
425 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
136 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
261 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
4answers
368 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 ...