Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

2
votes
1answer
65 views

Plotting a partial sum

I am given the Legendre expansion of the first kind. $$f(x)=\sum_{n=0}^{\infty}A_{n}P_{n}(x).$$ I have worked the coefficient to be $$A_{n}=\frac{1}{\left \| P_{n}(x) \right \|^{2}} ...
4
votes
0answers
41 views

Expanding PolyGamma function error

For large half-integer arguments of PolyGamma[] function the FunctionExpand[] or ...
1
vote
0answers
46 views

Legendre expansion of functions [on hold]

I am given the expansion for the Legendre function of the first kind. It is just $$f(x)=\sum_{n=0}^{\infty}A_{n}P_{n}(x).$$ I have found the expression for the coefficient $A_{n}$. However I would ...
3
votes
1answer
88 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 ...
3
votes
0answers
60 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 ...
1
vote
0answers
54 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) ...
4
votes
4answers
245 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 ...
1
vote
1answer
42 views

problem about Root and Hypergeometric2F1

See these example, why does the output is different? ...
2
votes
0answers
91 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 ...
5
votes
2answers
129 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: ...
0
votes
1answer
32 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 ...
1
vote
0answers
44 views
0
votes
0answers
33 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 ...
6
votes
2answers
157 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
1answer
86 views
5
votes
1answer
72 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 ...
9
votes
1answer
153 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
53 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
210 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$ ...
-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 ...
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 ...
0
votes
1answer
86 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 ...
3
votes
2answers
144 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 ...
1
vote
0answers
60 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; ...
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 ...
3
votes
3answers
136 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 ...
4
votes
2answers
110 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 ...
1
vote
2answers
117 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
57 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\}. $$ ...
3
votes
0answers
78 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] ...
14
votes
4answers
532 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
139 views

How can I find all solutions of a complicated expression in a desired interval with Mathematica?

I have an expression which is a combination of Bessel functions: ...
0
votes
0answers
58 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
63 views

PolyLog does not simplify

How come mathematica can not simplify a simple expression involving the PolyLog function? An example is: You know that ...
10
votes
2answers
249 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 ...
1
vote
0answers
80 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 ...
2
votes
1answer
103 views

Mathematica and special functions

How is it possible that Mathematica doesn't recognize it's own definitions of special functions ? I tried as input: ...
0
votes
0answers
70 views

Solving a Transcendental Function of two variables numerically

I have a transcendental function F of two variables, i.e. F(x,y). The function F contains ...
16
votes
2answers
295 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 ...
0
votes
1answer
74 views

What is a good way to check what kind of a function a function is?

i want a module or function or switch or If (whatever works), what asks if a function is an e function or a polynomial(a hole number one or fractorial), a exponential function, a logarithm function ...
0
votes
1answer
90 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 ...
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 ...
3
votes
1answer
118 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 ...
1
vote
0answers
65 views

Orthogonality relations of Hermite polynomials

The Hermite polynomials are orthogonal. $$ \int_{-\infty}^\infty H_m(x) H_n(x) e^{-x^2}\, \mathrm{d}x = \sqrt{ \pi} 2^n n! \delta_{nm} $$ Does Mathematica not use this relationship? Because running ...
-1
votes
1answer
127 views
0
votes
1answer
63 views

Problem with Erf near Infinity in NSolve

I have to solve a differential equation involving the error function Erf, but my code is not able to evaluate it in the range I need. I think I can formulate my ...
1
vote
1answer
87 views

Work around bugs in Summation, Hypergeometric function?

I'm so confused I don't even know how to phrase the question, so here's what's happening: I need an analytic form for this: ...