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