Sometime ago, we discussed unifying maxima's use of %j and bessel_j (and other Bessel functions) to use bessel_j and friends.
I've taken a look at this finally, and think I know how to do this. Do we want to keep backward compatibility and accept %j? Or just get rid of it completely?
There are a couple of other related issues. The docs for specint show that there's also %p (Legendre functions?) and %m (Whittaker's confluent hypergeometric function?). Do we want to replace these as well. I don't think maxima knows anything about %m. Maybe the specfun package knows about %p. (Sorry for the vagueness. I don't have the docs with me right now.)
Not to mention hstruve[n](x), li[n](x) and so on. These should, perhaps, be changed to hstruve(n,x) and li(n,x). Or something better for li.
Ray