Bug 718574 mentions that %j[n](x) is also used in other places to
denote the Bessel function, which we now call bessel_j(n,x).
I think these should be unified. Of course, now that I've looked a
bit more, it would have been a lot easier to modify the Bessel code to
use %j rather than modifying the code to use bessel_j in place of
%j. :-)
Does Macsyma support %j as a Bessel function? Here's a simple test to
tell:
(C7) specint(t^(1/2)*%j[1](2*a^(1/2)*t^(1/2))*%e^(-p*t),t);
- a/p
SQRT(a) %E
(D7) ---------------
2
p
There are also similar issues with other function. psi(x)
(logarithmic derivative of the Gamma function) is the same as
psi[0](x). And psi[n](x) is a known function. Should these be
changed to psi(n,x)? Should psi(x) be removed in favor of psi(0,x)?
What about li[n](x)? (Logarithmic integrals)
The specint package also has lots of other functions like %y[n](x),
%m, %h, etc.
Opinions?
Ray