dsolve(diff(y(x),x)+a*y(x)*(y(x)-x)-1 = 0,y(x))