dsolve(diff(diff(y(x),x),x)+diff(y(x),x)*x^(1/2)+(1/4/x^(1/2)+1/4*x-9)*y(x)-x*exp(-1/3*x^(3/2))=0,y(x))