Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

1
vote
2answers
119 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
125 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
61 views

NSolve can not handle a PolyGamma equation

I would like to solve the following equations numerically: ...
2
votes
1answer
151 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
102 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
80 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
68 views

Does Mathematica support Laguerre Polynomials of Matrix Argument?

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
192 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
214 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
133 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
103 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
231 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
215 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
138 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
335 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 ...
4
votes
0answers
94 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
906 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
259 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
141 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
261 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
44 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
305 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
118 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
164 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
190 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
235 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
122 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
207 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 ...
2
votes
1answer
71 views

Improving working precision of LegendreP[n,x]? [duplicate]

I was trying to evaluate N[LegendreP[5,0.1]] The cell gives me: N[LegendreP[5,0.1]]=0.178829 However I wanted more ...
5
votes
1answer
253 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$ , ...
8
votes
2answers
291 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 ...
2
votes
2answers
608 views

Find roots of a function involving Bessel functions

I'm trying to find the roots of a function involving Bessel functions. Here is my code ...
2
votes
1answer
261 views

Find solutions of equation involving Bessel functions

I'm new in Mathematica and I'm trying to find the solutions of this equation involving Bessel functions $$\eta \frac{ J_{n+1}(\eta a)}{J_n(\eta a)}+\chi \frac{ I_{n+1}(\chi a)}{I_n(\chi a)}=0$$ ...
9
votes
1answer
239 views

Mathieu function periodicity problem

According to the documentation, the Mathieu characteristic function generates parameter a: The characteristic value Subscript[a, r] gives the value of the parameter a in y′′+(a-2q cos(2z))y=0 ...
7
votes
2answers
238 views

Why does N[Re@f] give complex result?

Consider this code: ...
3
votes
3answers
305 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 ...
6
votes
1answer
211 views

Numerical error in Mathieu functions

Consider the MathieuCharacteristicA function, which is a piecewise function according to the documentation. The discontinuity happens at integer number. ...
4
votes
1answer
434 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 ...
3
votes
1answer
93 views

Hypergeometric function with large parameters [duplicate]

I need an efficient and accurate method to evaluate hypergeometric ratios of the form: $$\frac{_{2}F_{1}(a+1,b;c;x)}{_{2}F_{1}(a,b;c;x)}$$ for large positive values of ...
6
votes
1answer
169 views

Symbolic integration of elliptic functions

Is there a clever way to integrate products of elliptic functions like $\wp(z;g_2,g_3)$ or $\zeta(z;g_2,g_3)$ in Mathematica? Mathematica seems to be able to integrate functions like ...
9
votes
1answer
333 views

Expanding derivatives of hypergeometric functions

Sometimes Mathematica expresses results of integration or summation in terms of symbolic derivatives of Hypergeometric2F1 function, and cannot further simplify ...