Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

2
votes
1answer
121 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
2
votes
0answers
51 views

Eisenstein Series in Mathematica?

Mathematica doesn't seem to have built-in tools to deal with the Eisenstein Series: $E_{2}(\tau)= 1-24 \sum_{n=1}^{\infty} \frac{n e^{2 \pi i n \tau}}{1-e^{2 \pi i n \tau}}$ $E_{4}(\tau)= 1+240 ...
1
vote
0answers
46 views

Use Solve[] with Bessel, gamma, and hypergeometric functions?

I need to find values of {a,b,c} such that the 0th, 2nd, and 4th order moments of f[x]=Exp[-ax^4 - bx^2 - c] will equal respectively {1,2,10}. I didn't really expect this to work, and it didn't: ...
3
votes
2answers
203 views

Storing Variables in “Loops” and Point Plotting

Given the function $y=\sin x$ defined over the region $-\pi \leq x \leq \pi$, I need to implement a "do loop" such that I sweep over 100 or so points $-1 \leq y \leq 1$ and find precisely the two $x$ ...
14
votes
1answer
345 views

How to enlarge Mathematica's knowledge about certain functions?

I'm often troubled with the following task. I need to carry out symbolical computations involving certain special functions. Let me take as an example Barnes gamma-function. It is included in ...
-1
votes
1answer
35 views

Inverse of a Digamma/ Polygamma Function

How one can find the inverse of a digamma/ Polygamma function in Mathematica 10?
5
votes
5answers
2k views

Series expansion in terms of Hermite polynomials

I am trying to expand a polynomial in terms of orthogonal polynomials (in my case, Hermite). Maple has a nice built-in function for this, ChangeBasis. Is there a ...
6
votes
2answers
59 views

Simplify expression to Integer

I have the following function: f[n_]:=(n-1) E Gamma[n,1]; I know that this expression always evaluates to an Integer, but Mathematica doesn't always output the ...
0
votes
0answers
10 views

Airy function and its value at 0 [closed]

On the page http://mathworld.wolfram.com/AiryFunctions.html, the Airy function is defined as $$ \text{Ai}(z)=\frac{1}{2π}\int_{-\infty}^{+\infty} e^{izt+i\frac{t^3}{3}} dt, $$ a formula that can be ...
-2
votes
0answers
53 views

Solving an integral equation using Mathematica [closed]

Hi I have the following problem I want to solve numerically by Mathematica the integral equation $$v(x,t)=1-0.0000065*\int_0^t \theta_0(1/2 x| i\pi \tau)v^4(1,t-\tau)d\tau$$ where $\theta_0(1/2 x | ...
7
votes
2answers
203 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 ...
9
votes
2answers
173 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}] Mathematica 9.0 is ...
3
votes
1answer
261 views

plotting hypergeometric functions

Does anyone know why a plot of a hypergeometric function turns out differently in Mathematica than in Maple? The function I'm plotting between x=-30 and +30 is: ...
4
votes
1answer
402 views

Strange result for the analytic integration leads to Hypergeometric2F1

The integration result for Integrate[1/(r^2 Sqrt[x/r^(4 - 2 \[Gamma]) + 1]), r] is: ...
2
votes
1answer
111 views

Unexpected Weierstrass $\wp$-Function Behavior

In Mathematica, the Weierstrass $\wp$-function is expressed by WeierstrassP[u + I v, {g₂, g₃}] where g₂ and ...
14
votes
4answers
529 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 ...
2
votes
1answer
110 views

What are the terms of the sequence generated by Zeta(3s)/Zeta(s)?

The LiouvilleLambda function has Dirichlet generating function of Zeta[2s]/Zeta[s]. I am curious about an analogous function with Dirichlet generating function of Zeta[3s]/Zeta[s]. Can Mathematica ...
3
votes
1answer
87 views

Complex LogIntegral error

Going through Derbyshire's Prime Obsession & trying to take LogIntegral of 20^ZetaZero[1] & comes up with a value of ...
3
votes
1answer
127 views

Strange NSum behavior

If I do: NSum[(i + 1)/(i + 2) LegendreP[i, 0] LegendreP[i, 0], {i, 0, Infinity}] I get: 1.25216 If I do: ...
3
votes
2answers
131 views

ParametricPlot with 2 variables

I want to try plotting this: As you can see the first axis is v and the second is just dependent of the angle of $\Phi$ . The function is ...
0
votes
1answer
85 views

Why is the indefinite integral $e^\frac{x}{2}x^\frac{g}{2}dx$ the upper Gamma function in Mathematica? [closed]

I have an indefinite integral $e^\frac{x}{2}x^\frac{g}{2}dx$, I try to solve in Mathematica. Solve[Int[e^{x/2}x^{g/2-1}dx ,x]] (* Gamma[g/2,x/2] *) How does ...
2
votes
2answers
305 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
251 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
112 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
165 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
932 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
240 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
58 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; ...
10
votes
2answers
245 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 ...
3
votes
3answers
134 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
246 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
476 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
128 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
358 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
97 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
201 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
205 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
211 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
155 views

Definite Integral over Bessel Function

Hello I am interested in evaluating the following integral. ...
4
votes
1answer
298 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
103 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
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
417 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
521 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
89 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
77 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
103 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 ...