dsolve(diff(y(x),x) = 1/x^4*(2*x^2-4*x^3*y(x)+1+x^4*y(x)^2+x^6*y(x)^3-3*y(x)^2*x^5+3*y(x)*x^4-x^3),y(x))