x=var('x')
y=function('y',x)
desolve((x^2+1)*diff(y,x)+y^2+1,y)