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