Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

1
vote
1answer
23 views

How to solve an iterative integral equation for a function in Mathematica 8.0?

I would like to solve the following iterative integral equation in order to find the functional form of $P_{0}(k)$ numerically and then use it subsequently for the rest of my code: ...
0
votes
1answer
48 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 ...
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? ...
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
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 ...
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 ...
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 ...
5
votes
2answers
251 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 ...
2
votes
0answers
36 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
0answers
45 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) ...
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, ...
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. ...
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} ...
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 ...
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 ...
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 ...
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 ...
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 ...
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$, ...
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: ...
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 ...
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: ...
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
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.
1
vote
2answers
93 views

Asymptotic form of the “strange” function

I want to find the asymptotic form of this function ...
0
votes
0answers
60 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 = ...
1
vote
0answers
37 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}% ...
0
votes
0answers
37 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 ...
5
votes
0answers
76 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 ...
2
votes
1answer
49 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 ...
10
votes
1answer
143 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 ...
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. ...
3
votes
2answers
87 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: ...
1
vote
1answer
47 views

Unexpected behavior in symbolic integration with GenerateConditions->False

Consider the following two symbolic integrations: ...
4
votes
0answers
50 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 ...
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
96 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
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 ...
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 ...