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