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