dsolve(diff(y(x),x) = -1/4*x+1+y(x)^2+7/16*x^2*y(x)-1/2*x*y(x)+5/128*x^4-5/64*x^3+1/16*x^2+y(x)^3+3/8*x^2*y(x)^2-3/4*x*y(x)^2+3/64*y(x)*x^4-3/16*x^3*y(x)+1/512*x^6-3/256*x^5,y(x))