Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

3
votes
1answer
84 views

Range of summation in simple Plot seems off

I was trying to reproduce a picture in a book by Havil of the sum, $$s = \sum_{r=1}^{\infty}\frac{\mu(r)}{r}\left(Li(x^{\rho_k/r})+Li(x^{\rho_k*/r})\right) $$ using ...
1
vote
2answers
97 views

Numerical integration of modified bessel function

I need to compute the following integral: NIntegrate[ BesselI[-nu, k x]/x ,{x, r1, r}] in which nu=-(2m-1)/2 and I have to ...
3
votes
0answers
57 views

Refining a density plot of the Eisenstein series argument

Thanks to amazing code from "Guess who it is" here: Eisenstein Series in Mathematica? I'm able to make some nice plots using Eisenstein Series. What I'd like is a color plot of the argument of certain ...
9
votes
1answer
149 views

Eisenstein Series in Mathematica?

Mathematica doesn't seem to have built-in tools to deal with the Eisenstein series: $$\begin{align*} E_{2}(\tau)&= 1-24 \sum_{n=1}^{\infty} \frac{n e^{2 \pi i n \tau}}{1-e^{2 \pi i n \tau}}\\ ...
1
vote
0answers
51 views

Odd plotting/math issue (could be a precision problem)

I've got a pretty odd error on a project I'm working on and was hoping to enlist some advice to fix it. The goal of this notebook is to show that I can eliminate the non-normalizable (blowing up part) ...
5
votes
1answer
151 views

Inverse error function

I solved some equation in Mathematica and I obtained something like $$y(t)=\exp \left\lbrace \left[ \text{erf}^{-1} (\text{i}t) \right]^2\right\rbrace, (1)$$ where $\text{i}$ is imaginary unit and ...
2
votes
3answers
127 views

Solving equation containing Erf expressions

Given the equation below, how do I find the value of b so that the function is equal to 21. I tried solve but I keep getting an error message. ...
4
votes
4answers
244 views

How to stop a summation when a variable is small enough?

I meet with a problem. I hope to get an infinite summation of $f1(x)/f2(x)$ which converges to zero. So my code is Sum[f1[x]/f2[x],{x,Infinity}] or ...
22
votes
8answers
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 ...
6
votes
6answers
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 ...
2
votes
0answers
90 views

Nasty integral advice

I have a pretty ferocious integral to solve, and since it doesn't seem I'll be able to do much analytically, I've taken to Mathematica to get some information. Mainly, I want to see if there are any ...
2
votes
1answer
127 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 ...
1
vote
1answer
42 views

problem about Root and Hypergeometric2F1

See these example, why does the output is different? ...
5
votes
2answers
128 views

Plot of The RiemannSiegelZ Function

I would appreciate your help to visualize the of the following function RiemannSiegelZ[x] with this range. { 18154980120849865 , 18154980120849885 } I tried this: ...
6
votes
2answers
175 views

Number of divisors visualized with the QPochhammer function, how improve performance of code?

I have this code that is originally Jeffrey Stopple's code for the Riemann zeta function in the complex plane. Because I discovered yesterday that the number of divisors can be generated with the ...
0
votes
1answer
30 views

Sum of hypergeometric functions, variable number of arguments

How can I write in Mathematica an expression like this? $$\sum_{k=1}^n {_{k}F_{k}} (1,1,\dots,1; \,2,2,\dots,2; \,z) ~,$$ where ${_{p}F_{q}}$ is the generalized hypergeometric function. My problem ...
9
votes
1answer
184 views

Wrong Limit with LaguerreL

Bug introduced in 7.0 and fixed in 10.2.0 Limit[(n! LaguerreL[n,-1]/(n^(n+1/4)/Sqrt[2]/E^(n-2 Sqrt[n]+1/2))-1) Sqrt[n], n->∞] Mathematica (wrong) output ...
20
votes
1answer
541 views

Incorrect result from Integrate

