dsolve(diff(diff(y(x),x),x) = -((1-al1-bl1)*b1/(b1*x-a1)+(1-al2-bl2)*b2/(b2*x-a2)+(1-al3-bl3)*b3/(b3*x-a3))*diff(y(x),x)-1/(b1*x-a1)/(b2*x-a2)/(b3*x-a3)*(al1*bl1*(a1*b2-a2*b1)*(-a1*b3+a3*b1)/(b1*x-a1)+al2*bl2*(a2*b3-a3*b2)*(a1*b2-a2*b1)/(b2*x-a2)+al3*bl3*(-a1*b3+a3*b1)*(a2*b3-a3*b2)/(b3*x-a3))*y(x),y(x))