Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

7
votes
1answer
108 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
48 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
207 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
35 views

Inverse of a Digamma/ Polygamma Function

How one can find the inverse of a digamma/ Polygamma function in Mathematica 10?
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 ...
2
votes
1answer
112 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
132 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
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; ...
2
votes
2answers
308 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
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 ...
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 ...
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
56 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
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
130 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
54 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
56 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
246 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
76 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
100 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
61 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
290 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
86 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 ...
7
votes
1answer
153 views

Wrong Limit with LaguerreL

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
116 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
64 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
126 views
0
votes
1answer
62 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
83 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: ...
0
votes
1answer
68 views

Orthogonal Collocation Using Jacobi Polynomials

I'm trying to solve a PDE(diffusion-reaction in a spherical catalyst pellet) using Jacobi Orthogonal Collocation method. But at the stage of solving the system of ODEs(using ...
16
votes
0answers
197 views

Is MathieuC for moderately large imaginary arguments broken?

I'm trying to plot MathieuC[-3,0.3,I x] for $x\in[0,10]$, and here's what I get even with arbitrary precision arithmetic (here I use ...
0
votes
1answer
62 views

Using Assumptions in Expressions that Evaluate to be Elliptic Integrals

I have the following integral that I am trying to evaluate in Mathematica: $\int \sqrt{\alpha + m g l \cos(q)} dq$. If $\alpha > m g l$, then the result is a complete elliptic integral of the ...
0
votes
1answer
50 views

Having trouble interpreting the results DSolve gives for the Laguerre equation

I am trying to solve a second order ODE using Mathematica. Before I get into solving my (more complicated) problem, I am trying to use DSolve on known ODEs to check that the answer that Mathematica ...
3
votes
1answer
75 views

Why minimization does not work with symbolic array as arguments

If I try to minimize with constraints a function of several variables, with Gamma regularized function involved in the constraints it seems to works, as shown below (this is just a dummy example ...
7
votes
1answer
100 views

How to define $\operatorname{Mc}$ and $\operatorname{Ms}$ Mathieu functions in terms of MathieuC and MathieuS?

Reading multiple books about Mathieu functions, I always come across notation like $\operatorname{ce}_r(z,q)$, $\operatorname{se}_r(z,q)$ for angular functions and $\operatorname{Ce}_r(z,q)$, ...
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: ...
3
votes
1answer
87 views

Why is the Spherical Bessel Function acting strangely at this point?

I'm doing some computation that requires the use of Spherical Bessel Functions of the 1st kind, at high orders and values. So, I managed to find this, while running it over a wide range of values. I ...
1
vote
1answer
77 views

Computational Complexity of HypergeometricPFQ

What is the computational complexity of HypergeometricPFQ? I got a result of a product of a multinomial and HypergeometricPFQ and I was wondering if that would be considered a closed form solution ...
7
votes
1answer
93 views

Sum over Binomials and Gammas

Given the function, ...
18
votes
4answers
463 views

Expansion of a hypergeometric function takes ages with Mathematica 9 and 10 (regression?)

Mathematica 8 (Linux version) can evaluate AbsoluteTiming[Series[Hypergeometric2F1[1, 1 - eps/2, 3 - eps, 1/2], {eps, 0, 0}]] in no time. On one of the ...
3
votes
3answers
251 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 ...
-1
votes
2answers
76 views

Complete definition of special functions in Mathematica [closed]

How can one display the actual integral involved in special functions like this MarcumQ[2, a, b] in Mathematica?
2
votes
2answers
254 views

Perturbation theory with Mathematica: Definite integral of polynomial times exponential times hypergeometric function of imaginary argument

I would like to ask also Mathematica users about my question from the math forum. To expand, I'm adding the code which calculates the full double integral for $n=0$ and $\mu=0$ (the second in the ...
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 ...
1
vote
1answer
46 views

Defining functions that depend upon other functions which depend upon other functions and so on

The problem is that I have to evaluate functions, that take output from previous functions, and combine them to form new functions, and do some operations and evaluate the value at new function. But, ...
1
vote
1answer
110 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 ...