initialization x = [g(a),...,g(1)]', p1 = g(a), and i = 0