Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

18
votes
6answers
5k views

About multi-root search in Mathematica for transcendental equations

I have some questions for multiroot search for transcendental equations. Is there any clever solution to find all the roots for a transcendental equation in a specific range? Perhaps ...
3
votes
1answer
203 views

When does the real part of Zeta vanish on the critical line?

This seems to be a quite a simple problem but I cannot make it work. I am trying to find all values within a given range for which the real part of the Zeta function vanishes on the line: ...
6
votes
2answers
719 views

Why does Expand not work within a function?

I'm writing this fairly simple function: ...
40
votes
2answers
3k views

What is the difference between a few simplification techniques?

I am trying to understand the difference between Refine, Simplify and FullSimplify, and when ...
12
votes
4answers
2k views

Numerical underflow for a scaled error function

I calculate scaled error function defined as f[x_] := Erfc[x]*Exp[x^2] but it can not calculate f[30000.]. ...
2
votes
3answers
207 views

Why do these two different zetas produce the same value?

Zeta[-13] == Zeta[-1] == -1/12 Why do these two different zetas produce the same value?
17
votes
2answers
2k views

Visualizing vector spherical harmonics

I have painstakingly derived the vector-spherical harmonics $\mathbf{V}_{J,\,M}^\ell(\theta, \phi)$, which are the generalization of ordinary spherical harmonics $Y_\ell^m(\theta, \phi)$ to vector ...
8
votes
1answer
1k views

solve an integral equation numerically

