dsolve(diff(y(x),x) = 1/125*(150*x^3+125*x^(1/2)+125+125*y(x)^2-100*x^3*y(x)-500*y(x)*x^(1/2)+20*x^6+200*x^(7/2)+500*x+125*y(x)^3-150*x^3*y(x)^2-750*y(x)^2*x^(1/2)+60*y(x)*x^6+600*y(x)*x^(7/2)+1500*x*y(x)-8*x^9-120*x^(13/2)-600*x^4-1000*x^(3/2))/x,y(x))