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