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
2answers
60 views

Taking derivatives

This is my code: X := g[P] P = a*t R = t*X FOC = Dt[R, t] This is what I get for the last line: ...
1
vote
1answer
229 views

How to differentiate an integral on Mathematica?

I am so new on Mathematica. I try to find the first variation of this function according to $t$ on mathematica but I could not achieve. Here is the function ; ...
10
votes
4answers
443 views

Mathematica complaints that convergent integral diverges

Bug introduced in 10.0 and fixed in 10.0.2 Trying to do the integral ...
2
votes
0answers
193 views

Mathematica not evaluating q derivative of Jacobi theta function

Jacobi theta functions, $\theta_a(u,q)$ for $a=1,2,3,4$ are defined in the unit disk $|q|<1$. For some reason that I would like to understand, Mathematica does not evaluate numerically the $q$ ...
3
votes
3answers
202 views

Factoring a scalar function out of a vector

I've taken a first and second derivative: ...
2
votes
1answer
241 views

Definite Integral over Bessel Function

Hello I am interested in evaluating the following integral. ...
4
votes
1answer
184 views

Why can't the Bernstein Bears work normally?

Bug introduced in 7.0 and fixed in 9.0 I want to use the built-in BernsteinBasis[] to learn about Bezier curves. I tried the following code: ...
10
votes
2answers
262 views

Mathematica integration failure - new or old?

Bug introduced in 7.0.1 or earlier, persists through 10.1 Consider the following integral: Integrate[Log[a Cos[x]^2 + b Sin[x]^2], {x, 0, 2Pi}] This takes a ...
2
votes
2answers
523 views

Integration of piecewise function

I'm trying to understand why Mathematica is not evaluating a piecewise function, while it's able to evaluate each of the regions separately. This fails: ...
8
votes
3answers
651 views

Biot-Savart Law/Magnetostatics Solution

I am working on solving the Biot-Savart Law equation for the magnetic field around a charged ring of uniform current density. The expression that Mathematica gives is rather nasty, as expected. Notice ...
7
votes
5answers
2k views

Find function inverse

I'm trying to find the inverse of a function: (30*x^2 (1 - x)^2) (* where 0<x<1 *) I tried all the following options: 1. ...
3
votes
1answer
1k views

How to solve time dependent Optical Bloch Equations for a three level system?

I am trying to solve the optical block equations for three level system which are time dependent, which i generated using the "Atomic Density Matrix" package. When i use the command to solve this 9 ...
1
vote
1answer
63 views

NMaximize not substutiting parameters of objective function [duplicate]

NMaximize is not substituting parameters of its objective function. Is there a way to fix that? I have a three line program that exhibits the problem. ...
10
votes
2answers
516 views

Interpolating an Antiderivative

I'd like to be able to make InterpolatingFunctions for antiderivatives of functions that can't be integrated symbolically. However, the following code returns ...
1
vote
2answers
64 views

Convergence problem when I try to integrate a function of one variable and one symbolic parameter

I have an integral with a parameter a and a > 0 ...
2
votes
1answer
197 views

Double integral changing the limits

This integral does not evaluate: Integrate[ Integrate[1/(E^((a - c)^2 + (b - c)^2) ), {c, 0, 1}], {b, -Infinity, Infinity}] However, reversing the integrals ...
0
votes
3answers
644 views

Is there a way to rewrite integrals in Mathematica using u substitution?

In Mathematica, can I give it an integral and a few substitution rules and have it rewrite the integral in terms of those variables?
10
votes
3answers
462 views

Maximize violating constraints

I have Maximize[{(h*10)/(300*(100 - (l^.5 + d^.4 + H^.6))), (l + d + H + h) == 669, l > 0, d > 0, H > 0, h > 0}, {h, l, d, H}] I believe ...
1
vote
2answers
277 views

Integrating Squared of Interpolating Function with respect to one variable

I am interested in evaluating a two dimensional interpolating function produced by solving the wave equation. Here is the code that includes the resulting interpolating function. ...
4
votes
0answers
243 views

Strange behaviour of MMA in derivatives of some standard functions

There are some peculiar things to be discovered in derivatives of some standard functions in MMA: Strange behaviour Example 1: Abs We have ...
6
votes
2answers
418 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 ...
0
votes
0answers
55 views

Custom simplification of symbolic definite integral with regards to bounds

I would like to rewrite symbolic definite integrals such that the bounds are always positive. For example : $\int_{-\infty}^\infty f(x) \mathrm{d} x \rightarrow \int_0^\infty \left[ f(x) + f(-x) ...
5
votes
1answer
258 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$ , ...
2
votes
2answers
198 views

How to calculate the analytical result of this singular integral? [duplicate]

I'm trying to solve this integral: $I(a,b)=\displaystyle{\int_1^{\infty }\dfrac{e^{-ax}}{\sqrt{x^2-1} (1+x\sqrt{1-b^2}) (x-1/\sqrt{1-b^2})} \, dx} \,\, \,\,(a \;\textrm{real>0} \, , 0<b<1)$ ...
3
votes
3answers
185 views

When NSolve fails due to a differential situation?

c = 1.1111; y[x_] = x - c Sin[x] NSolve[y[x] == 0, x] The method has procured no result. Successive derivatives were plotted in an attempt to fix the problem. ...
26
votes
3answers
873 views

Symbolic integration error

fixed in 10.1 (windows) I'm running Mathematica 10.0.0 and encountered a disturbing error in the symbolic integration of a rather simple function ...
4
votes
1answer
215 views

Efficient Dyson series implementation

I'm trying to implement a Dyson series \begin{array}{lcl} U(x,x_0) & = & 1 + \int_{x_0}^{x}{dy_1V(y_1)}+\int_{x_0}^x{dy_1\int_{x_0}^{y_1}{dy_2V(y_1)V(y_2)}}+\cdots \\ & &{} + ...
1
vote
0answers
117 views

Giving hints to Integrate

I working with the integral ...
1
vote
1answer
62 views

Specifying independence of a variable in a function

I am dealing with a function of a finite number of variables and among the operations I wish to do is to differentiate it repeatedly with respect with any of the variables. Now this function happens ...
35
votes
6answers
2k views

Finding length of intersection of two surfaces

I would like to know how we find the length of the intersection of two surfaces. For instance, in the following example,a surface intersects with a plane: How do we find the length of intersection ...
5
votes
1answer
242 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 ...
1
vote
0answers
367 views

Differential equations with a complex variable

Is Mathematica able to handle ordinary differential equations where the variable itself is complex? I am looking for solutions of ODE systems of the form $$\left\{ \begin{align} i\frac{da_1}{dt} ...
6
votes
1answer
264 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]? ...
0
votes
1answer
236 views

Sum rule in integration not working [duplicate]

I am new to Mathematica and I just tried to play around with some integrals. But the sum rule in integration doesn't work. Here a simple example: Mathematica doesn't recognize that $\int_0^h (1-f(x)) ...
4
votes
2answers
348 views

Integral returns function of variable that was integrated over

Bug introduced in 8 and fixed in 10 When I enter ...
0
votes
0answers
39 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
160 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
696 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
393 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
295 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
79 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
229 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
111 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
172 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
400 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
130 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
514 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
572 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
396 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: ...