solve x * (1. / sqrt(1 + x**2) + C1)**2 - C2=0, x