analysis of anomalies
xai[:,i] = xbi[:,i] - (1/2) * K @ ObsOp(xbi[:,i])