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