a(n+1)=((n^2+n+1)a(n)+n^2+3n+1)/(na(n)+1), a(1)=1