Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

2
votes
2answers
294 views

How I can integrate $\int_0^{a} x e^{-\frac{b^2 x^2}{2 c}} J_0(n x) dx$?

How can I get a solution to the integral given below ? $\quad \quad \int_0^{a} x e^{-\frac{b^2 x^2}{2 c}} J_0(n x) dx$ where $a,\,b,\,c$ and $n$ are constants and $J_0$ is a Bessel function of the ...
3
votes
3answers
248 views

Hypergeometric function with a matrix argument

I am looking for the evaluation of a Hypergeometric function with a matrix argument as for example in Koev and Edelman or as showcased in this Wikipedia article. From what I understand from ...
0
votes
0answers
62 views

Why is PolyLog[] giving weird answers for ordinary values? [duplicate]

Possibly related to this question, but it seems slightly different: Strange behaviour of PolyLog Function Wikipedia says that for real s, z<1 should be real. So I was confused when MMa returned: ...
4
votes
2answers
107 views

Imaginary terms in the derivative of Jacobi theta function (2) on the real line

I am trying to calculate/plot the derivative of the second Jacobi theta function $d\theta_2(0, e^{-\pi t} )/dt$. Calculating or plotting the function itself works fine: ...
2
votes
2answers
159 views

Define PolyLog so that positive reals evaluate on upper edge of branch cut

Generalizing the previous question: Define Log so that negative reals evaluate on lower edge of branch For real positive values $x>1$, Mathematica's polylogarithm function ...
28
votes
1answer
907 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 ...
4
votes
1answer
236 views

Numerical errors/inaccuracies in ProductLog

Context In cosmology, a fairly accurate model to describe the gravitational potential, $\psi(r)$ of dark matter halos is given by $\psi( r)=\log(1+r)/r$. ...
1
vote
0answers
54 views

Improve speed of evaluating a sum over four indices

I am trying to implement the Fox-H function with several variables as sum of residues. I have arrived at the following function; ...
13
votes
4answers
483 views

How to plot Ramanujan's continued fraction in Mathematica?

I want to plot Ramanujan's continued fraction: $$R(q)=\cfrac{q^{1/5}}{1+\cfrac{q}{1+\cfrac{q^2}{1+\ddots}}}$$ but I do not know how to define this function in Mathematica. How do I define and plot ...
10
votes
2answers
239 views

Asymptotics of $\frac{\sum _{i=0}^{\lfloor n/2 \rfloor} {2(n-2i) \choose n-2i} {n \choose 2i} {4i \choose 2i}}{2^{3n - 1}}$

I am fairly sure that asymptotically $$\frac{\sum _{i=0}^{\lfloor n/2 \rfloor} {2(n-2i) \choose n-2i} {n \choose 2i} {4i \choose 2i}}{2^{3n - 1}} \sim \frac{2}{\pi n}.$$ I tried ...
2
votes
1answer
106 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
3
votes
3answers
133 views

Calculating $\lim_{x\to 1} \, \int -\frac{i \text{Li}_2\left(x-x^2\right)}{\sqrt{3} \left(x-\frac{1}{2}-\frac{i \sqrt{3}}{2}\right)} \, dx$

How can we force Mathematica to compute this limit? $$\lim_{x\to 1} \, \int -\frac{i \text{Li}_2\left(x-x^2\right)}{\sqrt{3} \left(x-\frac{1}{2}-\frac{i \sqrt{3}}{2}\right)} \, dx$$ It seems it ...
5
votes
2answers
245 views

When does the real part of Zeta vanish on the critical line?

This seems to be a quite a simple problem but I cannot make it work. I am trying to find all values within a given range for which the real part of the Zeta ...
1
vote
1answer
470 views

Plotting Fresnel function

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

How can I compute the real part of $\zeta^2$ numerically? [duplicate]

I want to compute and plot $\Re(\zeta(x+iy)^2)$ and $\Im(\zeta(x+iy)^2)$. How can I do that with Mathematica?
0
votes
1answer
344 views

