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