>> syms x y(x)
>> equ = diff(y, x) == x^2 + y^2
equ(x) =
D(y)(x) == x^2 + y(x)^2
>> dsolve(equ)
ans =
-((C2*besselj(1/4, x^2/2) + C3*bessely(1/4, x^2/2))/(2*x^(1/2)) + x^(1/2)*(C2*(x*besselj(-3/4, x^2/2) - besselj(1/4, x^2/2)/(2*x)) + C3*(x*bessely(-3/4, x^2/2) - bessely(1/4, x^2/2)/(2*x))))/(x^(1/2)*(C2*besselj(1/4, x^2/2) + C3*bessely(1/4, x^2/2)))