y'[x] == y[x]^2 + x^2 - 1, y[0] == 1