How to use `FindRoot` to solve an equation containing a parameter?

I'm trying to derive some of the results of the following paper: Electrodynamics of semiconductor-coated noble metal nanoshells, JT Manassah - Physical Review A In the paper there is matrix $\mathbf ...
4
votes
2answers
83 views

Evaluating a Series expansion of PolyLog function

I am trying to evaluate an expansion of the following integral $$ \int_0^\infty \frac{p^4}{1+\exp\left({\frac{p^2}{2mT} - \frac{\mu}{T}}\right)}\, dp = A_0(m,\mu) + A_2(m,\mu)T^2 + \ldots $$ in terms ...
0
votes
1answer
192 views

How to calculate the unknown quantity in an infinite series?

I'd like to calculate x value in this equation. Basically, I tried to 2 types of method which are FindRoot and NSolve. But, I have failed the calculation caused by these errors up to now. If ...
5
votes
1answer
198 views

How to plot the result of this singular integral?

Please I open a new post here after this one : http://mathematica.stackexchange.com/a/59203/10158 Now I want to plot the function $f(a,b)$ as a function of $b$ for different values of $a$ : $a=0.5$ , ...
1
vote
1answer
204 views

Plotting Integral equation

I want to plot the following indefinite integral : $C_l^{CC}=\int k^2\mathrm{d}k\: [e^{-2k^{2}}P_{Cl}^2(k\eta)|\dot{h}(\eta)|^2]$ with k from 0 to some large value (considered to be $\infty$), ...
2
votes
1answer
152 views

Definite Integral over Bessel Function

Hello I am interested in evaluating the following integral. ...
4
votes
1answer
287 views

Fredholm integral equation of the second kind with kernel containing Bessel and Struve functions

I need to solve this Fredholm integral equation of the second kind: f[s]+integrate[f[t] K[s,t],{t,0,1}]=s where ...
21
votes
7answers
5k 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 ...
1
vote
2answers
87 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 ...
0
votes
0answers
8 views

Scaling properties of Airy functions. [migrated]

I'm specifically intereseted in how to rewrite $(\alpha\in \mathbb{R})$ $$\mathrm{Ai}(\alpha\cdot z)$$ as some constant times $\mathrm{Ai}(z).$
3
votes
1answer
55 views

Generating terms of the Stirling series

The Stirling series starts as follows: $$ n!=\left(\frac{n}{e}\right)^{n}\sqrt{2\pi n} \Bigl\{1+\frac{1}{12n}+\frac{1}{288n^{2}}-\frac{139}{51840 n^{3}}-\frac{571}{2888380 n^{4}}+O(n^{-5})\Bigr\}. $$ ...
0
votes
1answer
404 views

Trying to solve a transcendental equation involving Bessel functions

I've never used Mathematica before and am trying to numerically solve equation (12) from this paper. Ideally I'd be able to find the smallest value of $x_{n\nu}$ for $\exp(-kr\pi)$ close to 1, and ...
6
votes
1answer
513 views

Integrating over Bessel Function erroreous? (Hankel Transform)

Bug introduced in 8.0.4 or earlier and persists through 10.0.2. The Hankel Transform is given by Integrate[f[x] x BesselJ[0,x t],{x,0,Infinity}] It is ...
0
votes
0answers
88 views

Symbolic representation of bessel series derivative

I want to get symbolic expression for BesselJ derivative, where BesselJ is represented as series: ...
3
votes
0answers
75 views

speed up evaluating a listable function

I apologize in advance for the vagueness, but I can't think of a more descriptive title. I am trying to find the average value of the square of BesselJ[0, x] ...
0
votes
0answers
99 views

Converting a hypergeometric function to a Bessel function

