N1[x_]=x-(x^3-x-1)/(3 x^2-1)
Q1[x_]=Factor[x-N1[N1[N1[x]]]]