dsolve(diff(diff(y(x),x),x) = -(x^2*((x^2-a1)*(x^2-a2)+(x^2-a2)*(x^2-a3)+(x^2-a3)*(x^2-a1))-(x^2-a1)*(x^2-a2)*(x^2-a3))/x/(x^2-a1)/(x^2-a2)/(x^2-a3)*diff(y(x),x)-(A*x^2+B)/x/(x^2-a1)/(x^2-a2)/(x^2-a3)*y(x),y(x))