Last week Bujanda et al published a paper  that gives new expansions for the confluent hypergeometric function. I’ll back up explain what that means before saying more about the new paper.
Hypergeometric functions are something of a “grand unified theory” of special functions. Many functions that come up in application are special cases of hypergeometric function, as shown in the diagram below. (Larger version here.)
I give a brief introduction to hypergeometric functions in these notes (4-page pdf).
Confluent hypergeometric functions
The confluent hypergeometric function corresponds to Hypergeometric 1F1 in the chart above . In Mathematica it goes by
Hypergeometric1F1 and in Python goes by
This function is important in probability because you can use it to compute the CDFs for the normal and gamma distributions.
Bujanda et al give series expansions of the confluent hypergeometric functions in terms of incomplete gamma functions. That may not sound like progress because the incomplete gamma functions are still special functions, but the confluent hypergeometric function is a function of three variables M(a, b; z) whereas the incomplete gamma function γ(s, z) is a function of two variables.
They also give expansions in terms of elementary functions which may be of more interest. The authors give the series expansion below.
The A functions are given recursively by
for n > 1.
The F functions are given by
The authors give approximation error bounds in . In the plot below we show that n = 3 makes a good approximation for M(3, 4.1, z). The error converges to zero uniformly as n increases.
 Blanca Bujanda, José L. López, and Pedro J. Pagola. Convergence expansions of the confluent hypergeometric functions in terms of elementary functions. Mathematics of computation. DOI https://doi.org/10.1090/mcom/3389
 “Confluent” is a slightly archaic name, predating a more systematic naming for hypergeometric functions. The name mFn means the hypergeometric function has m upper parameters and n lower parameters. The confluent hypergeometric function has one of each, so it corresponds to 1F1. For more details, see these notes.