dsolve(diff(y(x),x) = -1/32/y(x)*(16*x*y(x)^3-8*y(x)^3-8*y(x)+8*x*y(x)^2-2*x^2*y(x)^3-8+12*x*y(x)-6*x^2*y(x)^2+x^3*y(x)^3)/x,y(x))