x_{k+1}=x_k\left(1+\frac{1-ax_k^2}2\right)