matlab
dsolve('D2y=(1/x)*Dy+x*e^x','x')
ans =
1/2*x^2*C1+(-1+x*log(e))/log(e)^3*e^x+C2