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