这是个多值函数。
由matlab解得：
((- 4*x^3 + 12*x^2 + 1)^(1/2) - 1)/(2*x^2)
-((- 4*x^3 + 12*x^2 + 1)^(1/2) + 1)/(2*x^2)
代码：
syms x y
solve(x + y + x^2*y^2 - 3,y)
pretty(ans)

然后可用数值积分。MATLAB的符号积分，要花很长时间。