A while back I came upon a seemingly not-too-difficult programming exercise: Define a recurrence  by Compute . This isn't too hard to code up, using perhaps a recursive function to represent . With normal double-precision floats, as increases, the result converges neatly toward 100. Super! Unfortunately, 100 is not even close to the right answer. This … Continue reading Muller's Recurrence - roundoff gone wrong →