I am trying to find a numerical solution for an equation of the form: $$ f(t)=\int_{t_{min}}^{t} {\exp[2 (t^\prime - t)] E(t^\prime) f(t^\prime)} + \int_{t}^{0}{ \exp(t^\prime - t) E(t^\prime) ...
8
votes
3answers
1k views

Checking if the roots of a function are real

I'm trying to determine if the roots of a function are real. How would you do that? (In particular I'm interested in verifying that the roots of LegendreP[6, x] ...
-1
votes
1answer
401 views

Gamma Incomplete Function representation in Mathematica and Matlab

I am confused about Gamma Incomplete function calculation in Mathematica and MatLab: For example, in Mathematica: Gamma[5,3] = 19.56 But in MatLab: ...
15
votes
1answer
412 views

How to represent a continuous monotonic phase of Airy functions?

Note: In this question I am concerned only with real-valued variables and functions. DLMF, §9.8 Airy Functions, Modulus and Phase, formula $9.8.4$ defines the phase of Airy functions: ...
8
votes
2answers
247 views

How to simplify an expression with special functions to zero

The following is a well-known Bessel function identity: $$J_{-n}(z)=(-1)^n J_n(z),\qquad n\in\mathbb Z$$ To check this, I used the following code and the result is as what I expected. ...
6
votes
2answers
287 views

Cannot Get Numerical Results to Match

I try this numerical summation (in two parts) ...
18
votes
1answer
481 views

Incorrect result from Integrate

I attempted to calculate the following integral: ...
10
votes
1answer
706 views

Is my expression too complicated for FullSimplify or am I doing something wrong?

I have a messily defined function $v(h, w)$ with $h, w \in \mathbb{R}$ and with a removable singularity at $h=1/2$, and I am trying to prove some of its properties using Mathematica. In particular I ...
9
votes
4answers
706 views

Finding the roots of Hypergeometric1F1[]

I am trying to find the roots, λ, for this equation: Hypergeometric1F1[1/4 (2 -  λ /β), n + 1, β] for certain ...
2
votes
1answer
428 views

FindRoot giving false roots with Bessel Functions

I have read in some places about the errors associated with FindRoot, but the closest thing I can find on this website seems to be due to the imaginary unit. I am dealing with what should be a ...
0
votes
2answers
544 views

Using Mathematica to find poles of Gamma functions

I am concerned about the expression on the RHS of equation A.5 (page 19) in this paper: $$\int\frac{d^d ...
9
votes
2answers
2k views

problem with coloring spherical harmonics

I want to color a spherical harmonics. So I write as follows. ...
8
votes
1answer
696 views

Hankel Transform integrals won't work in Mathematica

I'm trying to do this integral, which is shown on the Wikipedia page on the Hankel transformation: $$\int_0^{2\pi}\mathrm d\varphi\;e^{\mathrm im\varphi}e^{\mathrm ikr\cos(\varphi)}$$ The answer is ...
4
votes
2answers
298 views

How do I evaluate a symbolic integral involving Hermite polynomials?

I want to test a difficult integral : Integral on all reals of some complicated function involving the Hermitian polynomials, exponentials, squares, factorials, and being general considering any ...
3
votes
2answers
323 views

Assigning an analytical approximation to the error function erf(x)

Working with some iterative integral equations, I have Gaussian density functions involved therein. Integrating the Gaussian function I obtain the error function. When I take the second integration, ...
31
votes
3answers
766 views

Is there a Mathematica API for the functions.wolfram site?

Is there a Mathematica API for the functions.wolfram site? If there's not, has anyone implemented a web scraper for it? For example it would be nice to be able to access ...
18
votes
1answer
803 views

Peirce's quincuncial projection

The Peirce quincuncial projection is the cartographic projection of a sphere onto a square. In short, I would like to see it implemented in Mathematica. Here is my code: ...
8
votes
1answer
177 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 ...
7
votes
1answer
574 views

Visualizing vector-spherical waves

This is a follow-up question to this one on visualizing vector-spherical harmonics. This time, I would like to visualize the vector spherical waves (including the radial dependence). The functions ...
8
votes
1answer
604 views

how to simplify large expression with lots of special functions in it (BesselY, Hypergeometric, MeijerG etc…)

I saw this DE in Maple forum. When solving it using Mathematica 9.01, even though the result was correct (both solutions gave the same numerical answer for some random values), Mathematica's answer ...
6
votes
0answers
202 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 ...
14
votes
1answer
568 views

Incorrect results for elementary integrals when using Integrate

There is a rather simple integral ($K_0$ is the 0-th order MacDonald function) $$\int_0^\infty e^{-x \cosh\xi}\, d\xi = K_0(x)$$ which mathematica cannot solve. This even though the documentation ...
7
votes
2answers
312 views

Find asymptotics of Sum[2^i*Binomial[n-i-1,2*n/3-1],{i,0,n/3}]

I have an expression 2^n / Sum[ 2^i Binomial[ n - i - 1, 2n/3 - 1], { i, 0, n/3}] ...
4
votes
1answer
92 views

Why do certain fractional values in TriangleWave not evaluate?

While answering another question I discovered that TriangleWave does not automatically evaluate when given certain fractional values, specifically fractions with a ...
3
votes
1answer
71 views

How to make the result of Piecewise be a closed interval?

Description This question comes from two questions. Namely Q1 and Q2 The defintion of B-Spline basis function as shown below: Let $\vec{U}=\{u_0,u_1,\ldots,u_m\}$ a nondecreasing sequence of real ...
3
votes
1answer
247 views

Strange behaviour of PolyLog Function

I discovered some strange behaviour of the PolyLog[] Function in Mathematica which seems to me like a bug in the function implementation. It looks like ...
3
votes
2answers
282 views
2
votes
0answers
143 views

Power::infy: error using FindFit [closed]

I have a 2D data to fit (the Airy type pattern from diffraction of Gaussian beam on the circular aperture). Please see the attached picture. I was going to use a function like ...
2
votes
1answer
215 views

Why FullSimplify doesn't work here?

Since the emphasis of this question is on finding a workaround, I decided to post this question with an emphasis on the explanation of the behavior of Mathematica. The Bessel function satisfies the ...
1
vote
1answer
400 views

Plotting Fresnel function

I am trying to plot the partial sums and the cesaro means of the function $\sqrt{|x|}$ and for $a_{n}$, I obtained the following code which contains FresnelS. ...