dsolve(diff(y(x),x)^2-4*y(x)^3+a*y(x)+b = 0,y(x))