dsolve(diff(y(x),x) = 1/25*(30*x^3+25*x^(1/2)+25*y(x)^2-20*x^3*y(x)-100*y(x)*x^(1/2)+4*x^6+40*x^(7/2)+100*x)/x,y(x))