solve y=x+x^3/(x^2+1) for x