Tagged Questions

Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

14
votes
1answer
357 views

How to represent a continuous monotonic phase of Airy functions?

Note: In this question I am concerned only with real-valued variables and functions. DLMF, §9.8 Airy Functions, Modulus and Phase, formula $9.8.4$ defines the phase of Airy functions: ...
7
votes
2answers
296 views

Find asymptotics of Sum[2^i*Binomial[n-i-1,2*n/3-1],{i,0,n/3}]

I have an expression 2^n / Sum[ 2^i Binomial[ n - i - 1, 2n/3 - 1], { i, 0, n/3}] ...
1
vote
0answers
137 views

How to prevent simplification of hypergeometric functions resulting from integrations?

Definite integrals from 0 to Infinity over a product of two hypergeometric (including exponential, trigonometric, hyperbolic, ...
3
votes
2answers
275 views

Integrating special functions

I would like to integrate the following function with Legendre polynomial and Gamma function. I am open to suggestions. ...
3
votes
1answer
177 views

Symbolic Integration of Special Functions

Sorry in advance if this formatting comes out strange, this is my first question! I can't find a way to integrate, e.g., a function of the Hermite polynomials for general (still integer) order. For ...
3
votes
0answers
54 views

SiegelTheta fails to evaluate when given proper arguments

SiegelTheta often returns error messages when I give it arguments that should be of the correct form. For instance, I have a numerical matrix like ...
5
votes
2answers
297 views

Numerical Integration with InverseErfc

I am trying to numerically integrate an equation that involves InverseErfc (embedded in the copula defined). The equation looks like the following: $$ \int_0^T \int_0^\infty \int_0^\infty ...
2
votes
0answers
151 views

Real integral evaluating as indeterminate

Mathematica evaluates the following integral as: ...
5
votes
1answer
209 views

Accurately evaluating the hypergeometric function

As part of another problem, I am working to evaluate hypergeometric functions such as Hypergeometric2F1[1, 1, n, -1] for large $n$. I am hoping to obtain at ...
7
votes
1answer
378 views

Numerical contour integrations in the complex plane - contour deformation gives different answer for analytic kernel

I am trying to do a contour integration in Mathematica numerically. In particular, I'm checking the identity: $$ H_m^{(1)}(z) =\frac{i^{-m}}{\pi}\int_{-\pi/2 + i \infty}^{\pi/2 - i \infty} \exp[i m ...
5
votes
1answer
148 views

How can I program the RiemannR function using the LogIntegral command?

I would like to program the RiemannR function using the LogIntegral command because I would like to later experiment with a ...
0
votes
1answer
361 views

Searching for roots of complex function

I'm searching for roots of complex function $$ 2\imath q \ln(-2\imath k)+\imath\pi-2\imath \Im(\ln(\Gamma(1+2\imath q)))+\ln(\frac{\Gamma(1+\imath q-\imath q x/k)}{\Gamma(1-\imath q-\imath q ...
13
votes
2answers
352 views

What kind of hypergeometric function is it?

I found a formula for an integral of a product of three Bessel functions at The Wolfram Functions Site: I cannot understand what kind of hypergeometric function it is. The Mathematica code given ...
7
votes
1answer
484 views

Visualizing vector-spherical waves

This is a follow-up question to this one on visualizing vector-spherical harmonics. This time, I would like to visualize the vector spherical waves (including the radial dependence). The functions ...
2
votes
3answers
200 views

Why do these two different zetas produce the same value?

Zeta[-13] == Zeta[-1] == -1/12 Why do these two different zetas produce the same value?
8
votes
2answers
199 views

How to simplify an expression with special functions to zero

The following is a well-known Bessel function identity: $$J_{-n}(z)=(-1)^n J_n(z),\qquad n\in\mathbb Z$$ To check this, I used the following code and the result is as what I expected. ...
15
votes
2answers
1k views

Visualizing vector spherical harmonics

I have painstakingly derived the vector-spherical harmonics $\mathbf{V}_{J,\,M}^\ell(\theta, \phi)$, which are the generalization of ordinary spherical harmonics $Y_\ell^m(\theta, \phi)$ to vector ...
9
votes
1answer
484 views

Why can't I change the value of MaxRecursion in NIntegrate when integrating BesselJ?

I am trying to evaluate this integral numerically $$ \int_0^{\infty } J_0(q R) \tanh(q) \, \mathrm{d}q $$ for large values of $R$. This makes the integrand oscillate more quickly and Mathematica ...
9
votes
1answer
277 views

How to calculate this integral? Integrate[BesselJ[0, x - BesselJZero[0, 1]]/x, {x, -Infinity, Infinity}]

I tried to calculate the following integral, but it returned unevaluated. ...
2
votes
0answers
97 views

RSolve and incomplete gamma function

I am interested in how Mathematica uses the incomplete gamma function to solve difference equations. For example, if we have this inhomogeneous, 2nd order equation and use ...
7
votes
2answers
302 views

Problem with NIntegrate when WorkingPrecision is specified

I am trying to evaluate this integral numerically: $$ \int_0^{\infty } m \exp (-m) J_1(m){}^2 \, dm $$ Everything is OK when only the integration method is specified: ...
3
votes
3answers
281 views

Mathematica not able to confirm its own solution to differential equation

I type the following into Mathematica: DSolve[q''[x] + 2 x/(x^2 - 1) q'[x] - 4*q[x]/(x^2 - 1) == 0, q[x], x] It gives me the result ...
17
votes
1answer
696 views

Peirce's quincuncial projection

The Peirce quincuncial projection is the cartographic projection of a sphere onto a square. In short, I would like to see it implemented in Mathematica. Here is my code: ...
0
votes
1answer
99 views

Reflection transform of function [duplicate]

I am trying to find the reflection function. Here is my function and its graph. ...
1
vote
2answers
1k views

How to take conjugate of a function?

Naïvely this is what happens and it obviously is not helpful! ...
3
votes
1answer
117 views

Find point at which equation stops having roots (if it exists)

I am interested in the roots of this function: f[M_, b_] := 1 - (2 M Gamma[2, 0, (1/M + b M)/Sqrt[b]])/(1/M + b M) for fixed values of b. In particular I want ...
5
votes
2answers
176 views

Why does Integrate return a solution that is not defined at a particular point when it actually is well defined at that point?

I am trying to compute Integrate[Sqrt[x^4 + (y - y^2)^2], {x, 0, y}] Mathematica 8 gives ...
5
votes
5answers
1k views

Series expansion in terms of Hermite polynomials

I am trying to expand a polynomial in terms of orthogonal polynomials (in my case, Hermite). Maple has a nice built-in function for this, ChangeBasis. Is there a ...
5
votes
2answers
294 views

Mathematica cannot calculate a limit

When I evaluate Limit[E^(-n)*Sum[n^k/(k!),{k,0,n}], n -> ∞] Mathematica gives me the result ...
1
vote
0answers
120 views

Getting poles of a Gamma functions [duplicate]

Why do the following 2 sequences give different answers? n = 1.5 Series[Gamma[0.5 - n - x], {x, 0, 2}] Series[Gamma[-1 - x], {x, 0, 2}] (..clearly the output from the second expression is ...
4
votes
1answer
230 views

Simplifying numerical expressions involving special functions

I've encountered the following problem. There is the identity that the special functions EllipticK[x] and EllipticE[x] satisfy: ...
3
votes
1answer
205 views

Strange behaviour of PolyLog Function

I discovered some strange behaviour of the PolyLog[] Function in Mathematica which seems to me like a bug in the function implementation. It looks like ...
1
vote
2answers
609 views

Irregular Confluent Hypergeometric Functions (Spherical Coulomb Wavefunctions)

I want to program in the regular and irregular spherical Coulomb wavefunctions $F_\ell(\gamma,kr)$ and $G_\ell(\gamma,kr)$, respectively, which are defined in terms of the regular and irregular ...
5
votes
0answers
103 views

Expansion of $E(i c \mid m)$ at $c\to\infty$?

Currently, I am using a Windows machine with Mathematica 8. I noticed a difference in a series expansion of the function EllipticE[] in comparison with a result ...
0
votes
2answers
491 views

Using Mathematica to find poles of Gamma functions

I am concerned about the expression on the RHS of equation A.5 (page 19) in this paper: $$\int\frac{d^d ...
6
votes
2answers
267 views

Cannot Get Numerical Results to Match

I try this numerical summation (in two parts) ...
6
votes
0answers
299 views

Reproducing the Integral Definition of the Modified Bessel function

I need to simplify some integral expressions in terms of special functions, such as the modified Bessel function of the first kind. See for example Eq. (5) on this page. Notice that the real ...
2
votes
1answer
258 views

Directional derivative of SiegelTheta

I'm working on a problem where I have to integrate both the Mathematica function SiegelTheta and some of its second order directional derivatives. Using the function works well but something goes ...
1
vote
1answer
327 views

Plotting Fresnel function

I am trying to plot the partial sums and the cesaro means of the function $\sqrt{|x|}$ and for $a_{n}$, I obtained the following code which contains FresnelS. ...
9
votes
2answers
2k views

problem with coloring spherical harmonics

I want to color a spherical harmonics. So I write as follows. ...
8
votes
1answer
1k views

solve an integral equation numerically

I am trying to find a numerical solution for an equation of the form: $$ f(t)=\int_{t_{min}}^{t} {\exp[2 (t^\prime - t)] E(t^\prime) f(t^\prime)} + \int_{t}^{0}{ \exp(t^\prime - t) E(t^\prime) ...
9
votes
3answers
770 views

Multi-Factorial and Series with Triple-factorial terms

Let $n!^{(k)}$ denote a multi-factorial which is defined by $$ n!^{(k)} = \begin{cases} 1 & n \leqslant 0, \\ n, & 0 < n < ...
3
votes
0answers
148 views

Why doesn't Log[Gamma[]] simplify to LogGamma[] where it could?

I have been playing with various equations involving amount of permutations in relatively large sets. Easiest way to look at these is something like Log[10, bignumber!] . Often expressions, even ...
8
votes
1answer
905 views

Implementation of Incomplete Fermi-Dirac Integral in Mathematica

I'm working on a special algorithm to implement a more accurate effective mass calculation for hole carriers in silicon in Mathematica. This rather involved algorithm uses incomplete Fermi-Dirac ...
4
votes
2answers
166 views

LevinRule and SphericalBessels

I'm currently looking at a simplified problem that approximates another problem I'm looking into. In this simplified problem I at least have an analytic integrand and can easily provide all info on ...
10
votes
1answer
604 views

Is my expression too complicated for FullSimplify or am I doing something wrong?

I have a messily defined function $v(h, w)$ with $h, w \in \mathbb{R}$ and with a removable singularity at $h=1/2$, and I am trying to prove some of its properties using Mathematica. In particular I ...
5
votes
2answers
242 views

Solve for $a$ as a function of $\beta$?

I am trying to solve this equation: $$\beta^{-a} \Gamma(a) \sin(a \pi) + e^\beta \beta^{2 a - 1} \Gamma(1 - a) \sin(a \pi) = 0$$ I tried the following: ...
18
votes
0answers
653 views

Fast Spherical Harmonics radiative transfer

This is a rather specific question and I apologize for spamming you with some lengthy code. But it could be interesting for some reader and maybe you can help out, so please bear with me. I am using ...
0
votes
1answer
287 views

Iterative way to find roots of confluent hypergeometric function

I am trying to find roots of confluent hypergeometric function and I wonder if I can choose the initial guess by the choice of $\beta$. ...