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)

4
votes
2answers
326 views

Integral returns function of variable that was integrated over

Bug introduced in 8 and fixed in 10 When I enter ...
0
votes
0answers
36 views

different values of integral with/without bogus assumptions

I do not understand why I get different answers for the following two identical integrals. ...
1
vote
1answer
152 views

Under what conditions does Mathematica give a natural antiderivative when integrating?

Natural antiderivative is defined as follows, using Fourier transform: $$f^{(-1)}(x)=\frac{i}{2\pi}\int_{-\infty}^{+\infty} \frac{e^{- i \omega x}}{\omega} \int_{-\infty}^{+\infty}f(t)e^{i\omega t}dt ...
1
vote
2answers
685 views

Computing 10-dimensional volume of a 9-sphere [closed]

I'm trying to compute 10-dimensional volume of a 9-sphere with radius r using Monte Carlo. ...
2
votes
3answers
382 views

How can I apply calculus to functions obtained from NDSolve?

Originally, I asked the question below, but the real underlying issue is as follows: When we solve an ODE numerically, I get the answer like this: ...
3
votes
2answers
293 views

Integral converging in M9 but not in M10

Background: Mac OSX 10.9.4, Mathematica 9.0.1 vs. Mathematica 10.0 (both versions are Student Editions). I have a notebook. It used to evaluate fine in Mathematica 9. I upgraded to 10 and, without ...
2
votes
0answers
98 views

Result of symbolic integration changed drastically by making assumptions

I would like to know the underlying reason for different outcomes for the two integration operation below. One of them includes a few assumptions, otherwise both have the same integrand: ...
1
vote
0answers
73 views

Getting long complex-valued integrals when simpler real-valued expressions exist

I have a long list of real-valued functions I'd like to integrate symbolically. For many of them, Mathematica gives me results with long complex-valued expressions involving weird functions such as ...
1
vote
1answer
212 views

Is it necessary to introduce ImplicitRegion and ParametricRegion in version 10?

In version 10.0, Mathematica introduced a new function about region plotting: ImplicitRegion. But I'm wondered that haven't this function been already existed in the older version? That is, ...
5
votes
0answers
110 views

Convoluting inverse square root with Gaussian

I would like to convolute the inverse square root on the interval [0,inf] with a Gaussian function, like so: ...
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 ...
7
votes
3answers
373 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 ...
2
votes
2answers
128 views

Integrate returns unexpected result

Consider the following function $$g(x,y):= \frac{1}{( (1+y)^2+x^2 )( 1+ax^2y^2 )^2}$$, where I assume that $y\geq 0$ and $a\in (0,1]$ is a parameter. When I try to evaluate the integral $\int ...
17
votes
5answers
481 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
483 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
346 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
161 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
123 views

a mysterious error in multiple symbolical integration

I am trying to do the following integration: ...
3
votes
0answers
491 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
207 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 ...
9
votes
1answer
293 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
383 views

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

Mathematica can't find a solution to this Expectation. ...
2
votes
1answer
117 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
95 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
314 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
1answer
80 views

Limits explanation in doc

Limit[expr,x->Subscript[x, 0],Direction->1] computes the limit as x approaches Subscript[x, 0] from smaller values. ...
2
votes
3answers
1k 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
112 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
97 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
165 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
89 views

How can I avoid a scoping problem when differentiating?

I simplify my real problem as follows: I define a rule function as Clear[rule]; rule[i_] := tt -> i; the real rule function is not this simple. I define a ...
0
votes
1answer
140 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
407 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
396 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
146 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
139 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
224 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
258 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
2k 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
1k 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
121 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
151 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 ...
13
votes
9answers
2k 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
106 views

Why does this integral have a complex result

I wanted to find the normalization of distribution F2 . ...
1
vote
1answer
116 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: ...
-3
votes
1answer
120 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
499 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
467 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
866 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 ...