Bug introduced in 8.0 and fixed in 10.0 I attempted to calculate the following integral: ...
7
votes
2answers
108 views

Erfc Not Returning Results Specified in Documentation

In the documentation for Erfc (under "Possible Issues"), the following command returns a number that is extremely close to 2: However, when I run this same command in a fresh kernel, I get: ...
8
votes
1answer
596 views

Numerical contour integrations in the complex plane - contour deformation gives different answer for analytic kernel

I am trying to do a contour integration in Mathematica numerically. In particular, I'm checking the identity: $$ H_m^{(1)}(z) =\frac{i^{-m}}{\pi}\int_{-\pi/2 + i \infty}^{\pi/2 - i \infty} \exp[i m ...
0
votes
0answers
32 views

Timing of associated Legendre polynomials

I encountered a strange issue with the associated Legendre polynomials implemented with LegendreP[l,m,z]. Quite simply, the time used for the numerical computation of those quantities depends on ...
4
votes
1answer
197 views

Problems with series of generalized hypergeometric functions

I have been trying to do the following series: Series[HypergeometricPFQ[{1, 4, 4}, {4 - Sqrt[3], 4 + Sqrt[3]}, z],{z,1,0}] Mathematica says that the result is ...
6
votes
2answers
156 views

How does Mathematica calculate LaguerreL

About the function LaguerreL[n,a,x], the helping documents in Mathematica only say that this function satisfies equation $xy^{\prime\prime}+(a+1-x)y^\prime+ny=0$. ...
3
votes
2answers
346 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 ...
5
votes
1answer
71 views

Spurious infinite limit of integration message in NIntegrate

Bug introduced in 7.0 and persisting through 10.2 NIntegrate returns an error complaining that ComplexInfinity is an invalid ...
6
votes
1answer
89 views

`FindSequenceFunction` for sum of hypergeometric terms?

The built-in function FindSequenceFunction is quite good at recognizing hypergeometric terms, i.e. terms $c_k$ for which $c_{k+1}/c_k$ is a rational function of ...
3
votes
0answers
79 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) } ...
1
vote
0answers
127 views

Mathematica not evaluating q derivative of Jacobi theta function

Jacobi theta functions, $\theta_a(u,q)$ for $a=1,2,3,4$ are defined in the unit disk $|q|<1$. For some reason that I would like to understand, Mathematica does not evaluate numerically the $q$ ...
12
votes
1answer
703 views

Numerically solve 2nd order differential equation with singularity

Consider a second order differential equation with a potential that diverges at some generic value in the variable. For example: $$-y^{\prime\prime}(s)+\frac1{\mathrm{cn}{(s\mid k^2)}}y(s)=0$$ where ...
3
votes
3answers
257 views

3D Gamma function

I try to plot the following function: Plot3D[Gamma[1+0.5*(n+m)]/Sqrt[Gamma[1+n]*Gamma[1+m]],{n,0,1000},{m,0,1000}] I expect that for m=n and near to it the value ...
28
votes
1answer
962 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 ...
7
votes
2answers
772 views

Mathematica 10 cannot solve definite integral [duplicate]

Bug introduced in 10.0 and fixed in 10.0.2 Mathematica 10 fails to solve the following integral, saying that it does not converge. ...
3
votes
1answer
116 views

Bug in associated Legendre Polynomials?

Mathematica's definition of the connection of associated Legendre polynomials with $m$ and $-m$ is: $P_l^{-m}=(-1)^m \frac{(l-m)!}{(l+m)!} P_l^m$. We also now that $|m|>l \Rightarrow P_l^m=0$. ...
2
votes
1answer
124 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
1
vote
0answers
51 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
208 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
355 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
42 views

Inverse of a Digamma/ Polygamma Function

How one can find the inverse of a digamma/ Polygamma function in Mathematica 10?
6
votes
2answers
61 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 ...
7
votes
2answers
210 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
175 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
265 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
412 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
114 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
531 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
116 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 ...
4
votes
1answer
101 views

Complex LogIntegral error

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