x=var('x')
y=function('y',x)
deq=y*diff(y,x,1)+25*x==0
a=desolve(deq,[y,x])
a
diff(a,x)