Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

0
votes
1answer
47 views

Solving and Plotting Infinite sums over solutions to transcendental equation involving Bessel functions

I am trying to use Mathematica to plot the following equation $E(q,\Delta)$ vs $qa$, for set values of $ D \Delta/a^2$ and $\bar\rho a/D$. I have looked at other stackexchange posts which deal with ...
2
votes
1answer
67 views

Round off in Mathematica Built-in functions

Is there a way to avoid Mathematica to replace Built-in functions to other functions? For instance, the Hypergeometric1F1[a,b,x] function has a exponential form when its firsts parameters are ...
0
votes
0answers
38 views

Get asymptotic form of a Legendre function expression

How to get the asymptotic expression of the first element of the following list with small x and very large n? ...
7
votes
2answers
297 views

Can your Mathematica do the following integral?

Backslide introduced in v10 and fixed in v10.4.1. I was trying to do the following integral in Mathematica: Integrate[(z-2) PolyLog[2,z]Log[1-z]/z^3,z] What I ...
0
votes
2answers
94 views

Integral over real valued function becomes complex

I've tried to solve the following integral, but I get a complex solution even tough it should only have a real part. $$f(x)=- \frac{10^{-20} x}{0.99005- e^{10^{-12}x}}$$ Now I want to calculate the ...
9
votes
1answer
227 views

Vastly incorrect answers obtained by increasing WorkingPrecision with modified Bessel functions

Bug introduced in 7.0 or earlier and persisting through 10.4.1 This is a follow-up to this question regarding numerical instabilities occurring with modified Bessel functions. In trying to explore ...
1
vote
0answers
30 views

`QPochhammer` function simplification?

Consider the function which in Mathematica is denoted as QPochhammer[a,q,n], and its infinite product cousine: which in Mathematica reads ...
0
votes
2answers
69 views

Wrong numerical results from LegendreP

