Lots of tiny bubbles.... read more >>

Three-term recurrence relations are the basis for computing Bessel functions. Contents A Familiar Three-Term Recurrence Friedrich Bessel Bessel Functions Miller's Algorithm ... read more >>

"ULP" stands for "unit in the last place." An ulps plot samples a fundamental math function such as $\sin{x}$, or a more esoteric function like a Bessel function. The samples are compared with more accurate values obtained from a higher precision computation. A plot of the accuracy, measured in ulps, reveals valuable information about the underlying algorithms.... read more >>

An incredible book, published in several editions from 1909 to 1933, by German mathematicians Eugene Jahnke and Fritz Emde, contains definitions and formulas for mathematical functions, hand-calculated tables of function values, and meticulous hand-drawn 2- and 3-dimensional graphs. An English edition was published by Dover in 1945.... read more >>

I first investigated the tower of powers... read more >>

The Lambert W function deserves to be better known. It pops up in all sorts of places. And our MATLAB function for evaluating the function is a beautiful use of the Halley method.... read more >>