Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

1
vote
3answers
229 views

Find all roots of a function with parabolic cylinder functions in a range of the variable

I want to find all roots of a function involving Parabolic Cylinder Functions. In what follows, I define 2 variables $\xi1$ and $\xi2$, which in turn depend on $\omega$. My function is then defined as ...
3
votes
1answer
69 views

Interval arithmetic for DawsonF (or other special functions)

I am currently trying to estimate a complicated expression involving DawsonF using interval arithmetic. The interval arithmetic is partially supported by ...
5
votes
0answers
61 views

Huge difference after changing a fraction to decimal

I have a limit to calculate. With[{p = 1/2, q = 0.1}, Limit[a^q Integrate[Sin[x]/x^p, {x, a, Infinity}], a -> Infinity]] gives correct result ...
3
votes
2answers
285 views

Roots of Whittaker W function

I am interested in finding the roots $u$ of the equation $$ W_{1,\imath b}(a)=0, $$ where $W_{\kappa,\mu}(z)$ denotes the Whittaker $W$ function, $a>0$ is a fixed parameter, $\imath=\sqrt{-1}$ and ...
5
votes
0answers
94 views

SiegelTheta throws errors from calling Range with complex arguments

Bug introduced in 10.4 or earlier. This may or may not be related to the bug reported in this question. I was trying to verify the results of this challenge over on Code Golf with the following ...
2
votes
0answers
59 views

Ramanujan's asymptotic formula for $p(n)$

The following expression is the exact formula of partition function. Ramanujan's asymptotic formula for $p(n)$ is following $$p(n)\sim\frac{1}{4n \sqrt{3}}e^{\pi \sqrt{2n/3}}$$ Can I use ...
1
vote
1answer
81 views

NIntegrate gives two results for two forms of the same function. Which one to trust?

