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