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