Questions on the special mathematical functions implemented in Mathematica.

learn more… | top users | synonyms

2
votes
1answer
118 views

Numerical evaluation of ChebyshevT

When I evaluate the following Chebyshev series of the first kind in two different ways, I get two very different results: ...
2
votes
2answers
201 views

Define PolyLog so that positive reals evaluate on upper edge of branch cut

Generalizing the previous question: Define Log so that negative reals evaluate on lower edge of branch For real positive values $x>1$, Mathematica's polylogarithm function ...
0
votes
1answer
241 views

How to calculate the unknown quantity in an infinite series?

I'd like to calculate x value in this equation. Basically, I tried to 2 types of method which are FindRoot and NSolve. But, I have failed the calculation caused by these errors up to now. If ...
0
votes
1answer
163 views

Cannot solve equation

I just started using mathematica and I'm facing a problem that I just can't solve. I want to solve the following equation: Code: ...
3
votes
1answer
635 views

HarmonicNumber problem

I am looking for all the roots of HarmonicNumber, in the domain [-30, 1] and [0, 6000] for the real and imaginary parts, respectively, and where the parameter ...
1
vote
2answers
248 views

A regularized hypergeometric function related question

I'm interested in finding a way (if possible) of expressing this specific value of the regularized hypergeometric function in terms of known constants. How might I use Mathematica to check this ...
0
votes
1answer
470 views

How to use `FindRoot` to solve an equation containing a parameter?

I'm trying to derive some of the results of the following paper: Electrodynamics of semiconductor-coated noble metal nanoshells, JT Manassah - Physical Review A In the paper there is matrix $\mathbf ...
4
votes
1answer
143 views

What is the formula Mathematica uses for ZetaZero?

What is the formula/algorithm Mathematica uses for the ZetaZero command?
3
votes
2answers
460 views

Bessel derivatives

I tried to expand BesselJ[k,x] function into a Taylor series with Series command. Here both ...
1
vote
1answer
199 views

How to choose a solution of a differential equation?

I am trying to solve the following differential equation: ...
1
vote
2answers
549 views

Are solid spherical harmonics implemented in Mathematica?

In certain applications, solid spherical harmonics can be very useful. They are essentially the usual, 'surface' spherical harmonics, with the appropriate power of the radius inserted: ...
3
votes
4answers
232 views

Evaluating Erf[x] in arbitrary precision

Is it possible to evaluate Erf[200] with arbitrary precision? I only get 1. as result but I would like to know if a arbitrary ...
1
vote
1answer
311 views

Is there a function for falling factorial in Mathematica

If I had to construct a function for falling factorial in mathematica I'd do something like that (hope I'm not mistaken): fallfact[x_,k_]:=$\prod_{j=0}^{k-1}(x-j)$ But is there a built-in ...
3
votes
1answer
196 views

Symbolically Expanding One-Variable q-hypergeometric Series

In general, how do I get Mathematica to symbolically expand various infinite series to an arbitrary degree of accuracy? I deal with q-hypergeometric series quite frequently, and specifically want to ...
0
votes
2answers
142 views

How can I compute the real part of $\zeta^2$ numerically? [duplicate]

I want to compute and plot $\Re(\zeta(x+iy)^2)$ and $\Im(\zeta(x+iy)^2)$. How can I do that with Mathematica?
3
votes
3answers
176 views

Symbolic derivatives of special functions yield incorrect results

When I evaluate with Mathematica this expression D[ Abs[ Zeta[x + I y]], {x, 2}] + D[ Abs[ Zeta[x + I y]], {y, 2}] 0 the ...
2
votes
1answer
148 views

Equation involving hypergeometric functions

I want to solve this equation but nor NSolve nor Solve are able to do this. ...
1
vote
1answer
190 views

How to calculate the following integral of the multiplication of two Bessel functions?

This integration has an analytical solution and its behavior is described by 1/r^2 function, but Mathematica gives some weird oscillating answer. Can anybody explain this and help me overcome this ...
1
vote
1answer
750 views

Trying to solve a transcendental equation involving Bessel functions

I've never used Mathematica before and am trying to numerically solve equation (12) from this paper. Ideally I'd be able to find the smallest value of $x_{n\nu}$ for $\exp(-kr\pi)$ close to 1, and ...
4
votes
0answers
317 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 ...
4
votes
1answer
109 views

Complex LogIntegral error

Going through Derbyshire's Prime Obsession & trying to take LogIntegral of 20^ZetaZero[1] & comes up with a value of ...
4
votes
1answer
227 views

Problems with series of generalized hypergeometric functions

I have been trying to do the following series: Series[HypergeometricPFQ[{1, 4, 4}, {4 - Sqrt[3], 4 + Sqrt[3]}, z],{z,1,0}] Mathematica says that the result is ...
3
votes
2answers
264 views

How to find the maximum of this function on the positive real line?

