initialization z = (0,...,0)', x1 = ΣQ and for j = 2...Nobs do xj = (0,...,0)'