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