Some corrections in code: now we work with current values mux, variancex, muy, variancey and updates mux1, variancex1, muy1, variancey1 (instead of storing all successive values in vectors).
Examples added.

Error in title: fixed.