Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

8
votes
1answer
127 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
133 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
106 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
109 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
113 views

Sum over Binomials and Gammas

Given the function, ...
20
votes
4answers
600 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
326 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
88 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
336 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
143 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
160 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 ...
2
votes
3answers
156 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. ...
1
vote
2answers
264 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 ...
1
vote
2answers
122 views

Function to Represent Recursive Integral

I'd like to represent the following recursive integral equation to evaluate/graph (for $n\leq3$): $$K_{i,n}\left(x\right)=\int_{x}^{\infty}K_{i,n-1}\left(y\right)dy$$ where ...
1
vote
1answer
126 views

How to find and verify relationships between functions?

MMa gave me a complicated result involving Hypergeometric0F1's. It was much less complicated after I discovered this identity: ...
0
votes
0answers
62 views

NSolve can not handle a PolyGamma equation

I would like to solve the following equations numerically: ...
2
votes
1answer
153 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 ...
0
votes
1answer
109 views

FactorialPower and Factorial

After some computation, I have obtained a function FactorialPower[1, n, -1]. Clearly, FactorialPower[1, n, -1] equals ...
1
vote
0answers
83 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
72 views

Does Mathematica support Laguerre Polynomials of Matrix Argument? [duplicate]

I am wondering whether Mathematica can also compute the Laguerre polynomials of matrix argument as the appear in multivariate statistics? For example, they appear in this paper here: ...
0
votes
2answers
195 views

NIntegrate and Integrate of a Spherical Bessel function

I am trying to integrate over a spherical bessel function. I have used both the Integrate and NIntegrate functions in Mathematica but the values given by each do not match. Any reason why this ...
2
votes
2answers
220 views

Probability of multivariate normal being positive on each coordinate

How can I find the probability that each coordinate of a specified multivariate normal distribution is positive? I tried the following, which I believed should work ...
3
votes
1answer
136 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$. ...
0
votes
1answer
57 views

EllipticF reduction

Is it possible to reduce the following Z to Legendre form Elliptic integral of the First kind? $$ \int \dfrac {\sec u\; du } {\sqrt{ (1- {(\nu \tan u)}^2 }} ..(1*)$$ With ...
0
votes
0answers
63 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: ...
1
vote
1answer
108 views

Strange timings of integrals involving Hermite's polynomials

I have used Mathematica to calculate tunneling for quantum harmonic oscillator. The code is simple: ...
1
vote
2answers
105 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 ...
8
votes
2answers
236 views

Number of divisors visualized with the QPochhammer function, how to 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 ...
2
votes
2answers
112 views

Using `Fold` to show stages of Euler product formula

I would like to recursively replicate each stage of the Euler product proof. This does it: ...
1
vote
1answer
109 views

Is there a LogBeta function like the LogGamma?

In a computation, I need Log[Beta[alpha+j,beta-j+n]] of some kind. Is there a LogBeta function built-in to avoid any under/over ...
0
votes
1answer
223 views

Get rid of Error Function: How to get rid of sequential appearances of error function?

We have a function as e[t_] :=(E^(-t^2)) Cos[0.1 t] and we must evaluate below integration (However I used the variable x ...
4
votes
2answers
142 views

Inverse of LogIntegral

I wanted the inverse logarithmic intgral, so I typed InverseFunction[LogIntegral] and received the expected symbolic answer. But when I try to integrate it or ...
5
votes
2answers
344 views

Lower branch of Lambert W function in mathematica

I am interested in values of the Lambert W function, which is defined as the solution to the equation $ z = W(z) e^{W(z)} . $ The solution is not, however, single-valued, but branches into two ...
0
votes
0answers
87 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 ...
5
votes
0answers
109 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) } ...
0
votes
2answers
112 views

Is this answer true?

I was using Mathematica to get the series solution for Legendre equation. But when I get the recurrence relation and use RSolve: ...
7
votes
2answers
909 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. ...
4
votes
1answer
266 views

Root finding: zeroes of Mathieu function

I am finding the roots of the Mathieu sine function, and find Mathematica and Maple do not agree on the solutions. For example, consider the solutions of ...
1
vote
2answers
145 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 ...
5
votes
1answer
279 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 ...
0
votes
0answers
45 views

Should this be equal to Gamma function or not?

I define the following function: Nat2[s_] := InverseFourierTransform[FourierTransform[1/t, t, w] (-I w)^(s - 1), w, x]/Cos[Pi (s - 1)] /. x -> 1 I expected ...
1
vote
1answer
316 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$), ...
6
votes
1answer
122 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 ...
4
votes
1answer
97 views

FreeQ and arguments of Hypergeometric2F1

I was trying to use FreeQ to test for the presence of Hypergeometric2F1 functions in my expressions. I encountered the following ...
6
votes
0answers
167 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 ...
2
votes
0answers
193 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$ ...
2
votes
1answer
240 views

Definite Integral over Bessel Function

Hello I am interested in evaluating the following integral. ...
4
votes
0answers
57 views

Example in Help File does not evaluate as claimed

In the help file, under BellB, I read at "Properties and Relations": Sum can give results involving ...
4
votes
1answer
123 views

Possible bug / numerical issues with HypergeometricU — any suggestions for a fast workaround?

I've encountered some problematic behaviour with HypergeometricU. I have a probability distribution on the positive integers that takes the following form after ...
0
votes
1answer
213 views

Solution of differential equation in terms of incomplete gamma function

I need help in solving equation 15 and 16 either manually or in Mathematica to get the solution in terms of the incomplete gamma function. This is what Mathematica tells me. I can't understand ...