x=var('x')
y=function('y',x)
deq1=diff(y,x)-y^2-1==0
desolve(deq1,[y,x])