I am interested in evaluation the following integral numerically (since apparently there is no analytical solution) $$\int dx \,x^3 \left(e^{2 i c x }-i \text{erfi}\left(\frac{x +i c ...
5
votes
1answer
70 views

ComplexInfinity for a convergent product

The infinite product involving the ratio of (n^2)! to its Stirling approximation ...
10
votes
1answer
134 views

Teaching Mathematica more about DiracDelta and KroneckerDelta

As the documentation and some experimentation indicates, Mathematica contains little information about representations of the DiracDelta and ...
3
votes
0answers
74 views

Bug in Series[Pochhammer[1+n, n], {n,Infinity,1}]?

The function Pochhammer[1 + n, n] tends to infinity. We have ...
2
votes
0answers
43 views

Unexpected imaginary term in the asymptotic expansion of DawsonF

FullSimplify[Series[2 DawsonF[x], {x, ∞, 8}]] (* -I Exp[-x^2] √π + (1/x + 1/(2 x^3) + 3/(4 x^5) + 15/(8 x^7) + O[1/x]^9) *) What is the reason the term ...
2
votes
1answer
90 views

Mathematica doesn't know about the absorption identity? [closed]

The well-known binomial absorption identity states that $$n\binom{n-1}{k-1} = k \binom{n}{k}$$ However, Mathematica gives ...
0
votes
1answer
95 views

Elliptic integral of the first kind [closed]

I want to plot the integral $$I(\phi) = \int_0^{\phi} \frac{\mathrm{d} \theta}{\sqrt{1 +\sin(\theta)^2}}$$ In Mathematica notation, it is a case of an elliptic integral of the first kind with $m=-1$, ...
3
votes
1answer
24 views

Asymptotics of Bessel function for real arguments

I am trying to calculate the following asymptotic behaviour: Normal@Series[BesselK[1, r \[CapitalLambda]] / BesselK[1, \[CapitalLambda]], {r, 0, 1}] but for ...
1
vote
2answers
95 views

How to convolve the unit box function and the modified Bessel function of the second kind in 2D?

In 1D the convolution of the unit box function and the modified Bessel function of the second kind $K_0(x)$ works very well. ...
12
votes
2answers
813 views

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

Bug introduced in 8.0.4 or earlier and persists through 10.4. 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$. ...
8
votes
1answer
294 views

Show how Mathematica defines a function

Is there a way to show how Mathematica defines a function, such as In: Something[Sqrt], Out: Sqrt[x] -> x^(1/2) As far as I understand it the command ...
1
vote
2answers
55 views

Finding roots of a function that includes Bessel functions [duplicate]

I'm fairly new to Mathematica so forgive any stupid mistakes. Here's my function: ...
5
votes
2answers
174 views

Definite integral closed-form expression

Is there a way to get Mathematica yield a closed-form expression (in terms of special functions) for the integral: $$ \int_{0}^{\infty} e^{-a t}\log(t)\log(1+t)\,dt, $$ where $a>0$? The obvious ...
1
vote
0answers
54 views

Output in FunctionExpand for function of Gamma

I used of this code: α Gamma[α] // FunctionExpand and get output: Gamma[1 + α] Also I used of this code: ...
1
vote
1answer
69 views

Macdonald-Koornwinder polynomials?

Does Mathematica have an internal implementation of the Macdonald-Koornwinder polynomials? (Also called Koornwinder polynomials.) I looked online but could not find it.
1
vote
2answers
90 views
0
votes
1answer
207 views

Solution of differential equation in terms of incomplete gamma function

I need help in solving equation 15 and 16 either manually or in Mathematica to get the solution in terms of the incomplete gamma function. This is what Mathematica tells me. I can't understand ...
0
votes
0answers
58 views

Another question about inverse Laplace transform

Given $r = 0.06;\quad \theta = 105;\quad \kappa = 1;\quad x_0 = 100;\quad K = 100;\quad \sigma = 0.10;\quad T = 0.25;$ Define $ \nu = -\kappa/\sigma^2 - 0.5;\quad p = \kappa*\theta/\sigma;\quad q = ...
8
votes
2answers
806 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 ...
1
vote
0answers
34 views

Derivative wrt to order of MacDonald function

I'm trying to get the following result confirmed in Mathematica: $$ \left.\frac{\partial\mathop{K_{\nu}}\nolimits\!\left(x\right)}{\partial\nu}% ...
4
votes
0answers
73 views

SiegelTheta gives misleading message when the dimensions don't match

Bug introduced in 6.0 and persisting through 10.4 SiegelTheta is new in 6.0 In order to test the SiegelTheta function, I ...
0
votes
0answers
36 views

Find Roots of expression with modified Bessel functions [duplicate]

I would like to solve an equation of this type $$A-B\,K_0(\kappa\,r)-C\,K_1(\kappa\,r)=0\ \ \quad \text{with}\ A,\,B,\,C,\,\kappa\in\mathbb{R},\ \text{known}$$ for $r$. I am not aware whether this ...
2
votes
1answer
48 views

Getting rid of sub-exponential terms in an asymptotic expansion for a modified Bessel function

I am trying to get Mathematica to produce suitable asymptotic expansions for some modified Bessel functions at large argument (more specifically, the expansion in the DLMF's eq. (10.40.1)), and I'm ...
3
votes
2answers
85 views

DifferenceRoot question

I was doing the following sum: $$\sum_{i=2}^k \frac{(-1)^i}{i-1} \binom{2k-i-1}{k-1}x^i$$ First, Mathematica simplifies it to some DifferenceRoot function: ...
5
votes
2answers
308 views

Implementation of a complex recurrence relation for polynomials

I would like to implement the recurrence relation for the polynomials $U_n(x)$ appearing in the large order asymptotics of the Bessel functions. The recurrence in question is: ...
5
votes
1answer
265 views

Find solution of nonlinear ODE in terms of JacobiCN

I am trying to find a specific solution for this differential equation: $-\frac{1}{2}\frac{d^2}{dx^2}\psi(x)-2k \; \psi(x)^3 + \frac{1}{2}k^2\; \psi(x)=0$ MMA gives me a solution in the form of a ...
4
votes
0answers
44 views

Integral Form of Modified Bessel Function of the Second Kind

Why can't Mathematica integrate r = Integrate[Exp[-x Cosh[t]], {t, 0, Infinity}]; r = Assuming[Element[x, Reals], Simplify[r]]; Together[r] From Wikipedia, it ...
8
votes
2answers
222 views

How can I get the solutions of $x^n-x-t=0$ in hypergeometric form?

$x^3-x-t=0$ has three roots that can be expressed in hypergeometric form, for example $$x_1=-1-\frac{t}{2}{_2F_1\left(\frac{1}{3},\frac{2}{3};\frac{3}{2};\frac{27t^2}{4} \right)}+\frac{3t^2}{8} ...
0
votes
0answers
61 views

Riemann Zeta function definition was expanded by Euler with an infinite product series [duplicate]

The Euler infinite product series definition for Riemann's zeta function requires that Mathematica use all prime numbers in the product series. Can anyone help me with the code that will give a ...
3
votes
2answers
168 views

Use Meijer-G function to represent elementary functions

I want to represent these elementary functions: $x^{2}\sqrt{x}$, $\sin{4x}$, and $x\ln{x}$ as cases of MeijerG. What arguments should I give to ...
0
votes
0answers
77 views

How to evaluate this integral

I am trying to evaluate the following integral in Mathematica: ...
8
votes
2answers
242 views

Can your Mathematica do the following integral?

Backslide introduced in v10 and persisting through v10.3.1. I was trying to do the following integral in Mathematica: ...
4
votes
1answer
81 views

Simplify Class invariant $G(25)$

How to simplify $$\frac{\sqrt[3]{\vartheta _3\left(0,e^{-5 \pi }\right)}}{\sqrt[12]{2} \sqrt[6]{\vartheta _2\left(0,e^{-5 \pi }\right) \vartheta _4\left(0,e^{-5 \pi }\right)}}$$ This is a ...
20
votes
4answers
588 views

Expansion of a hypergeometric function takes ages with Mathematica 9 and 10 (regression?)

Mathematica 8 (Linux version) can evaluate AbsoluteTiming[Series[Hypergeometric2F1[1, 1 - eps/2, 3 - eps, 1/2], {eps, 0, 0}]] in no time. On one of the ...
3
votes
1answer
91 views

Is the real spherical harmonic (l = 1, m = 0) really 'bigger' than (l = 1, m = 1)?

Using SphericalPlot3D to plot the real spherical harmonic with l = 1 and m = 0: ...
5
votes
1answer
101 views

Somewhat Irreproducible Integrate Results

Backslide introduced in v10 and persisting through v10.3.1. In the course of considering question 102922, I encountered erratic results from a particular integration. It is illustrated as follows. ...
34
votes
0answers
1k 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 ...
8
votes
2answers
288 views

On the definition of the associated Legendre polynomials

Mathematica computes for n = 1,2,...: (-1)^n (LegendreP[n, -1, -3]/Sqrt[2]) -I, -3 I, -11 I, -45 I, -197 I, ... Maple ...
2
votes
2answers
135 views

Why is LegendreQ[1/2,x] complex-valued for x>1?

Something is strange with $\sf LegendreQ$. Let $x>1$. I wonder why $\sf LegendreQ[\frac12,x]$ is complex-valued, and the following two codes do not give the same results: ...
1
vote
0answers
88 views

Product of two Meijer's Function

I want to evaluate an integral $I_1$ defined in $Eq.(1)$ as \begin{align} I_1=\int_{0}^{\infty}\frac{x\exp(-\beta x)K_1(\alpha x)}{1+x}dx\tag{1} \end{align} Where $\alpha\geq0$, $\beta\geq0$, and ...
0
votes
1answer
59 views

Plotting a function based on complicated integral

I have this function : f[Lambda_] := K Integrate[x Exp[- x^2-Lambda x] HypergeometricU[-Lambda,1/2,(x+ Lambda/2+2)^2],{x,0,Infinity}]; where ...
9
votes
3answers
276 views

Labeling solutions of an Eigenvalue equation involving Bessel functions

I'm solving the Schrödinger equation for a particle in an annular geometry with hard wall boundary conditions and I've reduced it to the following equation: $$J_m(k\,R_1)\,Y_m(k\,R_2) - ...