Questions tagged [special-functions]
Questions on the special mathematical functions implemented in Mathematica.
1,537 questions
3
votes
1
answer
307
views
How to express Meijer G-function into a form involving well-known functions?
I'm working on MAT 14.3, I am attempting to evaluate the integral int[n] for n=0,1,2,3 in symbolic form :
...
0
votes
0
answers
61
views
ComplexExpand on BesselJ
I know ComplexExpand assumes its variables are real by default. However, when applied to BesselJ function, it is different. For ...
3
votes
2
answers
204
views
Converting hypergeometric function to Struve form
In my question on MathOverflow, I was looking for a closed form result of the following sum:
$$\sum _{k=0}^n \frac{(-1)^{n-k} x^{2 k} (2 (n-k)-1)\text{!!}}{(2 k)\text{!!}}.$$
Someone suggested me to ...
1
vote
1
answer
149
views
How do we find an asymptotic approximation of the function $\mathbf{P}(r)=\left(t_1+\prod\limits_{k=1}^{r}(t_2+k^{s})\right)^n$?
Suppose we have the following function, where $s\in\mathbb{R}$ and $t_1,t_2,n\in\mathbb{N}\cup\{0\}$ are constants:
$$\mathbf{P}(r)=\left(t_1+\prod_{k=1}^{r}(t_2+k^{s})\right)^n$$
Question: What is ...
1
vote
1
answer
104
views
Convergence question about HypergeometricPFQ
Update in response to some questions and comments:
The book I mentioned is
Special Functions (Encyclopedia of Mathematics and its Applications, Series Number 71) by George Andrews , Richard Askey ...
0
votes
0
answers
28
views
2D $q$-Racah kernel $\tilde{K}$, the pre-limit barcode kernel $K_{L}^{\text{barcode}}$, and the conjectual limit $\mathcal{K}^{\text{barcode}}(s,t)$
Looking for Mathematica implementations of 2D $q$-Racah kernel $\tilde{K}$, the pre-limit barcode kernel $K_{L}^{\text{barcode}}$, and the conjectual limit $\mathcal{K}^{\text{barcode}}(s,t)$.
Does ...
1
vote
0
answers
62
views
Behaviour of Integrate with GenerateConditions flag
I am trying the following integral with Mathematica 13.3
...
3
votes
1
answer
229
views
Series expansion involving incomplete Beta function
I am trying to expand the following around u=1/2.
...
1
vote
1
answer
110
views
Help validating the solution to a nonlinear ODE: Does it solves other non-linear ODEs too? [closed]
I am working right now in this another question by doing examples and I got stuck in proving the results I found through Wolfram-Alpha are right.
The main differential equation is the following ODE:
$$...
3
votes
2
answers
167
views
Expanding the HypergeometricPFQ function gives complex values
I have a function
HypergeometricPFQ[{}, {2, 2}, -Log[n]]
which for real n > 0 gives a real values, as can be seen on the ...
0
votes
0
answers
73
views
How to improve the Mathematica implementations of q-Bessel polynomial?
I know Mathematica has implemented some QFunctions.
I have found this literature: q-Bessel polynomials, q-Laguerre polynomials, q-Charlier polynomials, etc.?
I have implemented q-Bessel polynomials in ...
1
vote
1
answer
117
views
Using `JacobiP` to approximate `BetaDistribution[a,b]`
I'm looking at generalized moment expansions of BetaDistribution. Below is an example using ChebychevT, I'm looking for help extending this to cover the general ...
1
vote
0
answers
42
views
0
votes
0
answers
122
views
Integrating Hypergeometric2F1 with large arguments
I am trying to integrate a function which contains two hypergeometric functions ${}_2F_{1}(a,b,c,z)$. In my particular case,
\begin{align}
a &= 1 + \ell - n,\\
b &= 1 + \ell + i \lambda x,\\
...
3
votes
1
answer
268
views
Demonstrating integral result involving Bessel function
On my laptop Mathematica does not take very long to determine that $\int_0^{\infty } \frac{\left(1-\exp \left(-a x^2\right)\right) J_0(b x)}{x} \, dx$ is equivalent to $\frac{1}{2} \Gamma \left(0,\...
1
vote
1
answer
254
views
How to speed up evaluation of a computationally expensive function?
The evaluation of the following code takes ages to finish:
...
1
vote
0
answers
64
views
Seek Mathematica implementation of q-analogs of Special functions
I know Mathematica has implemented some QFunctions.
What about Mathematica implementations of Elliptic gamma function, Hahn–Exton q-Bessel function, Jackson q-Bessel function,
q-exponential, q-gamma ...
6
votes
3
answers
815
views
Plotting spherical harmonics
I would like to recreate the following figure of spherical harmonics from Wikipedia.
I used the codes from Plotting real spherical harmonics with SphericalPlot3D, but the resulting spherical ...
0
votes
0
answers
134
views
How to properly use DirichletL?
I'm trying to construct in Mathematica the analytic continuation of the L-function associated with the number field Q(i). This was done, though opaquely, in the captivating YouTube series "RH ...
9
votes
3
answers
339
views
Obtaining symbolic series expansion of $(W_0(-z/e)+1)^{-1}$
Is there a way to massage Mathematica to give a symbolic expression for $s$th series coefficient of $(W_0(-z/e)+1)^{-1}$
...
2
votes
2
answers
243
views
FoxH function not working
I want to calculate the following Fox H-function:
$$
H^{1,1}_{1,2}
\left(
0.2
\,
\middle|
\begin{array}
\,
(-1, 1) \\
(-1, 1), \, (0, 0.5)
\end{array}
\right)
$$
So I tried ...
0
votes
0
answers
75
views
Conical functions with Imaginary argument
I am playing with a differential equation whose solutions are conical functions (Legendre functions with a complex order) with a pure imaginary argument
$
P_{i\nu -1/2}(i \sinh \xi)
$
or $Q_{i\nu -...
4
votes
0
answers
198
views
Why is Mathematica thinking that 1 is a variable?
I'm trying to solve an ODE and upon executing:
...
1
vote
2
answers
256
views
Simplifying expression with Erf and Erfc
I'm trying to validate the solution of an integral in a research paper, but I'm getting a different result. I get two terms of Erf and two terms of ...
2
votes
1
answer
275
views
MeijerG function with empty arguments
I am trying to find the integral from Mathematica:
...
0
votes
0
answers
18
views
`MellinConvolve` sucks
MellinConvolve doesn't evaluate to get expected result.
...
1
vote
2
answers
196
views
Why the series expansion of the complete elliptic integral of the third kind fails?
I am trying to expand the complete elliptic integral of the third kind:
...
1
vote
2
answers
216
views
I want to solve integrals in ellipsoidal coordinates?
I want to see the analytical expressions for these two integrals where s is ellipsoidal coordinates with "a" as a major axis and s>-c^2, where c is the minor axis for an ellipsoidal.
<...
0
votes
0
answers
171
views
Mistake in simplification of AppellF1
The following code
...
4
votes
3
answers
424
views
Integral in terms of 2F1(a,b;c;z) gives indeterminate answer, known well behaved solution. How do I get it to give well behaved answer? Report Issue?
I'm doing various definite integrals that I want analytic answers for (even in terms of special nonelementary functions), and I keep running into issues. Mathematica will give me an answer, but it is ...
0
votes
0
answers
175
views
Numerical Computation of Multivariate Fox-H function using given codes
This PDF is the expression of SNR of practical RIS-aided system: https://ieeexplore.ieee.org/document/9387559 (expression A.8)
Authors mentioned some sources of computation of multivariate Fox-H type ...
0
votes
1
answer
56
views
Derivative of QPochhammer
Related
Taylor expansion of a function containing QPochhammer[q, q, n]
My question
...
3
votes
2
answers
674
views
Double Integral Which is Clearly Convergent Not Converging
Mathematica has successfully computed this integral: $$\int_0^{\infty}\int_0^\infty{}dqdl(c+(q+l)^2)^{-s}$$ assuming $s>2,c>0$. The integrand comes from a first order approximation of Bessel ...
6
votes
1
answer
377
views
Did the computation of Bessel functions slow down in 14.2?
This is taking much longer to plot than I would expect:
...
0
votes
0
answers
46
views
How to get the correct result of the difference of two large numbers? [duplicate]
I found that when the variables of the Airy functions are complex numbers such as a=-36.3 - 63.1 I and b=-37.7- 65.4 I in the ...
2
votes
3
answers
264
views
Numerical integration involving Bessel function $e^{-x}I_n(x)$ with large $n$
I want to numerically integrate expressions involving $e^{-x}I_n(x)$, say
$$
\int_0^\infty \frac{1}{x\sqrt x}\Bigl(e^{-x}I_n(x)\Bigr)\Bigl(e^{-x}I_0(x)\Bigr)dx,
$$
with $I_n(x)$ the modified Bessel ...
7
votes
1
answer
311
views
Incorrect answer generated by Sum!
I tried on Mathematica on version: 13.3 and 14.1.0:
...
2
votes
1
answer
198
views
Changing the variables changes the formula result [closed]
I am intrigued why just changing the variables in an infinite series changes the result. Actually, the two different results are equivalent, but why is it affected by a mere change in variables?
The ...
0
votes
2
answers
168
views
Issue with Hypergeometric Summation in Mathematica: Discrepancy Between Sum and Analytical Expression
I am working on a symbolic summation in Mathematica involving a Hypergeometric function. Specifically, I define the following term:
...
5
votes
2
answers
411
views
Numerical approximation of implicit functions
I have an expression $t \equiv t(x)$,
\begin{equation}
t = 2^{-2+\frac{a}{2}} c^{1-\frac{a}{2}} \Gamma\left(-1+\frac{a}{2}, \frac{c}{2} x^2\right)
\end{equation}
where $\Gamma\left(-1+\frac{a}{2}, \...
8
votes
2
answers
521
views
NIntegrate cannot give high precision result for a well-behaved integral
I want to obtain value of following integral to a high precision (say 30 digits),
NIntegrate[DawsonF[Sqrt[t]]^2/t, {t, 0, Infinity}]
Graph of integrand looks like ...
0
votes
0
answers
73
views
How do I efficiently generate, store and read in Gaunt(-like) coefficients from a file?
I'm calculating the mode-mixing between various spin-weighted spherical harmonics (SWSHs) $_sY_{l}^{m}(\theta,\varphi)$. In what follows the indices can take the following values, $l_{max}$ being some ...
0
votes
0
answers
75
views
Ability of Integrate[ ] to try changes of variable on its own
I have a question about Mathematica's ability to try changes of variable when performing symbolic integration. My example is
$$ \int_0^\infty dx\ x^n \exp \left( -h \sqrt{x^2 + a^2} \right) = \frac{2^{...
0
votes
1
answer
155
views
Converting HurwitzZeta function to PolyGamma function
A generalization produces a result in terms of Zeta[s,a] function, which can be converted to PolyGamma[s-1,a] using the ...
0
votes
3
answers
372
views
How to compute the Jacobian matrix using Mathematica [duplicate]
Apologies if this is posted to the wrong forum. Have been using Maple and Maxima for 20+ years, and have only recently started using Mathematica (v. 14). Having a heck of a time getting Mathematica to ...
0
votes
1
answer
143
views
Issue with Reduction of Complete Elliptic Integral of the Second Kind
I am attempting to reduce the following equation:
y ==
I’ve entered it to be reduced as such, where L = 3.95:
By what means may this be properly reduced for y?
Both WolframAlpha and Desmos provide ...
0
votes
1
answer
118
views
Compute integrals in singular integral equation
I'm looking at this paper
https://arxiv.org/abs/2404.07307
and in particular I am interested in eqs. (16) (17), meaning I'd like to check the validity of (16) by insert it in (17).
So I'd like to ...
1
vote
2
answers
259
views
Solving ODE with DiracDelta on the right-hand side
I want to use Mathematica to validate the solution to this differential equation
$ (E+ \frac{\hbar^2}{2m} \frac{d^2}{dx^2}-\kappa x)G(x,y) = i \hbar \delta(x-y) $
The solution is supposed to be
$$ ...
2
votes
1
answer
127
views
How to force Mathematica to evaulate some values of LerchPhi function?
Mathematica cannot give LerchPhi[1,0,1]=-1/2, LerchPhi[1,1,1]=0, and LerchPhi[1,1,1/2]=0, as ...
0
votes
1
answer
143
views
A hypergeometric series function
Maple and Mathematica are very efficient to find closed form of a hypergeometric function as for example:
$$\sum _{k=0}^{\infty } \frac{(-4)^k z^k \Gamma \left(1 k+\frac{1}{6}\right) \Gamma \left(\...