Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

3
votes
1answer
38 views

Interval arithmetic for DawsonF (or other special functions)

I am currently trying to estimate a complicated expression involving DawsonF using interval arithmetic. The interval arithmetic is partially supported by ...
31
votes
0answers
1k 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 ...
16
votes
0answers
215 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 ...
8
votes
0answers
69 views

$\tt SphericalHarmonicY$ does not seem to be an eigenfunction of the spherical harmonic equation

I applied the spherical harmonic equation on the SphericalHarmonicY functions like this: ...
8
votes
0answers
729 views

Inverse Laplace transform not obtained

I can't seem to be able to compute the inverse Laplace transform of a Laplace transform: ...
6
votes
0answers
134 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 ...
6
votes
0answers
374 views

Reproducing the Integral Definition of the Modified Bessel function

I need to simplify some integral expressions in terms of special functions, such as the modified Bessel function of the first kind. See for example Eq. (5) on this page. Notice that the real ...
5
votes
0answers
115 views

Expansion of $E(i c \mid m)$ at $c\to\infty$?

Currently, I am using a Windows machine with Mathematica 8. I noticed a difference in a series expansion of the function EllipticE[] in comparison with a result ...
4
votes
0answers
89 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) } ...
4
votes
0answers
55 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
0answers
304 views

Spherical harmonic derivative

Consider the following substitution Derivative[2, 0][S][th, ph] /. S -> Function[{th, ph}, SphericalHarmonicY[3, 0, th, ph]] which gives correct answer. While ...
3
votes
0answers
71 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 ...
3
votes
0answers
77 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 ...
3
votes
0answers
171 views

Why doesn't Log[Gamma[]] simplify to LogGamma[] where it could?

I have been playing with various equations involving amount of permutations in relatively large sets. Easiest way to look at these is something like Log[10, bignumber!] . Often expressions, even ...
2
votes
0answers
96 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
0answers
169 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
0answers
239 views

How to prevent simplification of hypergeometric functions resulting from integrations?

Definite integrals from 0 to Infinity over a product of two hypergeometric (including exponential, trigonometric, hyperbolic, ...
2
votes
0answers
126 views

RSolve and incomplete gamma function

I am interested in how Mathematica uses the incomplete gamma function to solve difference equations. For example, if we have this inhomogeneous, 2nd order equation and use ...
1
vote
0answers
74 views

How to understand the symbolic integral in a result returned by DSolve?

I am trying to solve an ODE like $$\frac{dF}{dx}=\frac{a}{\ln^3 x+[c+d(ax+b)][2\ln x+c+d(ax+b)]\ln x},$$ where $a,b,c,d$ are constants. I guess it could contain a special function, say, the ...
1
vote
0answers
39 views

Problem on special functions

What does the symbol PolyLog^{(0,1)}(0,1/e) mean? I know the meaning of the Polylogarithm, but what is that exponent? It happens the same with the Lerch zeta function!! Thanks in advance.
1
vote
0answers
59 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) ...
1
vote
0answers
58 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: ...
1
vote
0answers
79 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; ...
1
vote
0answers
71 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 ...
1
vote
0answers
138 views

Converting a hypergeometric function to a Bessel function

Is there any way to get a Bessel function equivalent to a Hypergeometric function? Would it be possible in general? From Maple integration (as discussed in Integrating a BesselJ integrand to obtain ...
1
vote
0answers
207 views

Minimizing a functional takes forever

I need help with minimizing a functional in Mathematica. I have a function $V(\xi)=\sum_{i=1}^\infty C_{3_i}J_0(\xi/C_i)+C_{4_i}Y_0(\xi/C_i),~~~\Sigma\leq\xi\leq\xi_1$, and want to find such ...
0
votes
0answers
52 views

Solving Fredholm equation of the 2nd kind

While I was using the code from here to solve this integral equation: ...
0
votes
0answers
33 views

Determine class of special function from algebraic constraints?

Consider vectors x_i in arbitrary dimension. Let's say I have an expression in six variables F[x_1,x_2,x_3][x_4,x_5,x_6], which ...
0
votes
0answers
41 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 ...
0
votes
0answers
87 views

PolyLog does not simplify

How come mathematica can not simplify a simple expression involving the PolyLog function? An example is: You know that ...
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 ...
0
votes
0answers
58 views

NSolve can not handle a PolyGamma equation

I would like to solve the following equations numerically: ...
0
votes
0answers
63 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
0answers
81 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 ...
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 ...