dsolve(2*(x-a1)*(x-a2)*(x-a3)*diff(diff(diff(y(x),x),x),x)+(9*x^2-6*(a1+a2+a3)*x+3*a1*a2+3*a1*a3+3*a2*a3)*diff(diff(y(x),x),x)-2*((n^2+n-3)*x+b)*diff(y(x),x)-n*(n+1)*y(x)=0,y(x))