{Cos[Pi/180] // N, LegendreP[46, 0.9998476951563913`], LegendreP[46, Cos[Pi/180]] // N} give ...
5
votes
3answers
325 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: ...
11
votes
3answers
253 views

Integrating a BesselJ integrand to obtain the same result as Maple 16

I would like to check the following integration: Integrate[y Integrate[1/x BesselJ[1, x Exp[I π/4]] BesselJ[1, x Exp[-I π/4]], {x, 0, y}], {y, 0, r}] ...
33
votes
8answers
7k views

About multi-root search in Mathematica for transcendental equations

I have some questions for multiroot search for transcendental equations. Is there any clever solution to find all the roots for a transcendental equation in a specific range? Perhaps ...
3
votes
4answers
414 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 ...
5
votes
2answers
229 views

Symbolically Expanding One-Variable $q$-hypergeometric Series

In general, how do I get Mathematica to symbolically expand various infinite series to an arbitrary degree of accuracy? I deal with $q$-hypergeometric series quite frequently, and specifically want ...
9
votes
4answers
458 views

Numerical instability in cosh and sinh - integral functions [duplicate]

I'm trying to calculate the function: CoshIntegral[x] Sinh[x] - Cosh[x] SinhIntegral[x] Unfortunately Mathematica seems to hit a point (x~20) and things become ...
5
votes
2answers
250 views

Strange phenomenon occurring in analytic integration result involving Bessel functions

For the following integral, Integrate[x^2 Exp[-a x^2 - b x^4], {x, -∞, ∞}, Assumptions -> {a > 0, b > 0}] Mathematica gives the following analytic ...
0
votes
1answer
233 views

Asymptotic forms of Bessel function

I want to replace Bessel functions by asymptotic forms, so the question is: can I find the best ones with help of Mathematica? And if it's possible, how can I do it? Update How can I get with ...
5
votes
0answers
109 views

MacDonald formula for Modified Bessel Functions

How can I make Mathematica understand these two integrals? $$\int_0^{\infty} e^{-x \cosh{\xi}} d\xi = K_0(x)$$ $$\int_0^{\infty} e^{-\frac{1}{2} \Big( \frac{x y}{u} + u \frac{x^2+y^2}{x y} \Big) } ...
2
votes
0answers
35 views

Series expansion in Infinity issue with Zeta(s) function

With this code: Series[Zeta[s], {s, Infinity, #}] & /@ Range[10] // MatrixForm Series expansion for the Zeta(s) ...
2
votes
0answers
44 views

Solution of the Lagrange's minimal surface equation [closed]

Can I use Mathematica to solve the Lagrange's partial differential equation? $$(1+f_y^2)f_{xx}-2f_xf_yf_{xy}+(1+f_x^2)f_{yy}=0$$ I tried: ...
1
vote
2answers
114 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. ...
1
vote
1answer
232 views

Why Integral does not converge?

I am trying to solve for $\Omega$ this nonlinear integral equation: $$1+\dfrac{z}{k^2}-\dfrac{z^2}{K_2(z)} \dfrac{\Omega}{k^3} \displaystyle\int_{1}^{\infty} \gamma^2\, \text{ArcTanh} ...
1
vote
0answers
44 views

High Precision Plots of Eisenstein Series [closed]

When plotting the Eisenstein Series (great information here Eisenstein Series in Mathematica?) you observe highly non-trivial branch cut behavior close to the real axis. This makes the numerics break ...
3
votes
0answers
56 views

Are there any Mathematica programs on or about the LMFDB Archive?

I would like to explore the LMFDB Archive with L-functions using Mathematica. Is there any Mathematica sample code available to get me started?
0
votes
1answer
51 views

integral involving square of exponential integral

I'm trying to compute the integral $$ \int_{a}^{\infty}\frac{e^x}{x}[\mathrm{Ei}(-x)]^2\,dx, $$ where $\mathrm{Ei}$ is the exponential integral, and $a>0$. The obvious ...
1
vote
1answer
77 views

Solving two equations with modified Bessel functions

I am trying to solve two equations with Bessel functions in them, 1) C1*BesselK[0, 3.7268*10^-4*x] == 1.3*10^-6 2) ...
23
votes
2answers
1k 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: ...
1
vote
1answer
70 views

Finding roots of Bessel function $y=3J_1(x)+xJ_1'(x)$ is returning inaccurate roots. Not Kernel bug [closed]

I can't figure out why Mathematica is returning the incorrect roots. The first five should be 2.9496,5.84113,8.87273,11,9561, and 15.0624 according to my textbook. ...
3
votes
2answers
72 views

Integral invoving product of Whittaker functions

I'm trying to evaluate the integral $$ \int W_{0,a}(x)\,W_{1,b}(x)\,\frac{dx}{x}, $$ where $W$ denotes the Whittaker $W$ function, and $a,b\in[0,1/2]$. Using the Whittaker differential equations that ...
3
votes
2answers
150 views

How can I get LogIntegral[z] to be printed as “li[z]”?

I have no problem with the current formatting of the function, but for the sake of the reader less familiar with Mathematica functions, is there a way to define, say, ...
0
votes
1answer
56 views

NIntegrate[] Gamma Function

Gamma Function is known to be : Source first i plot the function z = 1; f[t_] := (t^(z - 1))/Exp[-t]; gamma = Plot[Gamma[g], {g, 1.0, 5.0}] then, i ...
1
vote
3answers
239 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
70 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
65 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 ...
4
votes
2answers
328 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 ...
6
votes
0answers
99 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
70 views

Ramanujan's asymptotic formula for $p(n)$ [closed]

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
82 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
75 views

ComplexInfinity for a convergent product

The infinite product involving the ratio of (n^2)! to its Stirling approximation ...
10
votes
1answer
142 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 ...
4
votes
0answers
82 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
44 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
95 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
104 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
27 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 ...
12
votes
2answers
817 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
299 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
62 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
181 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
57 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
74 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.