Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

4
votes
2answers
107 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: ...
28
votes
1answer
906 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 ...
2
votes
1answer
106 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
16
votes
0answers
192 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
640 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
108 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
217 views

Using FunctionExpand to evaluate symbolic derivatives

Some symbolic derivatives of certain special function are not expanded automatically, but FunctionExpand often helps to get a derivative-free closed form ...
6
votes
0answers
346 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
69 views

FindSequenceFunction for sum of hypergeometric terms

Mathematica's 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 ...
5
votes
0answers
109 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
93 views

Symbolic integration of elliptic functions

Is there some clever way to integrate products of elliptic functions $\wp(z;g_2,g_3)$ or $\zeta(z;g_2,g_3)$ in Mathematica? Mathematica seems to be able to integrate functions like ...
3
votes
0answers
74 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] ...
3
votes
0answers
58 views

SiegelTheta fails to evaluate when given proper arguments

SiegelTheta often returns error messages when I give it arguments that should be of the correct form. For instance, I have a numerical matrix like ...
3
votes
0answers
164 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
40 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 ...
2
votes
0answers
189 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
113 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
54 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
59 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
vote
0answers
50 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
57 views

MacDonald formula for Modified Bessel Functions

Mathematica seems to not know 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) } K_{i t}(u) ...
1
vote
0answers
205 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 ...
1
vote
0answers
169 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

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
55 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
56 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 ...
0
votes
0answers
55 views

NSolve can not handle a PolyGamma equation

I would like to solve the following equations numerically: ...
0
votes
0answers
84 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
0answers
55 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
62 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
76 views

Real or Imaginary result of spherical bessel and hankel functions of imaginary arguments

I am trying to calculate a rather complicate expression involving Spherical Bessel and Hankel functions. My problem is that somehow for pure imaginary arguments the functions are not pure real or ...
0
votes
0answers
43 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 ...
0
votes
0answers
27 views

Expanding Function Based on Congruence Classes

I'm trying to get Mathematica to expand $\quad \quad \beta_n=\sum_{r=0}^{n}\frac{\alpha_r}{(q^4;q^4)_{n-r}(q^4;q^4)_{n+r}},$ where $\alpha_r$ is an integer valued sequence that is given by ...
0
votes
0answers
46 views

JacobiCN inaccuracy

For my calculations I need JacobiCN[s,k] for s near 0 from the left and k near 1 from the ...
0
votes
0answers
99 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 ...
0
votes
0answers
88 views

Symbolic representation of bessel series derivative

I want to get symbolic expression for BesselJ derivative, where BesselJ is represented as series: ...