y=function('y',x)
deq=diff(y,x,1)-1-4*y^2==0
desolve(deq,[y,x])