Is there any way to get a Bessel function equivalent to a Hypergeometric function? Would it be possible in general? From Maple integration (as discussed in Integrating a BesselJ integrand to obtain ...
0
votes
1answer
266 views
1
vote
2answers
85 views

Integrating a compound expression

I have an integral of the form I[r]=∫(arExp[-r]-brSin[k(r-d)]Exp[-r])BesselJ[0,kr]dr where Besse1J[0,kr] is the modified ...
0
votes
0answers
62 views

Spherical Hankel and Bessel in Explicit form

This is probably an easy question, when i type the spherical Hankel (first kind) and the bessel function into WolframAlpha it gives back an explicit form, the one you would get if you were to do it by ...
11
votes
3answers
372 views
2
votes
1answer
114 views
0
votes
0answers
52 views

Plotting complex roots of equation involving special functions

I consider an equation of the kind $u(k,\lambda)=0$, where $u$ involves special functions, ParabolicCylinderD, and $k$ and $\lambda$ are real parameter. By implicit ...
0
votes
0answers
55 views

PolyLog does not simplify

How come mathematica can not simplify a simple expression involving the PolyLog function? An example is: You know that ...
16
votes
2answers
288 views

How to improve performance of BesselJ to the level of GSL?

Consider the following code: zs = N /@ Range[0, 12, 10^-5]; AbsoluteTiming[bessels = BesselJ[1, #] & /@ zs;] Length@zs I've tried to measure only computation ...
1
vote
2answers
96 views

Trouble using Solve and NSolve with functions involving Erf

I have the following functions: R[k_, x_, t_] := -.5*(k - x)*(1 + Erf[-(k - x)/t]) L[k_, c_, x_, t_] := .5*c*(k - x)*(1 + Erf[(k - x)/t]) I'm interested in ...
3
votes
1answer
243 views

Simplify doesn't simplify HypergeometricPFQ with exact arguments

Consider a series: $$\sum_{t=0}^\infty \frac{8^{-11-2t}(22+4t)!}{t!(11+t)!(11+2t)!(32+t)}$$ ...
4
votes
1answer
135 views

What is the formula Mathematica uses for ZetaZero?

What is the formula/algorithm Mathematica uses for the ZetaZero command?
1
vote
0answers
73 views

Bessel functions and ListContourPlot [closed]

I'm trying to plot some eigenfunctions (for a semicircle) in Mathematica. Functions are of the following form: $$g_{ms}=J_{2m+1}(y_{ms}x) \sin ((2m+1)\varphi)$$ For example, for $m=0$ and $s=1$ it ...
0
votes
0answers
18 views

Not Simplify automatic some expression(values)

I try to calculate $$\text{res}\left(h(z) (\psi ^{(0)}(z)+\gamma )^2,\{z,-4\}\right)$$ the result is $$\frac{1}{6} \left(6 h'(-4)-25 h(-4)\right)$$ the result it is possible to write as $$2 \left(6 ...
1
vote
1answer
103 views

Is there an easy way to let mathematica print out every Erfc and InverseErfc as F and F^{-1}

Mathematica uses complementary error function and its inverse as functions for example when integral of a Gaussian is taken. Therefore, all output expressions of Mathematica involve Erfc and ...
4
votes
1answer
289 views

Simplifying numerical expressions involving special functions

I've encountered the following problem. There is the identity (Legendre's relation) that the special functions EllipticK[x] and ...
6
votes
0answers
108 views

Fine tuning compiled code that computes dilogarithm function

As an exercise of writing a good Compile function, I want to do the simple task of coding a routine that outputs the real part of the dilogarithm function ...
1
vote
0answers
50 views

Plotting the Schwarz D minimal surface

With the aid of the following paper I'm trying to plot the Schwarz D minimal surface: paper. So far I have followed all the instructions to the best of my abilities and have come up with the following ...
0
votes
0answers
84 views

Computing Poincaré symbolic solution for an arbitrary integer order polynomial

In the 1880s, Poincaré created functions which give the solution to the nth order polynomial equation in finite form. These functions turned out to be "natural" generalizations of the elliptic ...