I need to maximize this function on the positive real line: $$ \frac{1}{\Gamma(x)^{14}}\cdot\frac{1}{{\frac{323.6}{14x}}^{14x}}\cdot(1.22578*10^{19})^{x-1}e^{-14x} $$ the correct answer should be ...
1
vote
1answer
199 views

I want change gamma function form into beta function form

I want to integrate below form ...
1
vote
1answer
456 views

Error Function Integral (Erf)

Any idea how to solve analytically this integral Integrate[(a Erf[a Sqrt[b/(a^2 + b)] c])/(a^2 + b)^(3/2), a] I tried substitution u=a^2 + b, but it didn't work. ...
3
votes
1answer
444 views

A Bessel & Struve functions related integral

I try to numerically compute this integral and I don't figure out why on earth Mathematica is not able to do it. Is my input correct? Does it possibly have a closed form? ...
2
votes
2answers
302 views

Solving Inequalities with Gamma Function

I am wondering why Mathematica outputs that the following system "cannot be solved with the methods available to Reduce". $\frac{\Gamma(\frac{1}{2}+n)}{n-1}<\frac{\Gamma(\frac{1}{2}+k)}{k-1}$ ...
6
votes
1answer
114 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 ...
1
vote
1answer
318 views

Integral over squared Hermite polynomial

I would like to calculate the uncertainty of the nth Eigenstate of a 1-dim harmonic oscillator. To obtain the result I have to solve the integral $$\int_{-\infty}^{\infty} \psi^* x^2 \psi \:dx$$ ...
4
votes
1answer
292 views

Ordinary differential equation invloving the function composition

I want to solve the following function: DSolve[(A1*Exp[B1*f[x]] + A2*Exp[B2*f[x]])*f'[x] == A1*Exp[B1*x] + A2*Exp[B2*x], f[x], x] And this is what I get as an ...
3
votes
1answer
365 views

Simplify doesn't simplify HypergeometricPFQ with exact arguments

Consider a series: $$\sum_{t=0}^\infty \frac{8^{-11-2t}(22+4t)!}{t!(11+t)!(11+2t)!(32+t)}$$ ...
6
votes
2answers
220 views

Gamma function computation efficiency?

I wonder what kind of algorithm is used to compute the values for the gamma function. Specifically, I am interested in how the computational load increases when the complexity of the input grows. So, ...
1
vote
0answers
157 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
225 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 ...
1
vote
1answer
168 views

Numerical integration of Hankel functions

I would like to know how to perform numeric integration for the following type of integrals in Mathematica. For the following integrand, we can not get the symbolic result. ...
9
votes
2answers
180 views

Integrating a BesselJ integrand to obtain the same result as Maple 16

I would like to check the following integration: Integrate[y*Integrate[1/x*BesselJ[1,x*Exp[I*π/4]]*BesselJ[1,x*Exp[-I*π/4]],{x,0,y}],{y,0,r}] Mathematica 9.0 is ...
0
votes
1answer
478 views

NDSolve + FindRoot for Bessel Zeros

I am trying to use a solution given by Michael E2 in this topic: ...
2
votes
1answer
315 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 ...
6
votes
1answer
334 views

Using FunctionExpand to evaluate symbolic derivatives

Some symbolic derivatives of certain special functions are not expanded automatically, but FunctionExpand often helps to get a derivative-free closed form ...
9
votes
1answer
853 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 ...
3
votes
1answer
619 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 ...
5
votes
1answer
487 views

Strange result for the analytic integration leads to Hypergeometric2F1

The integration result for Integrate[1/(r^2 Sqrt[x/r^(4 - 2 \[Gamma]) + 1]), r] is: ...
4
votes
1answer
307 views

Numerical errors/inaccuracies in ProductLog

Context In cosmology, a fairly accurate model to describe the gravitational potential, $\psi(r)$ of dark matter halos is given by $\psi( r)=\log(1+r)/r$. ...
20
votes
1answer
627 views

Incorrect result from Integrate

Bug introduced in 8.0 and fixed in 10.0 I attempted to calculate the following integral: ...
6
votes
1answer
505 views

How to express an expression with only ArcTan and ArcTanh?

I have an expression which is simply (j/k) x^(j/k) LerchPhi[x,1,j/k)] where 0 < j < k. Manually I have been able ...
-1
votes
1answer
697 views

Gamma Incomplete Function representation in Mathematica and Matlab

I am confused about incomplete gamma function calculation in Mathematica and MATLAB: For example, in Mathematica: Gamma[5,3] = 19.56 But in MATLAB: ...
5
votes
2answers
283 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 ...
3
votes
2answers
452 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, ...
4
votes
2answers
379 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 ...
1
vote
2answers
2k views

Q-Function representation in Mathematica

I observed that there is no Q-function representation in Mathematica. The definition of Q-Function is: \begin{align} Q(x) &= \frac{1}{\sqrt{2\pi}} \int_x^\infty e^{-\frac{u^2}{2}}du \\ ...