[TAForObvReasons]

    mu = mu + (x - mu) / (i+1)

This is more natural when written in C:

    mu += (x - mu) / (i+1)